Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,587,670 A→T 6.4% V22D (GTT→GAT) ‡ ydeQ ← putative fimbrial adhesin protein YdeQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,587,6700AT 6.4% 290.8 / 7.7 94V22D (GTT→GAT) ‡ydeQputative fimbrial adhesin protein YdeQ
Reads supporting (aligned to +/- strand):  ref base A (41/47);  new base T (4/2);  total (45/49)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01

CGTAATCATTCCAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACAT  >  NC_000913/1587524‑1587815
                                                                                                                                                  |                                                                                                                                                 
cgTAATCATTCCAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACtt                                                                                                                                                >  1:1001442/1‑150 (MQ=255)
 gTAATCATTCCAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACttt                                                                                                                                               <  2:474342/150‑1 (MQ=255)
      cATTCCAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTtgcctg                                                                                                                                          >  1:403170/1‑150 (MQ=255)
      cATTCCAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTtgcctg                                                                                                                                          >  2:1224896/1‑150 (MQ=255)
       aTTCCAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAGTTTTtgccttc                                                                                                                                         <  2:145271/150‑3 (MQ=255)
          ccAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCAt                                                                                                                                      >  1:1943194/1‑150 (MQ=255)
           cAGCAACTTATATTCTGAGACAAGTCTACAACCATATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCAtt                                                                                                                                     <  2:1295919/150‑1 (MQ=255)
           cAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATAGACATTACAACTGAAGGCAAAAACTTTGCCTGCCAtt                                                                                                                                     <  2:1040327/150‑1 (MQ=255)
            aGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTa                                                                                                                                    >  2:546071/1‑150 (MQ=255)
              cAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACGTAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAAt                                                                                                                                  <  2:1594819/150‑1 (MQ=255)
               aaCTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATa                                                                                                                                 <  2:490877/150‑1 (MQ=255)
                  ttATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATagaa                                                                                                                              <  2:115964/150‑1 (MQ=255)
                  ttATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATagaa                                                                                                                              <  1:639480/150‑1 (MQ=255)
                   tataTGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATagaag                                                                                                                             <  2:481806/150‑1 (MQ=255)
                    atatGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATagaaga                                                                                                                            <  1:1813095/150‑1 (MQ=255)
                          tGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATa                                                                                                                      >  2:189415/1‑150 (MQ=255)
                           gagaCAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATAc                                                                                                                     <  2:1943194/150‑1 (MQ=255)
                            tgaaaaGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATAcc                                                                                                                    <  1:1929258/146‑1 (MQ=255)
                            agaCAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATAcc                                                                                                                    <  1:1767837/150‑1 (MQ=255)
                               cAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGaa                                                                                                                 >  1:1501145/1‑150 (MQ=255)
                               cAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGaa                                                                                                                 >  2:285310/1‑150 (MQ=255)
                                   tCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAgg                                                                                                             <  2:1464804/150‑1 (MQ=255)
                                   tCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAgg                                                                                                             <  2:1346791/150‑1 (MQ=255)
                                      aCAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATATAAGATAAATACCGAACAGGAcc                                                                                                          >  2:610734/1‑150 (MQ=255)
                                       cAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCt                                                                                                         >  2:1068408/1‑150 (MQ=255)
                                       cAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCt                                                                                                         >  2:435441/1‑150 (MQ=255)
                                         aCCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCttt                                                                                                       <  2:595008/150‑1 (MQ=255)
                                         aCCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACAGAACAGGACCttt                                                                                                       >  2:1097758/1‑150 (MQ=255)
                                          ccAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTta                                                                                                      <  2:109001/150‑1 (MQ=255)
                                           cAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTtat                                                                                                     <  1:157271/150‑1 (MQ=255)
                                           cAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTtat                                                                                                     >  1:780859/1‑150 (MQ=255)
                                             gATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATAc                                                                                                   <  2:1001442/150‑1 (MQ=255)
                                             gATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATAc                                                                                                   <  1:1096747/150‑1 (MQ=255)
                                              aTTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACt                                                                                                  <  1:1261725/150‑1 (MQ=255)
                                                  tGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGAtt                                                                                              >  1:47391/1‑150 (MQ=255)
                                                  tGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGAtt                                                                                              <  1:537513/150‑1 (MQ=255)
                                                   ggCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTg                                                                                             >  2:1345822/1‑150 (MQ=255)
                                                      ccGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGttt                                                                                          >  1:1439625/1‑150 (MQ=255)
                                                        ggCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTa                                                                                        >  1:947094/1‑150 (MQ=255)
                                                         catGTATTACAGTGTCTAGATTAACATAAACATAGGTAGTACCTGCTCCAATACTTGACCATCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGACATTAATATAAGATAAATACCGAACATGACCTTTATACTGATTGTTTTAc                                                                                       <  2:500648/148‑1 (MQ=255)
                                                          cTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTAcc                                                                                      <  1:1235456/150‑1 (MQ=255)
                                                          cTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAACGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTAcc                                                                                      <  1:482349/150‑1 (MQ=255)
                                                             tATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCAt                                                                                   >  2:1363756/1‑150 (MQ=255)
                                                              aTTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATg                                                                                  >  1:1978148/1‑150 (MQ=255)
                                                                 aCAGGGTCGAGATTAACATAAACAGAGGTTGTACCGGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGata                                                                               <  1:1150024/150‑1 (MQ=255)
                                                                   aGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGatata                                                                             <  2:592229/150‑1 (MQ=255)
                                                                     ggTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTta                                                                                                      <  2:305251/123‑1 (MQ=255)
                                                                     ggTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTta                                                                                                      >  1:305251/1‑123 (MQ=255)
                                                                     ggTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATc                                                                           >  1:1005008/1‑150 (MQ=255)
                                                                      gTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATcc                                                                          >  1:1009664/1‑150 (MQ=255)
                                                                             ttAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGtt                                                                   >  1:1175394/1‑150 (MQ=255)
                                                                               aaCATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTaa                                                                 >  2:1997843/1‑150 (MQ=255)
                                                                               aaCATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTaa                                                                 <  1:1068408/150‑1 (MQ=255)
                                                                                   aaaaCAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGACAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAAt                                                             <  2:1987420/149‑1 (MQ=255)
                                                                                     aaCAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAAttt                                                           >  2:580855/1‑149 (MQ=255)
                                                                                       cAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGAt                                                         >  2:1568149/1‑150 (MQ=255)
                                                                                          agGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTa                                                      <  2:1300119/150‑1 (MQ=255)
                                                                                              tGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTaaaaa                                                  <  2:939183/150‑1 (MQ=255)
                                                                                                tACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAaaga                                                >  1:678891/1‑150 (MQ=255)
                                                                                                 aCCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAaagaa                                               >  1:1276234/1‑150 (MQ=255)
                                                                                                   cTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTTCCATTAATAGAATATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAgg                                             <  1:454325/150‑1 (MQ=255)
                                                                                                    tGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGc                                            >  2:1633578/1‑150 (MQ=255)
                                                                                                      cTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAaaga                                                >  1:937479/1‑144 (MQ=255)
                                                                                                      cTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAaaga                                                <  2:937479/144‑1 (MQ=255)
                                                                                                        ccAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGAATCTTTTACACATGATATATAATAAGGTTAAAAATTGATTTAAAAAGAAGGATaat                                        >  2:1931276/1‑149 (MQ=255)
                                                                                                         cAATACTTGACCCGCCATCGACATTACAACTGAATGCAAAAACTTTGCCTGCCATTAATAGAATATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGa                                       <  1:1991164/150‑1 (MQ=255)
                                                                                                         cAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGa                                       >  1:1013939/1‑150 (MQ=255)
                                                                                                           aTACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGaaa                                     <  1:435441/150‑1 (MQ=255)
                                                                                                            tACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGaaaa                                    >  1:716302/1‑150 (MQ=255)
                                                                                                             aCTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAAAAaat                                   >  2:1615128/1‑150 (MQ=255)
                                                                                                             aCTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCAATAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAaat                                   >  2:447644/1‑150 (MQ=255)
                                                                                                               ttGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAaatcg                                 <  2:403170/150‑1 (MQ=255)
                                                                                                                 gACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAaatcgaa                               <  2:1439625/150‑1 (MQ=255)
                                                                                                                    ccGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAaatcgaatcg                            >  1:208298/1‑150 (MQ=255)
                                                                                                                      gCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGAc                          >  2:579526/1‑150 (MQ=255)
                                                                                                                       ccATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACg                         <  1:546071/150‑1 (MQ=255)
                                                                                                                        cATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGt                        >  2:1055832/1‑150 (MQ=255)
                                                                                                                          tCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTa                      <  1:1568149/150‑1 (MQ=255)
                                                                                                                           cGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTAt                     <  2:1009664/150‑1 (MQ=255)
                                                                                                                            taaaTTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTAtt                    <  1:101408/147‑1 (MQ=255)
                                                                                                                                  aCAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAgg              <  2:678891/150‑1 (MQ=255)
                                                                                                                                  aCAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAgg              <  1:1997843/150‑1 (MQ=255)
                                                                                                                                   cAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGgt             >  1:1793085/1‑150 (MQ=255)
                                                                                                                                     aCTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGgtgt           <  2:345476/150‑1 (MQ=255)
                                                                                                                                      tggcaGGCAAAGTTTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAatcgaatcga                           <  1:1352077/129‑1 (MQ=255)
                                                                                                                                      tggcaGGCAAAGTTTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAatcgaatcga                           >  2:1352077/5‑133 (MQ=255)
                                                                                                                                      cTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTa          <  2:1013939/150‑1 (MQ=255)
                                                                                                                                        gcaGGCAAAGTTTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGa        >  1:383607/3‑150 (MQ=255)
                                                                                                                                        gcaGGCAAAGTTTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGa        >  2:1762162/3‑150 (MQ=255)
                                                                                                                                        gcaGGCAAAGTTTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGa        >  1:795770/3‑150 (MQ=255)
                                                                                                                                         aaGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGAt       <  1:1097758/150‑1 (MQ=255)
                                                                                                                                             cAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGgtgt           >  2:84247/1‑142 (MQ=255)
                                                                                                                                             cAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGgtgt           <  1:84247/142‑1 (MQ=255)
                                                                                                                                             cAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATcaca   >  1:1402447/1‑150 (MQ=255)
                                                                                                                                              aaaaaCTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCCAATCTACGTTATTGCCAGGTGTATATCACAt  >  2:1781270/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CGTAATCATTCCAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACAT  >  NC_000913/1587524‑1587815

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