Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,587,669 A→G 6.4% V22V (GTT→GTC) ‡ 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,6690AG 6.4% 284.1 / 7.3 94V22V (GTT→GTC) ‡ydeQputative fimbrial adhesin protein YdeQ
Reads supporting (aligned to +/- strand):  ref base A (41/47);  new base G (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 = 8.65e-01

CCGCCGTAATCATTCCAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACAT  >  NC_000913/1587520‑1587815
                                                                                                                                                     |                                                                                                                                                  
ccgccgTAATCATTCCAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCaaaa                                                                                                                                                    <  2:1472795/150‑1 (MQ=255)
    cgTAATCATTCCAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACtt                                                                                                                                                >  1:1001442/1‑150 (MQ=255)
     gTAATCATTCCAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACttt                                                                                                                                               <  2:474342/150‑1 (MQ=255)
          cATTCCAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTtgcctg                                                                                                                                          >  2:1224896/1‑150 (MQ=255)
          cATTCCAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTtgcctg                                                                                                                                          >  1:403170/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                                                                                                                              <  1:639480/150‑1 (MQ=255)
                      ttATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATagaa                                                                                                                              <  2:115964/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                                                                                                                 >  2:285310/1‑150 (MQ=255)
                                   cAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGaa                                                                                                                 >  1:1501145/1‑150 (MQ=255)
                                       tCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAgg                                                                                                             <  2:1346791/150‑1 (MQ=255)
                                       tCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAgg                                                                                                             <  2:1464804/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                                                                 <  1:1068408/150‑1 (MQ=255)
                                                                                   aaCATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTaa                                                                 >  2:1997843/1‑150 (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              <  1:1997843/150‑1 (MQ=255)
                                                                                                                                      aCAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAgg              <  2:678891/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:795770/3‑150 (MQ=255)
                                                                                                                                            gcaGGCAAAGTTTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGa        >  2:1762162/3‑150 (MQ=255)
                                                                                                                                            gcaGGCAAAGTTTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGa        >  1:383607/3‑150 (MQ=255)
                                                                                                                                             aaGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGAt       <  1:1097758/150‑1 (MQ=255)
                                                                                                                                                 cAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGgtgt           <  1:84247/142‑1 (MQ=255)
                                                                                                                                                 cAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGgtgt           >  2:84247/1‑142 (MQ=255)
                                                                                                                                                 cAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATcaca   >  1:1402447/1‑150 (MQ=255)
                                                                                                                                                  aaaaaCTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCCAATCTACGTTATTGCCAGGTGTATATCACAt  >  2:1781270/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCGCCGTAATCATTCCAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACAT  >  NC_000913/1587520‑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: