Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,641,252 T→C intergenic (‑667/+87) essQ ← / ← cspB Qin prophage; putative S lysis protein/Qin prophage; cold shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,641,2520TC98.5% 268.0 / ‑4.0 68intergenic (‑667/+87)essQ/cspBQin prophage; putative S lysis protein/Qin prophage; cold shock protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (34/33);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

TGCCATGCCAGCAGTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATTAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTTGCTGCTGCAGGACCTTTAGCACCACTCTCTATAGAGA  >  NC_000913/1641106‑1641397
                                                                                                                                                  |                                                                                                                                                 
tgccatgccaGCAGTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCaa                                                                                                                                                 >  2:1133982/1‑149 (MQ=255)
   catgccaGCAGTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTTCTTTTTATCAAAgaga                                                                                                                                            <  2:1201501/151‑1 (MQ=255)
   catgccaGCAGTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAag                                                                                                                                               >  2:228661/1‑148 (MQ=255)
      gccaGCAGTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATcc                                                                                                                                         <  2:1030887/151‑1 (MQ=255)
       ccaGCAGTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCa                                                                                                                                        <  2:332556/151‑1 (MQ=255)
         agcagTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAat                                                                                                                                      <  2:104507/151‑1 (MQ=255)
         agcagTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAat                                                                                                                                      <  2:615460/151‑1 (MQ=255)
          gcagTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAat                                                                                                                                      >  1:568282/1‑150 (MQ=255)
            agTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATAtt                                                                                                                                   >  1:84585/1‑151 (MQ=255)
              ttAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAat                                                                                                                                      >  2:730802/1‑146 (MQ=255)
              ttAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAat                                                                                                                                      <  1:730777/146‑1 (MQ=255)
                aGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCac                                                                                                                                <  1:339318/150‑1 (MQ=255)
                aGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCa                                                                                                                                 <  2:390654/149‑1 (MQ=255)
                   ttCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCactact                                                                                                                            >  2:603951/1‑151 (MQ=255)
                       gctgTAAAACTATTCATGCAGCAAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGGCCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTActct                                                                                                                          <  2:1090650/149‑1 (MQ=255)
                       gctgTAAAACTATTCATGCAGCAAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGt                                                                                                                        >  2:1213498/1‑151 (MQ=255)
                         tgTAAAACTATTCATGCAGCAAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTcc                                                                                                                      >  2:618273/1‑151 (MQ=255)
                                  aTTCATGCAGCAAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTAtctc                                                                                                               <  1:1146130/149‑1 (MQ=255)
                                  aTTCATGCAGCAAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATctct                                                                                                              >  2:473223/1‑150 (MQ=255)
                                      aTGCAGCAAACCTGCACTTCTCCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTCACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCa                                                                                                         <  2:894992/151‑1 (MQ=255)
                                        gcagcaAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAgg                                                                                                       >  2:1218329/1‑151 (MQ=255)
                                            caAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATc                                                                                                   >  2:1122669/1‑151 (MQ=255)
                                             aaaCCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCa                                                                                                  <  2:225429/151‑1 (MQ=255)
                                                  tGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCctt                                                                                             >  2:501888/1‑151 (MQ=255)
                                                  tGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCctt                                                                                             >  2:1192096/1‑151 (MQ=255)
                                                   gCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCct                                                                                              <  1:855754/149‑1 (MQ=255)
                                                      cTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCttcttc                                                                                         >  2:148625/1‑151 (MQ=255)
                                                         caccacAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGtt                                                                                      <  2:628646/151‑1 (MQ=255)
                                                          accacAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGtt                                                                                      >  2:260677/1‑150 (MQ=255)
                                                           ccacAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTAt                                                                                    >  2:446702/1‑151 (MQ=255)
                                                           ccacAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTAt                                                                                    <  1:832195/151‑1 (MQ=255)
                                                             acaATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCg                                                                                  >  2:158644/1‑151 (MQ=255)
                                                                aTGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGtat                                                                               >  2:300434/1‑151 (MQ=255)
                                                                   gTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATAc                                                                             >  2:424231/1‑150 (MQ=255)
                                                                      aGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGa                                                                          <  2:773704/150‑1 (MQ=255)
                                                                       gCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACg                                                                        >  2:379819/1‑151 (MQ=255)
                                                                         aTACTTTTCCTGATTAAGATTTTGCCAAATGTGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGCCATCAGCCTTCTACGTTATCGTATACAGACGAg                                                                      <  1:643357/151‑1 (MQ=255)
                                                                         aTACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAg                                                                      <  2:322598/151‑1 (MQ=255)
                                                                           aCTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCg                                                                    <  2:412134/151‑1 (MQ=255)
                                                                                tCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGa                                                                <  2:1035020/150‑1 (MQ=255)
                                                                                    gATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAAttt                                                            <  1:390227/150‑1 (MQ=255)
                                                                                    gATTAAGATTTTGCCAAATATGCTAGCCATGGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAAtttt                                                           <  2:682151/151‑1 (MQ=255)
                                                                                         aGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATc                                                       >  1:316542/1‑150 (MQ=255)
                                                                                          gATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATc                                                       >  1:694032/1‑149 (MQ=255)
                                                                                           aTTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGGATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGt                                                    <  1:459580/151‑1 (MQ=255)
                                                                                             tttGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTa                                                   <  2:1036583/150‑1 (MQ=255)
                                                                                              ttGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAat                                                 >  1:1083442/1‑151 (MQ=255)
                                                                                                gCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAatg                                                >  2:355377/1‑150 (MQ=255)
                                                                                                gCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAatg                                                <  1:117156/150‑1 (MQ=255)
                                                                                                gCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAatg                                                <  2:1009847/150‑1 (MQ=255)
                                                                                                  cAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAatgatg                                             <  1:276956/151‑1 (MQ=255)
                                                                                                   aaaTATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAatgatg                                             >  2:625479/1‑150 (MQ=255)
                                                                                                    aaTATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGAc                                           <  1:326365/151‑1 (MQ=255)
                                                                                                        tGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACAttt                                       >  1:320659/1‑151 (MQ=255)
                                                                                                                 ttGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTtgctgctgc                               >  1:216706/1‑150 (MQ=255)
                                                                                                                  tGTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTtgctgctgc                               >  1:66582/1‑149 (MQ=255)
                                                                                                                   gTTTCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTTGCTGCTGCAgg                            >  2:772779/1‑151 (MQ=255)
                                                                                                                    tttCATGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTTGCTGCTGCAgg                            >  2:2522/1‑150 (MQ=255)
                                                                                                                        aTGTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTTGCTGCTGCAGGACCtt                       <  2:572080/151‑1 (MQ=255)
                                                                                                                          gTATTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTTGCTGCTGCAGGACCttt                      <  2:96558/150‑1 (MQ=255)
                                                                                                                            aTTGGACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTTGCTGCTGCAGGACCTTTAGc                   >  2:795322/1‑151 (MQ=255)
                                                                                                                               ggACCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTTGCTGCTGCAGGACCTTTAGcacc                >  2:375371/1‑151 (MQ=255)
                                                                                                                                 aCCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTTGCTGCTGCAGGACCTTTAGcaccac              <  1:935160/151‑1 (MQ=255)
                                                                                                                                 aCCTCCTTACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTTGCTGCTGCAGGACCTTTAGcaccac              >  2:1205754/1‑151 (MQ=255)
                                                                                                                                    tcctTACTTTTTATTAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTTGCTGCTGCAGGACCTTTAGCACCActct           >  2:387266/1‑151 (MQ=255)
                                                                                                                                        tACTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTTGCTGCTGCAGGACCTTTAGCACCACTCTCta        <  1:1213431/150‑1 (MQ=255)
                                                                                                                                         aCTTTTTATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTTGCTGCTGCAGGACCTTTAGCACCACTCTCtata      <  1:505390/151‑1 (MQ=255)
                                                                                                                                           tttttATCAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTTGCTGCTGCAGGACCTTTAGCACCACTCTCtat       >  1:1016565/1‑148 (MQ=255)
                                                                                                                                              ttatcaAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTTGCTGCTGCAGGACCTTTAGCACCACTCTCTATAgaga  <  1:20089/150‑1 (MQ=255)
                                                                                                                                               tatcaAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTTGCTGCTGCAGGACCTTTAGCACCACTCTCTATAgaga  >  2:351868/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGCCATGCCAGCAGTTAGCTTCTGCTGTAAAACTATTCATGCAGCAAACCTGCACTTCACCACAATGGTTAGCATACTTTTCCTGATTAAGATTTTGCCAAATATGCTAGCCATTGTTTCATGTATTGGACCTCCTTACTTTTTATTAAAGAGATCCAATATTCACTACTCTGTCCGTATCTCTACTCAGGCATCAGCCTTCTTCGTTATCGTATACAGACGAGCGATGAATTTTAATCAGTAATGATGACATTTGCTGCTGCAGGACCTTTAGCACCACTCTCTATAGAGA  >  NC_000913/1641106‑1641397

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