Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,052,264 A→G 5.1% intergenic (+24/+25) 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,2640AG 5.1% 238.9 / 3.2 79intergenic (+24/+25)gnsA/yccMputative phosphatidylethanolamine synthesis regulator/putative 4Fe‑4S membrane protein
Reads supporting (aligned to +/- strand):  ref base A (30/44);  new base G (2/2);  total (32/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

TATCGCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGT  >  NC_000913/1052117‑1052412
                                                                                                                                                   |                                                                                                                                                    
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:703237/1‑150 (MQ=255)
                                                                                        cTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTg                                                            >  1:1269885/1‑150 (MQ=255)
                                                                                         ttGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGCCCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGc                                                           <  1:1236254/150‑1 (MQ=255)
                                                                                          tGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAAt                                                                                                                  <  1:143888/94‑1 (MQ=255)
                                                                                          tGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAAt                                                                                                                  >  2:143888/1‑94 (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                                        <  1:647988/150‑1 (MQ=255)
                                                                                                            aatcaaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTg                                        <  2:575651/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               >  2:1214534/1‑150 (MQ=255)
                                                                                                                                     cATCACCTCGCACGGGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAAt               >  1:403872/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     <  2:364796/150‑1 (MQ=255)
                                                                                                                                               tGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAAcc     <  1:1027876/150‑1 (MQ=255)
                                                                                                                                               tGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACAAATCACCGCAAACATAATCCCCAGAAcc     <  2:145864/150‑1 (MQ=255)
                                                                                                                                               ggCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGACACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAAcc     <  2:1039508/149‑1 (MQ=255)
                                                                                                                                                 cGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAg   >  1:791136/1‑150 (MQ=255)
                                                                                                                                                 cGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAg   <  1:750698/150‑1 (MQ=255)
                                                                                                                                                  gAGGTGATGTTCCGGTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGt  <  1:383241/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TATCGCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGT  >  NC_000913/1052117‑1052412

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: