Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,639,617 T→C G147G (GGA→GGG rrrQ ← Qin prophage; putative lysozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,639,6170TC99.2% 434.5 / ‑4.5 119G147G (GGA→GGGrrrQQin prophage; putative lysozyme
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (67/51);  total (67/52)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 121 total aligned reads displayed.

ACATTTCCGACGATAAATCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGTCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTTGTGGTTCG  >  NC_000913/1639469‑1639764
                                                                                                                                                    |                                                                                                                                                   
aCATTTCCGACGATAAATCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGccc                                                                                                                                                   <  2:1178188/151‑1 (MQ=255)
 cATTTCCGACGATAAATCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCg                                                                                                                                                  >  1:199561/1‑151 (MQ=255)
     tCCGACGATAAATCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCa                                                                                                                                               <  2:1040860/150‑1 (MQ=255)
      ccGACGATAAATCCTGACAACACGAGAAAAATCACCATGAATATCTGATACACTGTTCTATCCCCCAGAAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCaac                                                                                                                                             >  3:195523/1‑149 (MQ=255)
       cgacgaTAAATCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATcc                                                                                                                                            >  3:1110002/1‑151 (MQ=255)
          cgaTAAATCCTGACAACACGAGAAAAATCGCCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCtt                                                                                                                                          <  3:1070209/150‑1 (MQ=255)
           gaTAAATCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTa                                                                                                                                         >  3:759256/1‑150 (MQ=255)
           gaTAAATCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTa                                                                                                                                         >  1:346757/1‑150 (MQ=255)
              aaaTCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATcc                                                                                                                                     <  2:768501/151‑1 (MQ=255)
                  ccTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATccacc                                                                                                                                  >  2:589766/1‑150 (MQ=255)
                    tGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGc                                                                                                                               >  3:626807/1‑151 (MQ=255)
                        aacaCGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAAt                                                                                                                           <  1:189693/151‑1 (MQ=255)
                        aacaCGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAAt                                                                                                                           <  3:314903/151‑1 (MQ=255)
                          cacGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCg                                                                                                                         <  2:438845/151‑1 (MQ=255)
                               gaAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCg                                                                                                                    >  2:1182327/1‑151 (MQ=255)
                                 aaaaTCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGc                                                                                                                   >  3:1368363/1‑150 (MQ=255)
                                 aaaaTCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGc                                                                                                                   >  3:807288/1‑150 (MQ=255)
                                 aaaaTCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCa                                                                                                                  <  2:626781/151‑1 (MQ=255)
                                    aTCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGc                                                                                                               >  3:1120816/1‑151 (MQ=255)
                                     tCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGc                                                                                                               >  3:1397210/1‑150 (MQ=255)
                                      cACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCAc                                                                                                             <  3:1246723/151‑1 (MQ=255)
                                        ccATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCt                                                                                                           >  4:16277/1‑151 (MQ=255)
                                         cATGAATATCTTATGCACTGGTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCtt                                                                                                          <  3:1060410/151‑1 (MQ=255)
                                         cATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCtt                                                                                                          <  3:861192/151‑1 (MQ=255)
                                            gAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTa                                                                                                        <  3:1210126/150‑1 (MQ=255)
                                            gAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTAc                                                                                                       >  2:1088403/1‑151 (MQ=255)
                                            gAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTAc                                                                                                       >  2:652783/1‑151 (MQ=255)
                                             aaTATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACg                                                                                                      >  3:688690/1‑151 (MQ=255)
                                             aaTATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACg                                                                                                      >  2:962307/1‑151 (MQ=255)
                                               tatCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGAt                                                                                                    >  3:822331/1‑151 (MQ=255)
                                                atCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTTACATTTATTTGAACGAACGCGGCAACCGCGTCCGCCATCCTTAATCCACCACCTAACCCCTTCGCATGCCCCTTCACGac                                                                                                    <  3:659033/150‑2 (MQ=255)
                                                atCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATc                                                                                                   >  4:468845/1‑151 (MQ=255)
                                                  cTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCAc                                                                                                 >  2:565440/1‑151 (MQ=255)
                                                       tCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGCCGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGc                                                                                             <  3:775348/150‑1 (MQ=255)
                                                       tCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGc                                                                                             >  2:221837/1‑150 (MQ=255)
                                                        cACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCa                                                                                            <  4:427185/150‑1 (MQ=255)
                                                        cACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCAt                                                                                           <  2:157448/151‑1 (MQ=255)
                                                           tGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTc                                                                                         >  3:1014926/1‑150 (MQ=255)
                                                              tCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGc                                                                                      <  3:957996/150‑1 (MQ=255)
                                                                   cccccAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTta                                                                                <  3:1208980/151‑1 (MQ=255)
                                                                        agcagGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATaaa                                                                            >  1:363227/1‑150 (MQ=255)
                                                                         gcagGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAAc                                                                          >  1:266231/1‑151 (MQ=255)
                                                                          cagGTTCGTGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACg                                                                         <  2:266444/151‑1 (MQ=255)
                                                                           agGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGt                                                                        <  2:163853/151‑1 (MQ=255)
                                                                              ttAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCg                                                                      >  3:386501/1‑150 (MQ=255)
                                                                                   gcgcTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACggg                                                                 >  2:597500/1‑150 (MQ=255)
                                                                                    cgcTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGa                                                                >  1:344660/1‑150 (MQ=255)
                                                                                      ctctCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCa                                                                                        >  1:236095/1‑124 (MQ=255)
                                                                                      ctctCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCa                                                                                        <  4:236092/124‑1 (MQ=255)
                                                                                      ctctCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGaaa                                                              >  4:435247/1‑150 (MQ=255)
                                                                                       tctcCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAAc                                                                          >  2:1064696/1‑137 (MQ=255)
                                                                                       tctcCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAAc                                                                          <  3:1064742/137‑1 (MQ=255)
                                                                                           cTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACtt                                                         >  3:132782/1‑150 (MQ=255)
                                                                                            tGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTAc                                                       >  3:680710/1‑151 (MQ=255)
                                                                                             ggTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTAcc                                                      <  4:26039/151‑1 (MQ=255)
                                                                                                cACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCggg                                                   >  3:640212/1‑151 (MQ=255)
                                                                                                cACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCggg                                                   <  3:387437/151‑1 (MQ=255)
                                                                                                   gacgaATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGcc                                                <  2:807253/151‑1 (MQ=255)
                                                                                                    acgaATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCa                                               <  2:822296/151‑1 (MQ=255)
                                                                                                    acgaATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCa                                               <  1:350724/151‑1 (MQ=255)
                                                                                                     cgaATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCa                                               >  2:256296/1‑150 (MQ=255)
                                                                                                         tAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGtt                                          >  1:170333/1‑151 (MQ=255)
                                                                                                          aaCCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGtt                                          >  1:64443/1‑150 (MQ=255)
                                                                                                               gACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAggg                                    >  3:935206/1‑151 (MQ=255)
                                                                                                               gACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAggg                                    >  4:13925/1‑151 (MQ=255)
                                                                                                                 ccGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGa                                   <  4:160787/150‑1 (MQ=255)
                                                                                                                 ccGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGAc                                  >  2:349881/1‑151 (MQ=255)
                                                                                                                   gTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACaa                                <  3:245252/151‑1 (MQ=255)
                                                                                                                       cAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATg                            >  1:246421/1‑151 (MQ=255)
                                                                                                                       cAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATg                            >  4:287381/1‑151 (MQ=255)
                                                                                                                        aGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAAt                             <  2:367393/149‑1 (MQ=255)
                                                                                                                        aGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGa                           <  3:1315725/151‑1 (MQ=255)
                                                                                                                          ttATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGAcg                         >  3:171628/1‑151 (MQ=255)
                                                                                                                          ttATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGAcg                         >  2:534888/1‑151 (MQ=255)
                                                                                                                           tATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGAcg                         <  3:1026887/150‑1 (MQ=255)
                                                                                                                            aTTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGAcgcg                       >  3:23147/1‑151 (MQ=255)
                                                                                                                            aTTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGAcgcg                       <  2:1397139/151‑1 (MQ=255)
                                                                                                                             tttGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGAcgcg                       >  3:1398316/1‑150 (MQ=255)
                                                                                                                              ttGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGAt                     >  2:388704/1‑151 (MQ=255)
                                                                                                                               tGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGAt                     <  2:935170/150‑1 (MQ=255)
                                                                                                                                gAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATa                    <  2:225148/150‑1 (MQ=255)
                                                                                                                                 aaCGAATGCGGCAATCGCGCCCGCCAACCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATAcc                  <  2:195494/151‑1 (MQ=255)
                                                                                                                                  aCGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATAccc                 >  3:180764/1‑151 (MQ=255)
                                                                                                                                  aCGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATAccc                 >  2:118436/1‑151 (MQ=255)
                                                                                                                                   cGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCg                >  2:785510/1‑151 (MQ=255)
                                                                                                                                    gAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGc               >  2:632077/1‑151 (MQ=255)
                                                                                                                                     aaTGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCACATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCt              >  2:297217/1‑151 (MQ=255)
                                                                                                                                        gCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCtttt           <  1:82487/151‑1 (MQ=255)
                                                                                                                                         cGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCttttt          <  2:576798/151‑1 (MQ=255)
                                                                                                                                         cGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCttttt          >  4:35610/1‑151 (MQ=255)
                                                                                                                                           gCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCttttt          <  2:1368294/149‑1 (MQ=255)
                                                                                                                                            cAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTtgtg       <  2:816390/151‑1 (MQ=255)
                                                                                                                                            cAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTtgtg       <  3:340311/151‑1 (MQ=255)
                                                                                                                                             aaTCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTTGTgg      >  3:1265873/1‑151 (MQ=255)
                                                                                                                                             aaTCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTTGTgg      >  2:1182593/1‑151 (MQ=255)
                                                                                                                                              aTCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTTGTGGt     >  2:898846/1‑151 (MQ=255)
                                                                                                                                                cgcgCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTTGTGGtt    <  4:368084/150‑1 (MQ=255)
                                                                                                                                                cgcgCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTTGTGGTTc   <  2:688665/151‑1 (MQ=255)
                                                                                                                                                cgcgCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTTGTGGTTc   <  4:416557/151‑1 (MQ=255)
                                                                                                                                                 gcgcCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTTGTGGTTCg  >  2:1203920/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACATTTCCGACGATAAATCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGTCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTTGTGGTTCG  >  NC_000913/1639469‑1639764

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: