Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,052,260 T→C 5.3% intergenic (+20/+29) gnsA → / ← yccM putative phosphatidylethanolamine synthesis regulator/putative 4Fe‑4S membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,052,2600TC 5.3% 217.5 / 3.2 76intergenic (+20/+29)gnsA/yccMputative phosphatidylethanolamine synthesis regulator/putative 4Fe‑4S membrane protein
Reads supporting (aligned to +/- strand):  ref base T (33/39);  new base C (2/2);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01

GACTTTATCGCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACC  >  NC_000913/1052112‑1052409
                                                                                                                                                    |                                                                                                                                                     
gACTTTATCGCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTg                                                                                                                                                      >  1:875875/1‑150 (MQ=255)
 aCTTTATCGCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGAGTTTGTGAACATCCCCCCGTGc                                                                                                                                                     >  1:110951/1‑150 (MQ=255)
  cTTTATCGCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAATTTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACACCGTGCg                                                                                                                                                    >  2:1131480/1‑150 (MQ=255)
  cTTTATCGCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCg                                                                                                                                                    >  2:1337878/1‑150 (MQ=255)
  cTTTATCGCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCg                                                                                                                                                    <  1:1299821/150‑1 (MQ=255)
     tATCGCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAgg                                                                                                                                                 <  1:325215/150‑1 (MQ=255)
        cgcgCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGa                                                                                                                                              <  2:272323/150‑1 (MQ=255)
          cgcAAAACATCGCCGAGCTGAACAAGAATACAGGGTAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATg                                                                                                                                            <  1:635147/150‑1 (MQ=255)
          cgcAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATg                                                                                                                                            <  1:1192389/150‑1 (MQ=255)
            cAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCGCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGtt                                                                                                                                          <  1:1287167/150‑1 (MQ=255)
              aaaaaTCGCCGAGCTGCACAAGAATCCAGCGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAAAGACCGGGCGTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCTCGCACGGGGTGATGTTca                                                                                                                                        <  1:649946/150‑2 (MQ=255)
              aaaaaTCGCCGAGCTGAACAAGAATACAGGGAAATAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAAACAAAATAATGTGATTATGTGAACATCACCCCGTGCGAGGTGATGTTcc                                                                                                                                        <  2:1116903/150‑1 (MQ=255)
              aaaaaTCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTcc                                                                                                                                        >  1:1010716/1‑150 (MQ=255)
               aaaaTCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCTCGCACGGGGTGATGTTcac                                                                                                                                       <  1:633297/150‑3 (MQ=255)
                aaaTCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGTTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGc                                                                                                                                      <  2:712518/150‑1 (MQ=255)
                aaaTCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGc                                                                                                                                      <  2:1248396/150‑1 (MQ=255)
                  aTCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCtt                                                                                                                                    <  1:779855/150‑1 (MQ=255)
                   tCGACGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCttg                                                                                                                                   <  2:864004/150‑1 (MQ=255)
                       cgagaggAACAAGAATACAGGGAAAGAAGTCTATGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGc                                                                                                                               <  1:906984/144‑1 (MQ=255)
                       cGAGCGGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGc                                                                                                                               <  1:330497/150‑1 (MQ=255)
                         aGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTa                                                                                                                             <  2:1049357/150‑1 (MQ=255)
                            tGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAAtt                                                                                                                          >  2:832362/1‑150 (MQ=255)
                             gAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCAGCTTGTTGCTAAttt                                                                                                                         >  1:145864/1‑150 (MQ=255)
                               acacGAATACACGGAAAGAAGTCTCTGAAATTCTCTTCACCGCACGAGAAATAATGACCGGGCTTGAAAGTTATGAAGTCAAAATCAAAATAATGTGATTTTGTGAACCTCACCCCGTGCGAGGGGATGTTCCGCTTCTTGCTAATTTAg                                                                                                                       <  2:963123/150‑1 (MQ=255)
                               acacGAATAACGGGCACGACGTCCCTGAACTAAGCTGCACCGCACGAGAAAAACTGACCGGGCTTGCCCGTTATTAGGGCCACATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTCATTTAg                                                                                                                       <  2:958397/150‑1 (MQ=255)
                                 aaGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTg                                                                                                                     <  1:801489/150‑1 (MQ=255)
                                 aaGAATACAGCGAAAGAAGTCTCTGAAATTCGCATCACTGAAAGAGAAAACATGAGCGGGCTTGAACGTTATGATGTCAAAATCATATTAATTTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTg                                                                                                                     <  1:200105/150‑1 (MQ=255)
                                   gAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGAc                                                                                                                   <  2:340469/150‑1 (MQ=255)
                                         aGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCa                                                                                                             <  2:80121/150‑1 (MQ=255)
                                           ggAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACACCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCAtt                                                                                                           >  1:935876/1‑150 (MQ=255)
                                                 aagTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGTGGTGATGTTGCGCTTGTTGCTAATTTTGTGTCCAATCATTGGCGCt                                                                                                     <  2:879187/150‑1 (MQ=255)
                                                    tctcTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTtgt                                                                                                  >  2:1241249/1‑150 (MQ=255)
                                                      tctGAAATTCGCTTCACAGCACGAGAAAAACTGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAAGGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGCCAAATCATTGGAGCAGGTgg                                                                                                >  2:903186/1‑150 (MQ=255)
                                                         taaaTTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTTATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAAt                                                                                             <  2:983964/149‑1 (MQ=255)
                                                         gAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAAt                                                                                             <  2:725325/150‑1 (MQ=255)
                                                           aaTTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTa                                                                                           <  2:425093/150‑1 (MQ=255)
                                                                 cTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTc                                                                                     >  2:728828/1‑150 (MQ=255)
                                                                   tCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGCGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCgg                                                                                   >  1:286938/1‑150 (MQ=255)
                                                                    cACCGCACGAGAAAAAATGACCGGGCTTGAAAGCTATGATGTCAAAATCAAAATAATGCGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGt                                                                                  <  1:1349290/150‑1 (MQ=255)
                                                                        gCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACaa                                                                              >  1:852895/1‑150 (MQ=255)
                                                                         cACGAGAAAAAATGCCCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACACCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAAt                                                                             <  1:780800/150‑1 (MQ=255)
                                                                         cACGAGAAAAAATGACCGGGCTTGATAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATATCCACGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAAt                                                                             >  2:855186/1‑150 (MQ=255)
                                                                           cGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATAAAAATAATGTGATATTGTGAACATCACCACGTGCTAGGTGATATTCCGCTTGTTGATAATTTAGTGACAAAACATTGGCTCTATTGGAATTAAGCGTCGGTTCAATtc                                                                           >  1:964012/1‑150 (MQ=255)
                                                                            gagaAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTAGTTGCAAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATtcc                                                                          >  2:94706/1‑150 (MQ=255)
                                                                               aaaaaaTGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATtcctcc                                                                       <  1:1181053/150‑1 (MQ=255)
                                                                                aaaaaTGACCGGGCTTGATAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTTGCGCTTGTGGAATTAAGCGTCGGTACAATtcctcag                                                                      >  2:474114/1‑148 (MQ=255)
                                                                                 aaaaTGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCgg                                                                     <  1:1227408/150‑1 (MQ=255)
                                                                                        ccGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCggg                                                              >  1:668394/1‑150 (MQ=255)
                                                                                             cTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTg                                                         >  1:1269885/1‑150 (MQ=255)
                                                                                             cTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTg                                                         >  1:703237/1‑150 (MQ=255)
                                                                                              ttGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGCCCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGc                                                        <  1:1236254/150‑1 (MQ=255)
                                                                                               tGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAAt                                                                                                               >  2:143888/1‑94 (MQ=255)
                                                                                               tGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAAt                                                                                                               <  1:143888/94‑1 (MQ=255)
                                                                                               tGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGcc                                                       <  2:487407/150‑1 (MQ=255)
                                                                                                    gTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACt                                                  >  2:101240/1‑150 (MQ=255)
                                                                                                    gTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACt                                                  >  1:281747/1‑150 (MQ=255)
                                                                                                                 aatcaaaatAATGTTATTTTGTGAACATCACCCCTTGCGATGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTg                                     <  1:107691/150‑1 (MQ=255)
                                                                                                                 aatcaaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTg                                     <  2:575651/150‑1 (MQ=255)
                                                                                                                 aatcaaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTg                                     <  1:647988/150‑1 (MQ=255)
                                                                                                                   tcaatatAATGTGATTTTGTGAACATCAACCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCCTACTCCCGAATACATTGCg                                   >  2:1006718/1‑150 (MQ=255)
                                                                                                                     aaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTa                                 <  2:569435/150‑1 (MQ=255)
                                                                                                                     aaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGATAATTTAGTAACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTAATCCGGCACCGGGCATTGCCATACTACAGCATAAATTGCGTa                                 >  2:663594/1‑150 (MQ=255)
                                                                                                                     aaaatAATGTGATTTTGTGAACACCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTa                                 >  1:1296095/1‑150 (MQ=255)
                                                                                                                      aaatAATGTGATTTTGTGAACATCCCCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTACGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACAATGCGGaa                                >  1:831569/1‑150 (MQ=255)
                                                                                                                              tgATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCAc                        <  1:1337878/150‑1 (MQ=255)
                                                                                                                               gATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCAcc                       >  2:205485/1‑150 (MQ=255)
                                                                                                                               gATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTAACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCAcc                       >  1:163010/1‑150 (MQ=255)
                                                                                                                                  tttGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCa                    <  2:2447/150‑1 (MQ=255)
                                                                                                                                       gAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACAt               >  1:713031/1‑150 (MQ=255)
                                                                                                                                        aaCATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTACTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATa              >  2:975827/1‑150 (MQ=255)
                                                                                                                                         aCATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATaa             >  1:1122163/1‑150 (MQ=255)
                                                                                                                                          cATCACCTCGCACGGGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAAt            >  1:403872/1‑150 (MQ=255)
                                                                                                                                          cATCACCTCGCACGGGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAAt            >  2:1214534/1‑150 (MQ=255)
                                                                                                                                            tCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAATCATAATcc          >  2:861501/1‑150 (MQ=255)
                                                                                                                                             cACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATccc         <  1:710050/150‑1 (MQ=255)
                                                                                                                                              aCCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGTCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATcccc        <  1:832362/150‑1 (MQ=255)
                                                                                                                                                    tGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCATTCATTGGCGCTTTTGGAATTAAGCGTCGGTACAATTGCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAAcc  <  2:626105/150‑1 (MQ=255)
                                                                                                                                                    tGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAAcc  <  1:1027876/150‑1 (MQ=255)
                                                                                                                                                    tGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAAcc  <  2:364796/150‑1 (MQ=255)
                                                                                                                                                    tGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACAAATCACCGCAAACATAATCCCCAGAAcc  <  2:145864/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GACTTTATCGCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACC  >  NC_000913/1052112‑1052409

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: