Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,291,373 G→A L11L (TTG→TTA yraK → putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,291,3730GA98.2% 379.7 / ‑1.1 108L11L (TTG→TTAyraKputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (56/50);  total (58/50)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

Only 100 of 112 total aligned reads displayed.

TGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAGGACGGGAGTTTACCTCTGCAAC  >  NC_000913/3291225‑3291519
                                                                                                                                                    |                                                                                                                                                  
tGTACGCTGGGGAAATCAGGCAGATCAACAATGTAAGTTTGCTTTCAGTACACCGGATTCAGAACCAACAACCTCTGTATTAAAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTttatt                                                                                                                                                  >  1:860879/1‑151 (MQ=255)
 gTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTttattg                                                                                                                                                 <  2:269050/151‑1 (MQ=255)
 gTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTttatt                                                                                                                                                  <  1:116282/150‑1 (MQ=255)
 gTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTttatt                                                                                                                                                  <  1:1358613/150‑1 (MQ=255)
   aCGCTGGGGAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGa                                                                                                                                                <  2:439272/150‑1 (MQ=255)
    cGCTGGGGAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGat                                                                                                                                               >  2:764184/1‑150 (MQ=255)
    cGCTGGGGAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGat                                                                                                                                               <  2:39739/150‑1 (MQ=255)
     gCTGGGGAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGatat                                                                                                                                             >  2:193452/1‑151 (MQ=255)
      cTGGGGAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATc                                                                                                                                            >  2:917618/1‑151 (MQ=255)
      cTGGGGAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATc                                                                                                                                            >  1:442304/1‑151 (MQ=255)
       tGGGGAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATcc                                                                                                                                           >  2:700193/1‑151 (MQ=255)
         gggAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCcac                                                                                                                                         >  2:1208366/1‑151 (MQ=255)
          ggAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCcaca                                                                                                                                        <  2:1166631/151‑1 (MQ=255)
          ggAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCcaca                                                                                                                                        <  1:27847/151‑1 (MQ=255)
            aaaTCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATc                                                                                                                                      <  2:644683/151‑1 (MQ=255)
                cAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTgc                                                                                                                                  <  2:1274520/150‑1 (MQ=255)
                  ggCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTgcg                                                                                                                                 <  1:370450/150‑1 (MQ=255)
                    cAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCtt                                                                                                                              >  2:183035/1‑151 (MQ=255)
                    cAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCtt                                                                                                                              <  2:655186/151‑1 (MQ=255)
                       aTCAACAATGGCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAACGCCTCTTATAACAGGACTTTTATTGATATCCACATCCTGCGCTTAt                                                                                                                            >  1:6164/1‑150 (MQ=255)
                        tCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGc                                                                                                                          >  1:682524/1‑151 (MQ=255)
                            caaTGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGcctcc                                                                                                                      >  1:509779/1‑151 (MQ=255)
                             aaTGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGcctcct                                                                                                                     <  2:1349622/151‑1 (MQ=255)
                              aTGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCCACATCCTGCGCTTATGCctcctc                                                                                                                    >  2:145190/1‑151 (MQ=255)
                              aTGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCctcctc                                                                                                                    >  2:829207/1‑151 (MQ=255)
                                 tCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCagg                                                                                                                 >  2:148721/1‑151 (MQ=255)
                                  cAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCagg                                                                                                                 >  1:68237/1‑150 (MQ=255)
                                   aGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCaggag                                                                                                               <  2:75094/151‑1 (MQ=255)
                                    gcttGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGCCATAGCGCAGTGCCAATAAGGATAAAAAAATGAAAAGTGCGCCTCTAATAACAGGTCTTTTACTGATATCGACATCCTGCGTTTTTGCCTCCTCaggtgg                                                                                                              >  1:733750/3‑151 (MQ=255)
                                       tGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGgtg                                                                                                           >  2:360093/1‑151 (MQ=255)
                                       tGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGgt                                                                                                            <  2:1318332/150‑1 (MQ=255)
                                       tGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGgt                                                                                                            <  1:429364/150‑1 (MQ=255)
                                        gCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGgtgt                                                                                                          <  2:119371/151‑1 (MQ=255)
                                         cTTTCAGTACACCGGATTCAGAACATACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGtgtg                                                                                                         >  1:989908/1‑151 (MQ=255)
                                              aGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGcc                                                                                                    >  1:267856/1‑151 (MQ=255)
                                                  cacCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACa                                                                                                <  2:122054/151‑1 (MQ=255)
                                                   acCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGAc                                                                                                 >  1:246516/1‑149 (MQ=255)
                                                     cGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGc                                                                                              >  2:467379/1‑150 (MQ=255)
                                                     cGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGc                                                                                              >  2:647825/1‑150 (MQ=255)
                                                      ggATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCAc                                                                                            >  2:75635/1‑151 (MQ=255)
                                                       gATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACt                                                                                           <  2:209160/151‑1 (MQ=255)
                                                       gATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCACCTCAGGAGGGTGTGGAGCCGACAGCACt                                                                                           >  2:883294/1‑151 (MQ=255)
                                                            aGAACCTACAACCTCTGTATTACAAGGCAAAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGATGGTGTGGATCAGACAGCACTAGc                                                                                        <  1:4416/149‑1 (MQ=255)
                                                             gAACCCAAAACCCCTGGAATACCAGGGACCGCCCAGGTCCCTTTAAGAAAAAAAAAAGAAAAGGGCGCCCCTTTAAACCGGGCGTTTATTGATACCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCg                                                                                       <  2:906943/149‑1 (MQ=255)
                                                              aaCCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTg                                                                                    <  1:614637/151‑1 (MQ=255)
                                                              aaCCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTg                                                                                    >  2:1145674/1‑151 (MQ=255)
                                                                 cTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCg                                                                                  >  1:434980/1‑150 (MQ=255)
                                                                   aCAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCggc                                                                                <  1:682116/150‑3 (MQ=255)
                                                                   aCAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGAc                                                                                <  2:56077/150‑1 (MQ=255)
                                                                    cAACCTCTTTATCACAAGGCACAGCGCAGTGCCATTATGGATAAAAAAATGAAAAGAGCGCCTCTTATAACACGTCTATTATTGATATCCACATCCTGCCCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACaa                                                                              <  1:88601/151‑1 (MQ=255)
                                                                    cAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGa                                                                                 <  1:772178/148‑1 (MQ=255)
                                                                      aCCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAAt                                                                            >  1:410950/1‑151 (MQ=255)
                                                                       ccTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAAtt                                                                           >  2:642140/1‑151 (MQ=255)
                                                                        ctctGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAAtt                                                                           >  1:1346624/1‑150 (MQ=255)
                                                                         tctGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTAc                                                                         >  2:228277/1‑151 (MQ=255)
                                                                           tGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTAc                                                                         >  1:475135/1‑149 (MQ=255)
                                                                            gTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTAcagc                                                                      >  1:622176/1‑151 (MQ=255)
                                                                            gTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTAca                                                                        <  2:755911/149‑1 (MQ=255)
                                                                            gTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATAACCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTAcag                                                                       <  1:1376560/150‑1 (MQ=255)
                                                                              aTTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTAcagcag                                                                    <  1:208813/151‑1 (MQ=255)
                                                                                 aCAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAgtg                                                                  >  1:510161/1‑150 (MQ=255)
                                                                                   aaGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTgg                                                               >  2:522521/1‑151 (MQ=255)
                                                                                    aGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGt                                                              <  2:949578/151‑1 (MQ=255)
                                                                                       cacaGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTtga                                                           <  2:1061730/151‑1 (MQ=255)
                                                                                        acaGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTtgat                                                          <  2:673159/151‑1 (MQ=255)
                                                                                        acaGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTtgat                                                          <  2:1220901/151‑1 (MQ=255)
                                                                                         caGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTtgatg                                                         >  1:815822/1‑151 (MQ=255)
                                                                                         caGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTAGAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTtgatg                                                         >  2:1422419/1‑151 (MQ=255)
                                                                                             gcAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGt                                                     >  2:918844/1‑151 (MQ=255)
                                                                                               aGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGt                                                     <  2:898263/149‑1 (MQ=255)
                                                                                                       tAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAAc                                           >  2:1073278/1‑151 (MQ=255)
                                                                                                        aaGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAAcc                                          >  2:747295/1‑151 (MQ=255)
                                                                                                          ggATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCGCTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACcag                                        <  2:131010/151‑1 (MQ=255)
                                                                                                          ggATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACcag                                        <  2:1024144/151‑1 (MQ=255)
                                                                                                           gATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAAcc                                          >  2:188873/1‑148 (MQ=255)
                                                                                                              aaaaaaaTGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGCGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACcagac                                      >  1:263400/1‑149 (MQ=255)
                                                                                                              aaaaaaaTGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACcagaca                                     <  2:426225/150‑1 (MQ=255)
                                                                                                              aaaaaaaTGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACcagac                                      >  2:732017/1‑149 (MQ=255)
                                                                                                                 aaaaTGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGAt                                  >  1:1249925/1‑150 (MQ=255)
                                                                                                                 aaaaTGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATa                                 <  2:385066/151‑1 (MQ=255)
                                                                                                                  aaaTGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATaa                                <  1:1384415/151‑1 (MQ=255)
                                                                                                                  aaaTGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATaa                                <  1:498226/151‑1 (MQ=255)
                                                                                                                    aTGAAAAGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACg                              >  1:721026/1‑151 (MQ=255)
                                                                                                                      gAAAAGAGCGCATCTTATAACAGCTCTTTTATTGATATCAACATCCTGCGCCTACGCCACCCCACGAGGGGGTGGAGCCGACATCATTAGCGGTGCGACAAATTACAGCAGTGTGGTGGATGGTGTTACGCGGAACGAGACAGATAAGGTg                            >  1:105694/1‑151 (MQ=255)
                                                                                                                        aaaGAGCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGAc                          <  2:1279466/151‑1 (MQ=255)
                                                                                                                          agagCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAg                        <  2:331880/151‑1 (MQ=255)
                                                                                                                           gagCGCCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAgg                       <  2:62917/151‑1 (MQ=255)
                                                                                                                              cgcCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAGGAc                     >  1:706006/1‑150 (MQ=255)
                                                                                                                               gcCTCTTATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAGGACgg                   >  1:232636/1‑151 (MQ=255)
                                                                                                                                    ttATAGCAGGTCTTTTACTGATATCAACATCCTGCGCTTATGCGTCCTCAGGGGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATAACAGCAGGGTGGTTGATGAAGTTACGGTGAACAAGACAGATAACGTGACAGGACGGGAGtt              >  1:967190/1‑151 (MQ=255)
                                                                                                                                    ttATAACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAGGACGGGAGtt              >  2:1041811/1‑151 (MQ=255)
                                                                                                                                       taACAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAGGACGGGAGTTTAc           >  2:955850/1‑151 (MQ=255)
                                                                                                                                        aaCAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAGGACGGGAGTTTAcc          >  1:73650/1‑151 (MQ=255)
                                                                                                                                        aaCAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAGGACGGGAGTTTAcc          >  1:1297217/1‑151 (MQ=255)
                                                                                                                                          cAGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAGGACGGGAGTTTACctc        >  1:1349424/1‑151 (MQ=255)
                                                                                                                                           aGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAGGACGGGAGTTTACctct       <  1:360080/151‑1 (MQ=255)
                                                                                                                                           aGGTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAGGACGGGAGTTTACctct       <  1:259866/151‑1 (MQ=255)
                                                                                                                                            ggTCTTTTATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAGGACGGGAGTTTACCTCTg      <  1:1091524/151‑1 (MQ=255)
                                                                                                                                                ttttATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAGGACGGGAGTTTACCTCTGCAAc  >  2:1383180/1‑151 (MQ=255)
                                                                                                                                                 tttATTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAGGACGGGAGTTTACCTCTGCAAc  <  2:461725/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGTACGCTGGGGAAATCAGGCAGATCAACAATGTCAGTTTGCTTTCAGTACACCGGATTCAGAACCTACAACCTCTGTATTACAAGGCACAGCGCAGTGCCATTAAGGATAAAAAAATGAAAAGAGCGCCTCTTATAACAGGTCTTTTGTTGATATCCACATCCTGCGCTTATGCCTCCTCAGGAGGGTGTGGAGCCGACAGCACTAGCGGTGCGACAAATTACAGCAGTGTGGTTGATGATGTTACGGTGAACCAGACAGATAACGTGACAGGACGGGAGTTTACCTCTGCAAC  >  NC_000913/3291225‑3291519

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