Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 998,451 T→A 6.4% intergenic (+44/‑39) elfA → / → elfD laminin‑binding fimbrin subunit/putative periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913998,4510TA 6.4% 349.6 / 9.3 124intergenic (+44/‑39)elfA/elfDlaminin‑binding fimbrin subunit/putative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base T (65/50);  new base A (3/5);  total (69/55)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

Only 100 of 124 total aligned reads displayed.

TGGAGATCAGACGTTGACTTTTAAAGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCATGGGCGGCAGGTAAAGGCGGGATTGGACTTGCAGC  >  NC_000913/998305‑998596
                                                                                                                                                  |                                                                                                                                                 
tGGAGATCAGACGTTGACTTTTAAAGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTg                                                                                                                                                <  2:163392/150‑1 (MQ=255)
 ggAGATCAGACGTTGACTTTTAAAGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTgg                                                                                                                                               <  2:1433384/150‑1 (MQ=255)
  gagaTCAGACGTTGACTTTTAAAGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGa                                                                                                                                              >  2:31749/1‑150 (MQ=255)
    gaTCAGACGTTGACTTTTAAAGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCaggagg                                                                                                                                            <  1:947242/150‑1 (MQ=255)
     aTCAGACGTTGACTTTTAAAGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAggg                                                                                                                                           <  1:457119/150‑1 (MQ=255)
      tCAGACGTTGACTTTTAAAGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGAACACTGGAgggg                                                                                                                                          <  1:785828/150‑1 (MQ=255)
       cAGACGTTGACTTTTAAAGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGAACACTGGAggggt                                                                                                                                         <  2:694271/150‑2 (MQ=255)
        aGACGTTGACTTTTAAAGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGat                                                                                                                                        <  1:1242232/150‑1 (MQ=255)
        aGACGTTGACTTTTAAAGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGAACACTGGAggggta                                                                                                                                        <  2:1469287/150‑3 (MQ=255)
         gACGTTGACTTTTAAAGTTGATTTTATCGCTACGGCGACAGCGACTTCATGTAATGTGAATGCGGTAACAAATTTCCATATTAACAATTATTAATAGAACTCATAAATATTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGata                                                                                                                                       <  1:350416/150‑1 (MQ=255)
         gACGTTGACTTTTAAAGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGata                                                                                                                                       >  1:1449772/1‑150 (MQ=255)
         gACGTTGACTTTTAAAGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGata                                                                                                                                       >  1:1398501/1‑150 (MQ=255)
         gACGTTGACTTTTAAAGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGAACACTGGAggggtac                                                                                                                                       <  2:1236506/150‑4 (MQ=255)
           cGTTGACTTTTAAAGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATAtt                                                                                                                                     >  2:942650/1‑150 (MQ=255)
               gACTTTTAAAGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTTATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCata                                                                                                                                 >  2:102870/1‑150 (MQ=255)
                aCTTTTAAAGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCatat                                                                                                                                >  2:425071/1‑150 (MQ=255)
                 cTTTTAAAGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATAtt                                                                                                                               <  1:502374/150‑1 (MQ=255)
                 cTTTTAAAGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATAtt                                                                                                                               >  2:568126/1‑150 (MQ=255)
                  ttttAAAGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATAttt                                                                                                                              <  2:1604789/150‑1 (MQ=255)
                   tttAAAGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATAtttt                                                                                                                             >  1:1073032/1‑150 (MQ=255)
                    ttAAAGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATAttttt                                                                                                                            <  1:987248/150‑1 (MQ=255)
                    ttAAAGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTCCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATAttttt                                                                                                                            >  1:1398700/1‑150 (MQ=255)
                     tAAAGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATAtttttt                                                                                                                           <  2:855652/150‑1 (MQ=255)
                       aaGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTaa                                                                                                                         <  1:1469998/150‑1 (MQ=255)
                        aGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAag                                                                                                                        <  2:1060598/150‑1 (MQ=255)
                        aGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAAGTTCCATATTAACTATTATTAATAGAACTCAGTATTTGGTTTATTAATTAGTACCCCTCCAGTGTGCTGGAGGGGATATTCAGATTTTTTAag                                                                                                                        >  2:1320877/1‑150 (MQ=255)
                           tGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGt                                                                                                                     >  1:1437380/1‑150 (MQ=255)
                           tGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGt                                                                                                                     >  1:1432670/1‑150 (MQ=255)
                                  atCGCTACGGGGAAAGCGACTACAGGGAATGTGAATGCGGTAACAACTTTACATATTAAATATTATAAATACAACACATCAATAGTGTTATTAAAAAGTACCACTACAGAGTTCTGGAGGGGATATTCATATTTTTTAAGACTGAGTAtt                                                                                                              >  1:1309138/1‑150 (MQ=255)
                                      cTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATg                                                                                                          >  2:1450087/1‑150 (MQ=255)
                                      cTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATg                                                                                                          <  2:484616/150‑1 (MQ=255)
                                        aCGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGaa                                                                                                        <  1:1380289/150‑1 (MQ=255)
                                        aCGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGaa                                                                                                        <  1:31749/150‑1 (MQ=255)
                                           gggAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAAc                                                                                                     >  1:1660429/1‑150 (MQ=255)
                                           gggAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAAc                                                                                                     <  1:1030878/150‑1 (MQ=255)
                                           gggAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAAc                                                                                                     >  2:1023409/1‑150 (MQ=255)
                                                aGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCa                                                                                                >  2:1290042/1‑150 (MQ=255)
                                                aGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCa                                                                                                >  1:1690621/1‑150 (MQ=255)
                                                  cGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATa                                                                                              >  2:116478/1‑150 (MQ=255)
                                                   gACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATaa                                                                                             >  1:1267416/1‑150 (MQ=255)
                                                    aCTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAAc                                                                                            >  2:1013802/1‑150 (MQ=255)
                                                     cTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACa                                                                                           <  2:1073032/150‑1 (MQ=255)
                                                     attCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAATAGTGACTATTTATGAAAACTTGCATAACa                                                                                           <  2:1291689/149‑1 (MQ=255)
                                                       tCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACaaa                                                                                         >  1:1452770/1‑150 (MQ=255)
                                                        taGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGATCTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTACTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAg                                                                                        <  2:240636/149‑1 (MQ=255)
                                                        aatgTAATGTTAATGCGGTAACAAATTTCAATATTAACTATTATTAATAGAAATTATTAATTTTTTAATTAATTAGTACCCATCCAGTGTTCTGGAGTTGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAg                                                                                        <  1:292774/147‑1 (MQ=255)
                                                             aaTGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAATAAAGGGAAt                                                                                   <  2:1660429/150‑1 (MQ=255)
                                                               tgtgAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAAttt                                                                                 >  1:1533954/1‑149 (MQ=255)
                                                               tgtgAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATtg                                                                                 <  2:1147215/150‑1 (MQ=255)
                                                               tgtgAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATtg                                                                                 >  2:158644/1‑150 (MQ=255)
                                                                gtgAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATtgt                                                                                >  2:690698/1‑150 (MQ=255)
                                                                  gAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGa                                                                              >  1:577527/1‑150 (MQ=255)
                                                                     tGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCg                                                                           >  2:1421889/1‑150 (MQ=255)
                                                                        ggTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTaa                                                                        >  2:1655922/1‑150 (MQ=255)
                                                                         gTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAg                                                                       <  1:1290042/150‑1 (MQ=255)
                                                                         gTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAg                                                                       >  2:1356275/1‑150 (MQ=255)
                                                                             aaaaTTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTa                                                                   <  1:1672958/149‑1 (MQ=255)
                                                                               aaTTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTTATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAAc                                                                 >  2:182304/1‑150 (MQ=255)
                                                                               aaTTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAAc                                                                 <  1:116478/150‑1 (MQ=255)
                                                                                 tttCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACgg                                                               <  1:1155502/150‑1 (MQ=255)
                                                                                  ttCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGc                                                              <  2:105920/150‑1 (MQ=255)
                                                                                  ttCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGc                                                              >  1:1453810/1‑150 (MQ=255)
                                                                                      atatTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAata                                                          <  1:1655922/150‑1 (MQ=255)
                                                                                      atatTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAata                                                          >  1:490461/1‑150 (MQ=255)
                                                                                       tatTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAatat                                                         >  1:1150323/1‑150 (MQ=255)
                                                                                        atTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATAtt                                                        >  2:472091/1‑150 (MQ=255)
                                                                                         ttAACTATTATTAATAGAACTCATTAATAGTTTTATTAATTAGTACCCCTCCAGTCTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGATTAAATGCAataata                                                       >  2:1290914/1‑150 (MQ=255)
                                                                                           aaCTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACt                                                     >  2:1544113/1‑150 (MQ=255)
                                                                                             cTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCa                                                   <  1:1023409/150‑1 (MQ=255)
                                                                                               attattAATAGAACTCATTAATTGTTTTATTAAATAGTACCCCTCCATTGTTCTGGAGTGTATATTCATATTTTTTAAGATTGACTATTTATGAAAACTTTCATAACAAAGGGAATTGTGAACGAAAGTTTAATGGAAATATTACTCAgc                                                 >  1:314276/1‑150 (MQ=255)
                                                                                                   ttaATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAgctgct                                             <  2:1398501/150‑1 (MQ=255)
                                                                                                   ttaATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAgctgct                                             <  1:1222796/150‑1 (MQ=255)
                                                                                                     aaTAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCa                                           >  2:888490/1‑150 (MQ=255)
                                                                                                        aGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATcc                                        >  1:1381506/1‑150 (MQ=255)
                                                                                                        aGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATcc                                        >  1:1434509/1‑150 (MQ=255)
                                                                                                         gAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCg                                       >  2:294701/1‑150 (MQ=255)
                                                                                                            cTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCAt                                    >  2:60891/1‑150 (MQ=255)
                                                                                                            cTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCAt                                    >  2:1082485/1‑150 (MQ=255)
                                                                                                            cTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCAt                                    <  2:216636/150‑1 (MQ=255)
                                                                                                                 tAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCATGggcg                               >  2:1150041/1‑150 (MQ=255)
                                                                                                                     tGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCATGGGCGGCAg                           >  2:1611929/1‑150 (MQ=255)
                                                                                                                      gTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCATGGGCGGCAgg                          >  2:1512782/1‑150 (MQ=255)
                                                                                                                      gTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCATGGGCGGCAgg                          >  1:1084050/1‑150 (MQ=255)
                                                                                                                             taattaGTACCCCTCCAGTGTTCTGGAGGGGATATTCATAATTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGAGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCATGGGCGGCAGGTTAAGGc                   >  2:74202/1‑150 (MQ=255)
                                                                                                                               attaGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGAGCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCATGGGCGGCAGGTAAAGGCgg                 >  2:1093068/1‑150 (MQ=255)
                                                                                                                                ttaGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCATGGGCGGCAGGTAAAGGCggg                <  2:1150323/150‑1 (MQ=255)
                                                                                                                                    tACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCATGGGCGGCAGGTAAAGGCGGGATTg            <  2:429435/150‑1 (MQ=255)
                                                                                                                                    tACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCATGGGCGGCAGGTAAAGGCGGGATTg            <  2:1306506/150‑1 (MQ=255)
                                                                                                                                    atccccTCCAGAACACTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCATGGGCGGCAGGTAAAGGCGGGATTg            >  1:614511/3‑150 (MQ=255)
                                                                                                                                    atccccTCCAGAACACTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCATGGGCGGCAGGTAAAGGCGGGATTg            >  1:667794/3‑150 (MQ=255)
                                                                                                                                     aCCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCATGGGCGGCAGGTAAAGGCGGGATTgg           >  2:192221/1‑150 (MQ=255)
                                                                                                                                      ccccTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCATGGGCGGCAGGTAAAGGCGGGATTggt          <  1:102870/150‑2 (MQ=255)
                                                                                                                                      ccccTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCATGGGCGGCAGGTAAAGGCGGGATTGGa          <  2:1690621/150‑1 (MQ=255)
                                                                                                                                      ccccTCCAGAACACTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCATGGGCGGCAGGTAAAGGCGGGATTGGa          >  1:900569/1‑150 (MQ=255)
                                                                                                                                       cccTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCATGGGCGGCAGGTAAAGGCGGGATTGGAc         <  2:47831/150‑1 (MQ=255)
                                                                                                                                       cccTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCATGGGCGGCAGGTAAAGGCGGGATTGGAc         <  1:714444/150‑1 (MQ=255)
                                                                                                                                        ccTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCATGGGCGGCAGGTAAAGGCGGGATTGGACt        >  1:791184/1‑150 (MQ=255)
                                                                                                                                          tCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCATGGGCGGCAGGTAAAGGCGGGATTGGACTTg      <  1:158644/150‑1 (MQ=255)
                                                                                                                                           ccAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGAGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCATGGGCGGCAGGTAAAGGCGGGATTGGACTTgc     >  2:525600/1‑150 (MQ=255)
                                                                                                                                              gtgtTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCATGGGCGGCAGGTAAAGGCGGGATTGGACTTgcagc  >  1:739777/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGGAGATCAGACGTTGACTTTTAAAGTTGATTATATCGCTACGGGGAAAGCGACTTCAGGTAATGTGAATGCGGTAACAAATTTCCATATTAACTATTATTAATAGAACTCATTAATTGTTTTATTAATTAGTACCCCTCCAGTGTTCTGGAGGGGATATTCATATTTTTTAAGAGTGACTATTTATGAAAACTTGCATAACAAAGGGAATTGTGACCGTAAGTTTAACGGCAATATTACTCAGCTGCTCATCCGCATGGGCGGCAGGTAAAGGCGGGATTGGACTTGCAGC  >  NC_000913/998305‑998596

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: