Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 819,215 C→T 5.3% N56N (AAC→AAT moaD → molybdopterin synthase, small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913819,2150CT 5.3% 243.5 / 2.3 76N56N (AAC→AATmoaDmolybdopterin synthase, small subunit
Reads supporting (aligned to +/- strand):  ref base C (32/39);  new base T (2/2);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

GCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGCAGCCGTTCAGCGTAGGAGAAGAGTACCCGTGGCTGG  >  NC_000913/819069‑819361
                                                                                                                                                  |                                                                                                                                                  
gCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAAc                                                                                                                                                    >  7:321032/1‑147 (MQ=255)
gCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCa                                                                                                                                                  <  2:233587/149‑1 (MQ=255)
  ccAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGa                                                                                                                                                >  6:264137/1‑149 (MQ=255)
    aGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACg                                                                                                                                              >  5:1744/1‑149 (MQ=255)
      gTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCATCTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCAGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCt                                                                                                                                            >  2:124467/1‑149 (MQ=255)
       tGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCCTCAACCAGACGCTg                                                                                                                                           >  6:163534/1‑149 (MQ=255)
         cgcgAGTTGGTGGGAACAGATGCAACCGAAGGGGCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGt                                                                                                                                         <  6:152814/149‑1 (MQ=255)
              gTTGGTGGGAAGAGATGCAACCGAAGTGGCTGCGGATTTCCCTACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGtt                                                                                                                                    <  7:167084/149‑1 (MQ=255)
              gTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGtt                                                                                                                                    <  8:304138/149‑1 (MQ=255)
              gTTGGTGGGAACAGATGCAACCGAAGCGGCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGTTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGtt                                                                                                                                    <  8:141965/149‑1 (MQ=255)
                  gtggGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTg                                                                                                                                 >  4:227063/1‑148 (MQ=255)
                   tggGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGAc                                                                                                                               <  1:18830/149‑1 (MQ=255)
                    gggAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAg                                                                                                                                      <  3:256902/141‑1 (MQ=255)
                    gggAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAg                                                                                                                                      >  4:256902/1‑141 (MQ=255)
                    gggAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGAc                                                                                                                               >  8:314506/1‑148 (MQ=255)
                      gAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGTCAAATTACTGGCTGCAGTCAACCAGACGCTGGTGAGTTTTGACAAt                                                                                                                            >  8:154913/1‑149 (MQ=255)
                          aGATGCAACCGAAGTGGCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGc                                                                                                                        >  1:267309/1‑149 (MQ=255)
                            aTGCAACCGAAGTGGCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTCACCATCCGct                                                                                                                       >  6:194216/1‑148 (MQ=255)
                              gCAACCGAAGTGGCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGAGAGTTTTGACCATCCGctgac                                                                                                                    >  8:312240/1‑149 (MQ=255)
                                 aCCGAAGTGGCTGCGGATTTCCCAACTGGTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGctgactga                                                                                                                 <  5:141818/149‑1 (MQ=255)
                                  ccGAAGTGGCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCtgactgac                                                                                                                >  7:60325/1‑149 (MQ=255)
                                  ccGAAGTGGCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCtgactgac                                                                                                                >  5:198144/1‑149 (MQ=255)
                                   cGAAGTGGCTGCGGATTTCCCAACTGTTGAAGCGTTACGTCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACg                                                                                                               <  3:263876/149‑1 (MQ=255)
                                    gAAGTGGCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTAGACCATCCGCTGACTGACgg                                                                                                              >  3:117425/1‑149 (MQ=255)
                                       gTGGCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGcga                                                                                                           <  5:264137/149‑1 (MQ=255)
                                       gTGGCTGCGGATTTCCCAACTGTGGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGcga                                                                                                           <  1:188278/149‑1 (MQ=255)
                                          gCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGtt                                                                                                                                    >  4:85925/1‑121 (MQ=255)
                                          gCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGcgacg                                                                                                         <  6:1744/148‑1 (MQ=255)
                                              cGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGt                                                                                                     <  6:198144/148‑1 (MQ=255)
                                              cGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTa                                                                                                    <  7:304826/149‑1 (MQ=255)
                                              cGGAGTTCCCAACTGTTGAAGCGGTACGCGAGCACATGGCTGTGCAGAGCGATCTCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAg                                                                                                      <  5:300178/147‑1 (MQ=255)
                                                     cccAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAATCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTttct                                                                                             >  2:373456/1‑149 (MQ=255)
                                                     cccAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAATCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTttct                                                                                             >  2:368261/1‑149 (MQ=255)
                                                     cccAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTttc                                                                                              <  7:350565/148‑1 (MQ=255)
                                                      ccAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGCCGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTttctt                                                                                            >  7:314589/1‑149 (MQ=255)
                                                      ccAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTttctt                                                                                            >  6:286596/1‑149 (MQ=255)
                                                      ccAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTttctt                                                                                            >  7:291147/1‑149 (MQ=255)
                                                       cAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTttcttc                                                                                           <  6:245916/149‑1 (MQ=255)
                                                        aaCTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTcc                                                                                          >  6:128835/1‑149 (MQ=255)
                                                             ttGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGc                                                                                                  <  2:226533/136‑1 (MQ=255)
                                                             ttGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGc                                                                                                  >  1:226533/1‑136 (MQ=255)
                                                             ttGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCccgccg                                                                                     >  3:93536/1‑149 (MQ=255)
                                                                   cGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTATCTTCCCGCCGGTAAcc                                                                               <  2:43265/149‑1 (MQ=255)
                                                                     ttACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCg                                                                              >  3:342281/1‑148 (MQ=255)
                                                                       aCGCCAGCAGATGGCTTCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACTAAGTAGCTTTCTTCCCGCCGGTAACCGGa                                                                            <  7:154913/148‑1 (MQ=255)
                                                                       aCGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAg                                                                           >  2:98382/1‑149 (MQ=255)
                                                                         gCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGt                                                                         <  8:38114/149‑1 (MQ=255)
                                                                          ccAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGtt                                                                        <  7:312240/149‑1 (MQ=255)
                                                                          ccAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGtt                                                                        <  8:60325/149‑1 (MQ=255)
                                                                          ccAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGtt                                                                        <  4:93536/149‑1 (MQ=255)
                                                                                 aTGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATgg                                                                 <  4:117425/149‑1 (MQ=255)
                                                                                  tGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATgg                                                                 >  5:74406/1‑148 (MQ=255)
                                                                                    gCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCa                                                               <  3:248499/148‑1 (MQ=255)
                                                                                    gCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAg                                                              <  6:74406/149‑1 (MQ=255)
                                                                                      tGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTATTGACAATCCGCTGACTGACGACGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGa                                                             >  8:110763/1‑148 (MQ=255)
                                                                                        cgcAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAAc                                                          <  8:269542/149‑1 (MQ=255)
                                                                                                  tCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAAttgtt                                                <  5:128835/149‑1 (MQ=255)
                                                                                                    gCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAAttgttgt                                              <  1:355463/149‑1 (MQ=255)
                                                                                                    gCTGGGCGCAGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAAttgttg                                               <  1:124467/148‑1 (MQ=255)
                                                                                                        ggcgctggcgctGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCTGTAACCGGAGGTTAAGATGGCAGAAACCAAAAttgttgttg                                            <  5:180533/147‑1 (MQ=255)
                                                                                                          cgctggcgctggAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTcc                                        <  8:321032/149‑1 (MQ=255)
                                                                                                          cgctggcgctggAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTcc                                        <  7:314506/149‑1 (MQ=255)
                                                                                                               gcgctggAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGCa                                     <  6:126628/147‑1 (MQ=255)
                                                                                                                   tggAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAAACGGAGGTTAAGATGTCAGAAACCAAAATTGTTGTTGGTCCGCAgaagtt                               >  1:87591/1‑144 (MQ=255)
                                                                                                                     gAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGCAGCCGTTCa                             <  2:267309/149‑1 (MQ=255)
                                                                                                                        gATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGCAGCCGTTCAGCg                          <  5:286596/149‑1 (MQ=255)
                                                                                                                          tGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGCAGCCGTTCAGCg                          >  5:16470/1‑147 (MQ=255)
                                                                                                                           ggCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGCAGCCGTTCAGCGTa                        <  4:342281/148‑1 (MQ=255)
                                                                                                                           ggCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGCAGCCGTTCAGCGGAg                       >  3:257436/1‑149 (MQ=255)
                                                                                                                           ggCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGACAGAAACCAAAATTGTTGTTGGTCCGCAGCCGTTCAGCGTAg                       >  7:257911/1‑149 (MQ=255)
                                                                                                                              aaaTTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGCAGCCGTTCAGCGTAg                       >  8:25233/1‑146 (MQ=255)
                                                                                                                              aaaTTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACAATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACAGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGCAGCCGTTCAGCGTAg                       >  2:26875/1‑146 (MQ=255)
                                                                                                                                   aCTGGCTGCCGTCAATCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGCAGCCGTTCAGCGTAGGagaaga                <  1:373456/148‑1 (MQ=255)
                                                                                                                                   aCTGGCTGCCGTCAATCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGCAGCCGTTCAGCGTAGGagaaga                <  1:368261/148‑1 (MQ=255)
                                                                                                                                   aCTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGACAGAAACCAAAATTGTTGTTGGTCCGCAGCCGTTCAGCGTAGGagaaga                <  8:257911/148‑1 (MQ=255)
                                                                                                                                     tGGCTGCCGTCAAACAGACGCTGGTGAGATTTGACCATCCGCTGACTGACGGCGACTAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGCCGCCGTTCAGCGTATGAGAAGCGTa             >  8:314484/1‑149 (MQ=255)
                                                                                                                                                aaCCAGACGCTGGTGAGTTTTGACTATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGCAGCCGTTCAGCGTAGGAGAAGAGTACCCGtggctgg  >  8:262862/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCCAACTGTTGAAGCGTTACGCCAGCACATGGCTGCGCAGAGCGATCGCTGGGCGCTGGCGCTGGAAGATGGCAAATTACTGGCTGCCGTCAACCAGACGCTGGTGAGTTTTGACCATCCGCTGACTGACGGCGACGAAGTAGCTTTCTTCCCGCCGGTAACCGGAGGTTAAGATGGCAGAAACCAAAATTGTTGTTGGTCCGCAGCCGTTCAGCGTAGGAGAAGAGTACCCGTGGCTGG  >  NC_000913/819069‑819361

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: