Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 177,409 G→C 6.4% S72R (AGC→AGG) ‡ yadS ← conserved inner membrane protein YadS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913177,4090GC 6.4% 257.6 / 6.2 78S72R (AGC→AGG) ‡yadSconserved inner membrane protein YadS
Reads supporting (aligned to +/- strand):  ref base G (40/33);  new base C (2/3);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

TACAAACCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAGGGTCCATACGCAATT  >  NC_000913/177261‑177554
                                                                                                                                                    |                                                                                                                                                 
tACAAACCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGc                                                                                                                                                  >  2:1556355/1‑150 (MQ=255)
tACAAACCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGc                                                                                                                                                  <  1:404640/150‑1 (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                                                                                            <  2:1452308/148‑1 (MQ=255)
                                                        ccGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGc                                                                                            >  1:1452308/1‑148 (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                                        <  2:490150/147‑1 (MQ=255)
                                                                                                             tAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCAc                                        >  1:490150/1‑147 (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:1283345/1‑150 (MQ=255)
                                                                                                                      tggctggcGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAg                            <  1:5026/150‑1 (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           >  1:544091/1‑150 (MQ=255)
                                                                                                                                       cGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAGGGTCCa           <  2:739387/150‑1 (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)
                                                                                                                                                    |                                                                                                                                                 
TACAAACCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACACCAAAAGGGTCCATACGCAATT  >  NC_000913/177261‑177554

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: