Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,807,318 A→G 5.1% intergenic (+19/‑82) pfkB → / → ydiZ 6‑phosphofructokinase II/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,807,3180AG 5.1% 456.2 / 8.0 135intergenic (+19/‑82)pfkB/ydiZ6‑phosphofructokinase II/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (72/56);  new base G (4/3);  total (76/59)
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.97e-01

Only 100 of 139 total aligned reads displayed.

AGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCATGTTGCATGAGGATACATTGACTGAAA  >  NC_000913/1807170‑1807466
                                                                                                                                                    |                                                                                                                                                    
 gAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCa                                                                                                                                                      >  1:244517/1‑148 (MQ=255)
 gAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCa                                                                                                                                                      >  2:599125/1‑148 (MQ=255)
     aTGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTgg                                                                                                                                                  <  1:599125/148‑1 (MQ=255)
          tctctTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGTAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAat                                                                                                                                            >  4:824316/1‑149 (MQ=255)
          tctctTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGaa                                                                                                                                             >  2:339388/1‑148 (MQ=255)
          tctctTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAATGCTGGGGGaa                                                                                                                                             >  1:854656/1‑148 (MQ=255)
           ctctTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAat                                                                                                                                            >  1:495331/1‑148 (MQ=255)
            tctTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAatca                                                                                                                                          >  3:739029/1‑149 (MQ=255)
             ctTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAATGCTGGGGGAatgt                                                                                                                                          <  2:243956/148‑3 (MQ=255)
               tGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAAtcatca                                                                                                                                       <  4:687347/149‑1 (MQ=255)
                gAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAATGCTGGGGGAatgtttt                                                                                                                                       <  2:345823/148‑6 (MQ=255)
                 aaGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCAc                                                                                                                                      >  1:120842/1‑148 (MQ=255)
                  agagATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCa                                                                                                                                    >  4:604208/1‑149 (MQ=255)
                    agaTGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAAc                                                                                                                                  <  3:862616/149‑1 (MQ=255)
                     gaTGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAAc                                                                                                                                  >  2:484811/1‑148 (MQ=255)
                      aTGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAAcc                                                                                                                                 >  1:750325/1‑148 (MQ=255)
                      aTGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCt                                                                                                                                >  3:655359/1‑149 (MQ=255)
                          ttCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTggg                                                                                                                                                 <  4:426807/128‑1 (MQ=255)
                          ttCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTggg                                                                                                                                                 >  3:426807/1‑128 (MQ=255)
                          ttCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTg                                                                                                                                                   <  2:543955/126‑1 (MQ=255)
                          ttCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTg                                                                                                                                                   >  1:543955/1‑126 (MQ=255)
                          ttCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCg                                                                                                                            >  4:578222/1‑149 (MQ=255)
                          ttCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCg                                                                                                                            >  3:569014/1‑149 (MQ=255)
                              ttttGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCa                                                                                                                         >  2:230616/1‑148 (MQ=255)
                                   gCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGt                                                                                                                   <  3:578222/149‑1 (MQ=255)
                                    cGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGt                                                                                                                   >  1:126439/1‑148 (MQ=255)
                                                   tGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGCAACGCGTTTCTCCGACTATGCTCa                                                                                                  >  4:428930/1‑149 (MQ=255)
                                                    gCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCa                                                                                                  >  4:634407/1‑149 (MQ=255)
                                                       gCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCaaa                                                                                                <  2:126439/148‑1 (MQ=255)
                                                          acacTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTc                                                                                            >  4:352872/1‑149 (MQ=255)
                                                             cTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCAtg                                                                                         >  3:728413/1‑149 (MQ=255)
                                                             cTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCAt                                                                                          >  1:158876/1‑148 (MQ=255)
                                                              tCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCAtgt                                                                                        >  4:178092/1‑149 (MQ=255)
                                                              tCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCAtg                                                                                         <  1:694939/148‑1 (MQ=255)
                                                               cAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTgg                                                                                                                                                  <  2:373549/90‑1 (MQ=255)
                                                               cAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTgg                                                                                                                                                  >  1:373549/1‑90 (MQ=255)
                                                               cAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCAtg                                                                                         >  1:681962/1‑147 (MQ=255)
                                                               cAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCAtg                                                                                         <  2:681962/147‑1 (MQ=255)
                                                                 aTCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGa                                                                                      >  1:3461/1‑148 (MQ=255)
                                                                  tCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCaaa                                                                                                <  1:695003/137‑1 (MQ=255)
                                                                  tCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCaaa                                                                                                >  2:695003/1‑137 (MQ=255)
                                                                  tCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGAt                                                                                     <  1:339388/148‑1 (MQ=255)
                                                                   cAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATa                                                                                    >  2:36995/1‑148 (MQ=255)
                                                                     gggAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAATGCTGGGGGAat                                                                                                                                            >  4:272999/1‑90 (MQ=255)
                                                                      ggAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAatcatc                                                                                                                                        <  4:507507/93‑1 (MQ=255)
                                                                       gAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACaa                                                                                >  1:269900/1‑148 (MQ=255)
                                                                        aaCACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACaaa                                                                               <  3:26857/148‑1 (MQ=255)
                                                                        aaCACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACaaa                                                                               >  4:26857/1‑148 (MQ=255)
                                                                        aaCACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACaaa                                                                               >  1:639325/1‑148 (MQ=255)
                                                                        aaCACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAg                                                                              <  3:152114/149‑1 (MQ=255)
                                                                        aaCACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAg                                                                              <  3:604208/149‑1 (MQ=255)
                                                                         acacGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAgg                                                                             <  4:655359/149‑1 (MQ=255)
                                                                               cTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGa                                                                        >  1:94214/1‑148 (MQ=255)
                                                                                   gCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTa                                                                    <  1:36995/148‑1 (MQ=255)
                                                                                     tCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATgg                                                                 >  4:656601/1‑149 (MQ=255)
                                                                                     tCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATg                                                                  <  3:773815/148‑1 (MQ=255)
                                                                                     tCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATg                                                                  >  4:773815/1‑148 (MQ=255)
                                                                                      cccATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCaaa                                                                                                <  2:737202/117‑1 (MQ=255)
                                                                                      cccATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCaaa                                                                                                >  1:737202/1‑117 (MQ=255)
                                                                                      cccATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGc                                                                >  3:565272/1‑149 (MQ=255)
                                                                                           gACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTg                                                           >  3:156530/1‑149 (MQ=255)
                                                                                                tACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGa                                                       <  1:33031/148‑1 (MQ=255)
                                                                                                 aCGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGAt                                                      <  2:120842/148‑1 (MQ=255)
                                                                                                  cGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATc                                                     >  2:571782/1‑148 (MQ=255)
                                                                                                   gCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCtt                                                   >  3:568742/1‑149 (MQ=255)
                                                                                                      aaaaaTTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTc                                                >  3:851128/1‑149 (MQ=255)
                                                                                                      aaaaaTTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTc                                                >  3:275398/1‑149 (MQ=255)
                                                                                                              aCGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATg                                         <  2:3461/148‑1 (MQ=255)
                                                                                                                gCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCa                                      >  3:556304/1‑149 (MQ=255)
                                                                                                                   tACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATa                                    >  1:745220/1‑148 (MQ=255)
                                                                                                                     ccTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCa                                                                                           <  2:4255/91‑1 (MQ=255)
                                                                                                                     ccTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCa                                                                                           >  1:4255/1‑91 (MQ=255)
                                                                                                                     ccTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAAcc                                 <  3:489522/149‑1 (MQ=255)
                                                                                                                      cTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCAtgt                                                                                        >  4:514307/1‑93 (MQ=255)
                                                                                                                      cTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCAtgt                                                                                        <  3:514307/93‑1 (MQ=255)
                                                                                                                         tCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCa                             >  4:220280/1‑149 (MQ=255)
                                                                                                                             gCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCAt                            >  3:338035/1‑146 (MQ=255)
                                                                                                                              cTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCATGtt                         >  1:193320/1‑148 (MQ=255)
                                                                                                                                   aaaaaCATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCATGTTGCATg                    <  1:519896/148‑1 (MQ=255)
                                                                                                                                   aaaaaCATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCATGTTGCATg                    <  2:94214/148‑1 (MQ=255)
                                                                                                                                   aaaaaCATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCATGTTGCATGa                   >  3:32595/1‑149 (MQ=255)
                                                                                                                                   aaaaaCATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCATGTTGCATGa                   <  4:851128/149‑1 (MQ=255)
                                                                                                                                   aaaaaCATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCATGTTACATGa                   >  3:576703/1‑149 (MQ=255)
                                                                                                                                    aaaaCATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCATGTTGCATGa                   >  2:625639/1‑148 (MQ=255)
                                                                                                                                    aaaaCATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCATGTTGCATGa                   >  2:92126/1‑148 (MQ=255)
                                                                                                                                     aaaCATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCATGTTGCATGAgg                 <  3:178092/149‑1 (MQ=255)
                                                                                                                                     aaaCATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCATGTTGCATGAg                  >  1:424741/1‑148 (MQ=255)
                                                                                                                                     aaaCATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAACTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCATGTTGCATGAgg                 <  3:614397/149‑1 (MQ=255)
                                                                                                                                        cATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCATGTTGCATGAGGATa              >  4:434027/1‑149 (MQ=255)
                                                                                                                                        cATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCATGTTGCATGAGGATa              <  3:656601/149‑1 (MQ=255)
                                                                                                                                            cccccAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCATGTTGCATGAGGATACAtt          <  4:275398/149‑1 (MQ=255)
                                                                                                                                              cgcAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCATGTTGCATGAGGATACATTGa        <  4:382257/147‑1 (MQ=255)
                                                                                                                                              cccAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCATGTTGCATGAGGATACATTg         <  2:558055/148‑1 (MQ=255)
                                                                                                                                              cccAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCATGTTGCATGAGGATACATTGa        <  4:160949/149‑1 (MQ=255)
                                                                                                                                               ccAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCATGTTGCATGAGGATACATTGa        >  2:454952/1‑148 (MQ=255)
                                                                                                                                               ccAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCATGTTGCATGAGGATACATTGAc       <  3:253386/149‑1 (MQ=255)
                                                                                                                                               ccAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCATGTTGCATGAGGATACATTGAc       >  3:6033/1‑149 (MQ=255)
                                                                                                                                               ccAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCATGTTGCATGAGGATACATTGAc       >  3:549249/1‑149 (MQ=255)
                                                                                                                                                cAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATa                                                                                    <  1:866455/71‑1 (MQ=255)
                                                                                                                                                    aTTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCATGTTGCATGAGGATACATTGACTGaaa  <  4:739029/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGAAAATGCCTCTCTTGAAGAGATGGTTCGTTTTGGCGTAGCTGCGGGGAGTGCAGCCACACTCAATCAGGGAACACGTCTGTGCTCCCATGACGATACGCAAAAAATTTACGCTTACCTTTCCCGCTAACAAAAACATTCCCCCAGCATTGGGGGAATCATCACCAACCTGTCGGCAACGCGTTTCTCCGACTATGCTCAAAAGTCATGTGATAACAAAGGGGTGAACTATGGCCAGTGGCGATCTTGTCCGTTATGTCATAACCGTCATGTTGCATGAGGATACATTGACTGAAA  >  NC_000913/1807170‑1807466

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: