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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,587,6700AT 5.2% 349.7 / 7.2 115V22D (GTT→GAT) ‡ydeQputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base A (52/57);  new base T (3/3);  total (55/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

Only 100 of 123 total aligned reads displayed.

CCGTAATCATTCCAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATT  >  NC_000913/1587523‑1587816
                                                                                                                                                   |                                                                                                                                                  
  gTAATCATTCCAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACtt                                                                                                                                                 <  1:1464477/149‑1 (MQ=255)
  gTAATCATTCCAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACtt                                                                                                                                                 <  1:951062/149‑1 (MQ=255)
     aTCATTCCAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCaaaaa                                                                                                                                                    <  1:1131151/143‑1 (MQ=255)
       ctatCCATCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGTGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTtgcct                                                                                                                                            <  2:120807/146‑1 (MQ=255)
         ttCCAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTtgcctgc                                                                                                                                          <  2:572716/149‑1 (MQ=255)
         ttCCAGCAACTTATATGCTGAGACAAGTCTACAACCAGAGTCTGTCCCGGCTGTATTACAGGGTCGAGATTAACATAATCAGAGGTTGTACCTTCTCCAATAATTGAACCGCGATCGACATTTCAACTGAAGGCAAAAACTTagcctgc                                                                                                                                          >  1:248554/1‑149 (MQ=255)
           ccAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAGTTTTtgccttcag                                                                                                                                        <  1:346309/149‑5 (MQ=255)
           ccAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAGTTTTtgccttcag                                                                                                                                        <  2:1153148/149‑5 (MQ=255)
            cAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCAt                                                                                                                                       <  2:1590751/149‑1 (MQ=255)
            cAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCAt                                                                                                                                       >  2:1336224/1‑149 (MQ=255)
            cAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCAt                                                                                                                                       <  1:1339701/149‑1 (MQ=255)
             aGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCAtt                                                                                                                                      >  1:1002874/1‑149 (MQ=255)
               cAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTaa                                                                                                                                    >  2:969646/1‑149 (MQ=255)
                aaCTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAAt                                                                                                                                   <  1:897854/149‑1 (MQ=255)
                aaCTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAAAACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCAATTAAt                                                                                                                                   >  2:1316347/1‑149 (MQ=255)
                 aCTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATa                                                                                                                                  <  1:310614/149‑1 (MQ=255)
                       atGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGAt                                                                                                                            >  1:1462153/1‑149 (MQ=255)
                        tGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATa                                                                                                                           >  2:869877/1‑149 (MQ=255)
                          cTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATaaa                                                                                                                         <  2:918749/149‑1 (MQ=255)
                             agaCAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATAc                                                                                                                      >  1:382592/1‑149 (MQ=255)
                              gaCAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATAcc                                                                                                                     >  2:624542/1‑149 (MQ=255)
                               aCAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCg                                                                                                                    >  1:1278714/1‑149 (MQ=255)
                                 aaGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTtgcctg                                                                                                                                           >  1:1454982/1‑124 (MQ=255)
                                 aaGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTtgcctg                                                                                                                                           <  2:1454982/124‑1 (MQ=255)
                                 aaGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGaa                                                                                                                  >  1:279249/1‑149 (MQ=255)
                                 aaGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGaa                                                                                                                  >  2:230858/1‑149 (MQ=255)
                                    tCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAg                                                                                                               >  1:1497684/1‑149 (MQ=255)
                                     gatctaCCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAgg                                                                                                              <  1:393747/144‑1 (MQ=255)
                                     cTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATATCGAACAgg                                                                                                              <  1:1336224/149‑1 (MQ=255)
                                      tACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGa                                                                                                             >  1:1168855/1‑149 (MQ=255)
                                       aCAAGCAGATTCTGGCCCGGTTGTATAACAGTGTCGAGATTAAAATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCAGCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGAc                                                                                                            <  1:302589/149‑1 (MQ=255)
                                         aaCCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCt                                                                                                          <  2:382592/149‑1 (MQ=255)
                                         aaCCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACAGACCAGGACCt                                                                                                          >  1:1209704/1‑149 (MQ=255)
                                          aCCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCtt                                                                                                         <  1:941508/149‑1 (MQ=255)
                                           ccAGATTATGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCAGTAATAGAAGATAAATACCGAACAGGACCttt                                                                                                        >  2:415788/1‑149 (MQ=255)
                                            cAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTta                                                                                                       >  1:438430/1‑149 (MQ=255)
                                            cAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTta                                                                                                       <  2:1416506/149‑1 (MQ=255)
                                              gATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGTCATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTtata                                                                                                     <  2:932726/149‑1 (MQ=255)
                                                ttCTGGACCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGAAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACt                                                                                                   <  1:192891/149‑1 (MQ=255)
                                                  cTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGa                                                                                                 <  2:1497684/149‑1 (MQ=255)
                                                       ccGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGtt                                                                                            >  2:1157388/1‑149 (MQ=255)
                                                       ccGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGtt                                                                                            <  2:1002874/149‑1 (MQ=255)
                                                         ggCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGtttt                                                                                          >  2:1036708/1‑149 (MQ=255)
                                                            tGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTATAcc                                                                                       >  2:209239/1‑149 (MQ=255)
                                                            tGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGAACTTTATACTGATTGTTTTAcc                                                                                       >  2:26541/1‑149 (MQ=255)
                                                              tATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGCCCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCa                                                                                     >  2:590800/1‑149 (MQ=255)
                                                                ttACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATg                                                                                   <  2:1543518/149‑1 (MQ=255)
                                                                ttACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTAACCAtt                                                                                   >  2:105484/1‑148 (MQ=255)
                                                                    aGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGatat                                                                               <  2:1278714/149‑1 (MQ=255)
                                                                      ggTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGatatat                                                                             <  1:1330287/149‑1 (MQ=255)
                                                                               tAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGtt                                                                    <  1:1551469/149‑1 (MQ=255)
                                                                                   aTAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTaaaa                                                                >  1:760853/1‑149 (MQ=255)
                                                                                   aTAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTaaaa                                                                >  1:717747/1‑149 (MQ=255)
                                                                                   aTAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTaaaa                                                                >  1:529556/1‑149 (MQ=255)
                                                                                   aTAAACAGAGGTTGTACATGCTCCAATAATTGACCCTCCATCGACATTAAAACTTAAGGCAAAAACTTTGACTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGATTTACCCATGCTATATCATAAGGTTaaaa                                                                >  1:1217961/1‑149 (MQ=255)
                                                                                    tAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTaaaaa                                                               >  1:1272345/1‑149 (MQ=255)
                                                                                     aaaCAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAAt                                                              <  2:1306040/149‑1 (MQ=255)
                                                                                     aaaCAGAGGTTGTACCTGCTCCAATACTTGACCAGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTAACCCATAATATATCATAAGGTTAAAAAt                                                              >  1:1071500/1‑149 (MQ=255)
                                                                                       aCAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATCATATATCCTAAGGTTAAAAAttt                                                            >  2:875269/1‑148 (MQ=255)
                                                                                        cagattttGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGa                                                           <  2:1228237/143‑1 (MQ=255)
                                                                                        cAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGa                                                           <  1:570707/149‑1 (MQ=255)
                                                                                          cagGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGTACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGAtt                                                         <  2:1168855/148‑1 (MQ=255)
                                                                                           agGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGAttt                                                        >  1:557967/1‑149 (MQ=255)
                                                                                              ttGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTaaa                                                     <  2:1410569/149‑1 (MQ=255)
                                                                                                  aCCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAaaga                                                 >  1:1048268/1‑149 (MQ=255)
                                                                                                   ccTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAaagaa                                                <  1:969646/149‑1 (MQ=255)
                                                                                                   ccTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTAATTTAAAaagaa                                                >  2:1508927/1‑149 (MQ=255)
                                                                                                             tACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGaaa                                      <  2:1355221/149‑1 (MQ=255)
                                                                                                                   aCCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAaatcgaa                                <  2:378776/149‑1 (MQ=255)
                                                                                                                   aCCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAaatcgaa                                <  2:57242/149‑1 (MQ=255)
                                                                                                                   aCCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAaatcgaa                                <  2:198677/149‑1 (MQ=255)
                                                                                                                    cccGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATAAATAGAAGATAAAAACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAaatcgaat                               <  2:623471/149‑1 (MQ=255)
                                                                                                                     ccGCCATCGACATTACAACTGAAGTCAAAAAATTTGCATGAAATTAATAGAATATAAATACCGAACAGTACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAaatcgaatc                              <  1:312028/149‑1 (MQ=255)
                                                                                                                     ccGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAaatcgaata                              >  2:1033155/1‑148 (MQ=255)
                                                                                                                        ccATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGAc                           >  1:1287384/1‑149 (MQ=255)
                                                                                                                             gACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTTATTTAAAAAGAAGGCTAAGAAAATCGAATCTACGTTAt                      >  2:956639/1‑149 (MQ=255)
                                                                                                                             gACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCTAATCGACGTTAt                      >  2:1004895/1‑149 (MQ=255)
                                                                                                                              aCATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTAtt                     <  1:1036708/149‑1 (MQ=255)
                                                                                                                              aCATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTCATTTAAAAAGAAGGCTAAGAAAATCGAATCGAAGTTAtt                     >  2:420423/1‑149 (MQ=255)
                                                                                                                                 ttACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGttt                                                                                           >  1:956081/1‑76 (MQ=255)
                                                                                                                                 ttACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGttt                                                                                           <  2:956081/76‑1 (MQ=255)
                                                                                                                                  tACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCa                 >  1:646683/1‑149 (MQ=255)
                                                                                                                                   aCAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAg                >  2:313147/1‑149 (MQ=255)
                                                                                                                                   aCAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAg                >  1:1384580/1‑149 (MQ=255)
                                                                                                                                   aCAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAg                <  1:590800/149‑1 (MQ=255)
                                                                                                                                     aaCTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGgt              >  1:1238048/1‑149 (MQ=255)
                                                                                                                                      atggcaGGCAAAGTTTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGgtg             >  1:1257570/6‑149 (MQ=255)
                                                                                                                                        tGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTTTTTTACCCATGATATATCCTAAGGTTAAAATTTGATTTAAAAAGAAGGCTAAGAAAATCTAATCGACGTTATTGCCAGGTGTa           <  2:529556/149‑1 (MQ=255)
                                                                                                                                          caGGCAAAGTTTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGa         >  2:787147/2‑149 (MQ=255)
                                                                                                                                          aaGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGa         <  1:122978/149‑1 (MQ=255)
                                                                                                                                          aaGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGa         <  2:646683/149‑1 (MQ=255)
                                                                                                                                          aaGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGa         <  2:450947/149‑1 (MQ=255)
                                                                                                                                           agtcAATAACTTTGCATGCAATTAATAGAAGATAAATACCGATCAGTACCTTTATACTGATTTTTTTACCCTTGAAATATCCTAAGGTTAAAAATTGATTTAAAAAGATGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGAt        <  1:105484/146‑1 (MQ=255)
                                                                                                                                           aGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGAt        <  2:279249/149‑1 (MQ=255)
                                                                                                                                           aGGCAAAAACTTTGCATGCAATTAATAGAATATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGAt        <  2:716714/149‑1 (MQ=255)
                                                                                                                                            ggCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATc       <  2:557967/149‑1 (MQ=255)
                                                                                                                                                aaaaCTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAATAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACAt   >  2:1509880/1‑149 (MQ=255)
                                                                                                                                                aaaaCTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTAAACTGATTGTTTTACACATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTATATCACAt   >  1:844714/1‑149 (MQ=255)
                                                                                                                                                 aaaCTTTGCCTGCCATTAATAGAAGATAAATACCGAACATGACCTTTATACTTATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTTATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACAtt  <  1:415788/149‑1 (MQ=255)
                                                                                                                                                 aaaCTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACAtt  <  2:438430/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCGTAATCATTCCAGCAACTTATATGCTGAGACAAGTCTACAACCAGATTCTGGCCCGGCTGTATTACAGGGTCGAGATTAACATAAACAGAGGTTGTACCTGCTCCAATACTTGACCCGCCATCGACATTACAACTGAAGGCAAAAACTTTGCCTGCCATTAATAGAAGATAAATACCGAACAGGACCTTTATACTGATTGTTTTACCCATGATATATCCTAAGGTTAAAAATTGATTTAAAAAGAAGGCTAAGAAAATCGAATCGACGTTATTGCCAGGTGTAGATCACATT  >  NC_000913/1587523‑1587816

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