Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 177,410 C→A 6.6% S72I (AGC→ATC) ‡ yadS ← conserved inner membrane protein YadS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913177,4100CA 6.6% 254.3 / 5.6 76S72I (AGC→ATC) ‡yadSconserved inner membrane protein YadS
Reads supporting (aligned to +/- strand):  ref base C (39/32);  new base A (2/3);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

TACAAACCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAGGGTCCATACGCAATTTTCCG  >  NC_000913/177261‑177559
                                                                                                                                                     |                                                                                                                                                     
tACAAACCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGc                                                                                                                                                       <  1:404640/150‑1 (MQ=255)
tACAAACCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGc                                                                                                                                                       >  2:1556355/1‑150 (MQ=255)
   aaaCCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCATTTTTGATAAAGATAGCCAGACCAACGTCTTCCAGAACATTCTACATGCATTTTGGTAAGGTTGTTGTTTGGATCACCAGCACGATGGTCAGCATGCTgg                                                                                                                                                    <  2:652267/150‑1 (MQ=255)
     aCCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTg                                                                                                                                                  >  1:274152/1‑150 (MQ=255)
         cgATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCa                                                                                                                                              <  1:755374/150‑1 (MQ=255)
         cgATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGACAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCa                                                                                                                                              <  2:1151639/150‑1 (MQ=255)
            ttAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCACCAGCATGCTGACCAtcg                                                                                                                                           <  1:773950/150‑3 (MQ=255)
              aaCGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCGAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCACCAGCATGCTGACCAtcgtg                                                                                                                                         <  1:1583263/150‑5 (MQ=255)
                 ggACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACg                                                                                                                                      <  1:1117513/150‑1 (MQ=255)
                 ggACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCACCAGCATGCTGACCAtcgtgctg                                                                                                                                      <  1:1594462/150‑8 (MQ=255)
                    ccGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGAcc                                                                                                                                   <  1:1544041/150‑1 (MQ=255)
                    ccGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGACTTGGCTGGCGCACCAGCACGATCGTCAGCATGCTGGTGACCATTGCAACGAcc                                                                                                                                   >  1:1142404/1‑150 (MQ=255)
                            cGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGt                                                                                                                           >  1:1196960/1‑150 (MQ=255)
                                ttAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGa                                                                                                                       <  2:1598775/150‑1 (MQ=255)
                                 tAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGAt                                                                                                                      >  2:1073453/1‑150 (MQ=255)
                                  aaaGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATc                                                                                                                     <  2:1214081/150‑1 (MQ=255)
                                  aaaGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATc                                                                                                                     >  2:570892/1‑150 (MQ=255)
                                     ggCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCttt                                                                                                                  >  1:963187/1‑150 (MQ=255)
                                      gCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTc                                                                                                                 >  1:904552/1‑150 (MQ=255)
                                          tATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCAccc                                                                                                             >  2:801409/1‑150 (MQ=255)
                                           aTTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCa                                                                                                            >  1:887247/1‑150 (MQ=255)
                                           aTTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGAACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGCCCAGATCGGTGGGATCTTTCACCCa                                                                                                            >  2:710159/1‑150 (MQ=255)
                                            ttCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCaa                                                                                                           >  2:1617941/1‑150 (MQ=255)
                                                  ccAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGGCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATAcc                                                                                                     <  2:1068208/150‑1 (MQ=255)
                                                    aaTGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCgg                                                                                                   <  1:1277928/150‑1 (MQ=255)
                                                       gCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGcc                                                                                                <  2:887247/150‑1 (MQ=255)
                                                       gCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCg                                                                                               >  1:36609/1‑150 (MQ=255)
                                                        ccGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGc                                                                                                 >  1:1452308/1‑148 (MQ=255)
                                                        ccGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGc                                                                                                 <  2:1452308/148‑1 (MQ=255)
                                                        ccGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCg                                                                                               >  2:571561/1‑150 (MQ=255)
                                                            cAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGAt                                                                                           >  2:1233831/1‑150 (MQ=255)
                                                              aaCACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATcc                                                                                         >  2:55543/1‑150 (MQ=255)
                                                              aaCACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTCGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTAACCCAAAATACCGGGCCGTGATcc                                                                                         >  1:209508/1‑150 (MQ=255)
                                                              aaCACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCAATTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATcc                                                                                         >  1:942451/1‑150 (MQ=255)
                                                                   cGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGcc                                                                                   <  2:36609/150‑1 (MQ=255)
                                                                           cAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTcgc                                                                            >  1:918687/1‑150 (MQ=255)
                                                                                   tCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTc                                                                    <  1:55543/150‑1 (MQ=255)
                                                                                     cAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTccc                                                                  >  1:877121/1‑150 (MQ=255)
                                                                                      aGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCCCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCccg                                                                 >  2:174964/1‑150 (MQ=255)
                                                                                       gCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCt                                                                                                                    <  2:338358/98‑1 (MQ=255)
                                                                                       gCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCt                                                                                                                    >  1:338358/1‑98 (MQ=255)
                                                                                           cGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCgccgcc                                                            <  1:610020/150‑1 (MQ=255)
                                                                                             gCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTa                                                          >  1:739387/1‑150 (MQ=255)
                                                                                                 cATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGc                                                      <  1:1590268/150‑1 (MQ=255)
                                                                                                  aTCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCg                                                     >  1:1315162/1‑150 (MQ=255)
                                                                                                   tCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCgg                                                    <  1:570892/150‑1 (MQ=255)
                                                                                                      aTTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTa                                                 <  1:860651/150‑1 (MQ=255)
                                                                                                      aTTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTa                                                 <  1:801409/150‑1 (MQ=255)
                                                                                                       ttttGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTAc                                                <  1:1014247/150‑1 (MQ=255)
                                                                                                            gTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGc                                           <  1:830753/150‑1 (MQ=255)
                                                                                                             tAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCAc                                             >  1:490150/1‑147 (MQ=255)
                                                                                                             tAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCAc                                             <  2:490150/147‑1 (MQ=255)
                                                                                                              aaGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGccc                                         <  1:674022/150‑1 (MQ=255)
                                                                                                                 ccatcgtgctggtGCGCCAGCACGATGGTCAGCATGCTGGTGCCCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCCGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGt                                      >  1:1476190/8‑150 (MQ=255)
                                                                                                                  catcgtgctggtGCGCCAGCACGATGGTCAGCATGCTGGTGACCATTGCACAGACCAGATCGGAGGGATCTTTCACCCAAAATCCCGGGCCGTGATACAGCGCCATGTCGCGAATTGTCGCGCCGCCTACTGCGGTTACCACGCCAAGTa                                     >  1:847489/7‑150 (MQ=255)
                                                                                                                   tCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTAc                                    <  1:710159/150‑1 (MQ=255)
                                                                                                                   atcgtgctggtGCGCCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTAc                                    >  1:1046836/6‑150 (MQ=255)
                                                                                                                      tggctggcGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAg                                 >  1:1274677/1‑150 (MQ=255)
                                                                                                                      tggctggcGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAg                                 <  1:5026/150‑1 (MQ=255)
                                                                                                                      tggctggcGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAg                                 >  1:1283345/1‑150 (MQ=255)
                                                                                                                          tggtGCGCCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAaca                             <  2:1046836/150‑1 (MQ=255)
                                                                                                                           ggcGCACCAGCACGATGGGCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAacac                            <  1:1077514/150‑1 (MQ=255)
                                                                                                                            gcgcACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCCGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACAcc                           <  2:1196960/150‑1 (MQ=255)
                                                                                                                             cgcACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGGCTACTGCGGTTACCACGCCCAGTACCAGAACACCa                          <  2:1142404/150‑1 (MQ=255)
                                                                                                                              tcACCAGCATGCTGACCATCGTGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCaa                         >  1:1519605/2‑150 (MQ=255)
                                                                                                                               cACCAGCATGCTGACCATCGTGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTAACCAAAAATACCGGGCCGTGATCCAGCGCCATGTAGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCaaa                        >  2:1104226/1‑150 (MQ=255)
                                                                                                                                 ccAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAg                      <  1:866906/150‑1 (MQ=255)
                                                                                                                                   aGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAggg                    >  2:916037/1‑150 (MQ=255)
                                                                                                                                      gcaaTGGTCACCAGCATGCTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGCTCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAGAGGGTcc                 >  2:589747/4‑150 (MQ=255)
                                                                                                                                      gcaaTGGTCACCAGCATGCTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAGGGTcc                 >  1:438963/4‑150 (MQ=255)
                                                                                                                                       cGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCTCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAGGGTCCa                >  2:1321244/1‑150 (MQ=255)
                                                                                                                                       cGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAGGGTCCa                <  2:739387/150‑1 (MQ=255)
                                                                                                                                       cGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAGGGTCCa                >  1:544091/1‑150 (MQ=255)
                                                                                                                                         aTGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAGGGTCCATa              <  1:1617941/150‑1 (MQ=255)
                                                                                                                                          tGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACCCCAAAAGGGTCCATAc             >  2:1601193/1‑150 (MQ=255)
                                                                                                                                            gTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAGGGTCCATACGc           >  1:1286597/1‑150 (MQ=255)
                                                                                                                                              cAGCATGCTGGTGACCATTTCAACGACCAGATCGGTGGGATCTTTAACACAAAATACCGGGCCGTTATACAGCGCCATGTATCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACAACAACACCAAAAGGGTCAATACGCaa         >  2:83159/1‑150 (MQ=255)
                                                                                                                                                gCATGCTGGTGACCATTGCAACGACCAGATCGCTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAGGGTCCATACGCAAtt       >  1:887971/1‑150 (MQ=255)
                                                                                                                                                     cTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAGGGTCCATACGCAATTTTCCg  >  1:133512/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACAAACCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAGGGTCCATACGCAATTTTCCG  >  NC_000913/177261‑177559

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: