Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,768,620 G→A 16.3% intergenic (+47/‑45) yfjP → / → yfjQ CP4‑57 prophage; 50S ribosome‑binding GTPase family protein/CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,768,6200GA16.3% 511.1 / 70.6 204intergenic (+47/‑45)yfjP/yfjQCP4‑57 prophage; 50S ribosome‑binding GTPase family protein/CP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (81/88);  new base A (18/15);  total (101/103)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01

Only 100 of 208 total aligned reads displayed.

TGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTC  >  NC_000913/2768521‑2768719
                                                                                                   |                                                                                                   
 gCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGc                                                                                                    <  2:734807/100‑1 (MQ=255)
   gTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAg                                                                                                  >  2:3249860/1‑100 (MQ=255)
    ttCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAgg                                                                                                 <  1:466606/100‑1 (MQ=25)
    ttCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAgg                                                                                                 <  2:2204878/100‑1 (MQ=25)
     tCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTTTTCTCCGAGCCCTGTCGCAAGTGCGGCAggg                                                                                                >  1:2372630/1‑100 (MQ=255)
     tCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAggg                                                                                                <  2:1014729/100‑1 (MQ=255)
     tCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAggg                                                                                                <  2:2719104/100‑1 (MQ=255)
     tCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAggg                                                                                                <  2:4214603/100‑1 (MQ=25)
     tCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAggg                                                                                                <  2:3518837/100‑1 (MQ=25)
     tCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCGTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAggg                                                                                                <  2:2123716/100‑1 (MQ=255)
      ccTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGc                                                                                               <  2:889474/100‑1 (MQ=255)
      ccTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCAAGCCCTGCCGCACTTGCGACAGGGc                                                                                               <  1:1645515/100‑1 (MQ=21)
       cTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAGGGct                                                                                              <  2:1846873/100‑1 (MQ=25)
       cTCGGTGGTGTCAGTGGCCCGCGCGGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGACCCCTGCCGCACTTGCGACAggg                                                                                                <  1:679061/98‑1 (MQ=21)
        tCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGTCCTGCCGCACTTGCGACAGGCctc                                                                                             <  1:1995574/100‑2 (MQ=21)
        tCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctt                                                                                             <  2:2459694/100‑1 (MQ=255)
          ggtggtGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttc                                                                                           >  2:2988795/1‑100 (MQ=255)
          ggtggtGGCAGTGGCCCGCGCTGGCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGCCGCACTTGCGACAggg                                                                                                <  1:1526051/95‑1 (MQ=25)
           gtggtgTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttct                                                                                          <  2:3749150/100‑1 (MQ=255)
           gtggtgTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGctttct                                                                                          <  1:2232701/100‑1 (MQ=255)
              gtgtCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTa                                                                                       <  2:1628867/100‑1 (MQ=255)
               tgtCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTAt                                                                                      >  2:4084902/1‑100 (MQ=255)
               tgtCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTAt                                                                                      >  1:2477719/1‑100 (MQ=255)
                gtCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTAtt                                                                                     <  1:2519917/100‑1 (MQ=255)
                  cAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTg                                                                                   <  2:3350677/100‑1 (MQ=255)
                  cAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTg                                                                                   <  1:3243314/100‑1 (MQ=255)
                   aGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGt                                                                                  >  2:2388173/1‑100 (MQ=255)
                      ggCCCGCGCTGTCTGGGATTTCTTCTTCTGCGTGTTTAATTTTGCTCTATTATTCTCCGAGCCCTGTCGCAAGTTCGTCTCAGCTTTCTTTATTTGtttt                                                                               >  1:1954034/1‑100 (MQ=37)
                      ggCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGtttt                                                                               <  2:992460/100‑1 (MQ=255)
                       gCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTa                                                                              >  1:1174917/1‑100 (MQ=255)
                         ccGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTAcc                                                                            <  1:2112416/100‑1 (MQ=255)
                         ccGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTAcc                                                                            <  1:4018721/100‑1 (MQ=255)
                         ccGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTAcc                                                                            <  1:1101294/100‑1 (MQ=255)
                          cgcgCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTAccc                                                                           <  1:3305079/100‑1 (MQ=255)
                          cgcgCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTAccc                                                                           <  2:328489/100‑1 (MQ=255)
                               tGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACt                                                                      >  1:801063/1‑100 (MQ=255)
                                gTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTg                                                                     >  1:1223760/1‑100 (MQ=255)
                                  cTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgag                                                                   <  2:2494669/100‑1 (MQ=255)
                                  cTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgag                                                                   >  2:2582924/1‑100 (MQ=255)
                                    gggATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgagga                                                                 >  2:147279/1‑100 (MQ=255)
                                    gggATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTgagg                                                                  <  1:428158/99‑1 (MQ=255)
                                      gATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGt                                                               >  1:1061934/1‑100 (MQ=255)
                                       aTTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTc                                                              <  2:1911815/100‑1 (MQ=255)
                                           cttcttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTCAGGAGTCCGCt                                                          >  2:1494732/1‑100 (MQ=255)
                                           cttcttctGAGTGTTTAATTCTGCCCTATTATCCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCt                                                          >  1:2900701/1‑100 (MQ=255)
                                            ttcttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCtt                                                         >  1:415434/1‑100 (MQ=255)
                                             tcttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTa                                                        <  1:4265680/100‑1 (MQ=255)
                                             tcttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTa                                                        <  1:1325251/100‑1 (MQ=255)
                                             tcttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTa                                                        <  2:690265/100‑1 (MQ=255)
                                               ttctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATg                                                      <  1:2347208/100‑1 (MQ=255)
                                                 ctGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGAc                                                    >  2:1347256/1‑100 (MQ=255)
                                                  tGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGa                                                     <  2:3882972/98‑1 (MQ=255)
                                                     gtgtTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGt                                                >  1:3792343/1‑100 (MQ=255)
                                                       gtTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCt                                              >  2:2235802/1‑100 (MQ=255)
                                                        tttAATGCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTg                                             <  1:1694053/100‑1 (MQ=255)
                                                          tAATTCTGCCCTATTATTCTCCGAGCCCTGTAACAAGTTCGGCAGGGCTTTCTTTAGTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTgg                                            >  1:2154831/1‑99 (MQ=38)
                                                             ttCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTc                                        >  1:3266842/1‑100 (MQ=255)
                                                                tGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTcgcg                                     >  2:3808527/1‑100 (MQ=255)
                                                                 gCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCgcgc                                    >  2:4401253/1‑100 (MQ=255)
                                                                  cccTATTATTCTCCGTGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCt                                   >  2:1781436/1‑100 (MQ=255)
                                                                   ccTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCtt                                  >  1:1013342/1‑100 (MQ=255)
                                                                    cTATTATGCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCttn                                 <  2:3461316/100‑2 (MQ=255)
                                                                       ttattCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGgc                              <  2:3682056/100‑1 (MQ=255)
                                                                        tattCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGgcg                             >  1:417193/1‑100 (MQ=255)
                                                                        tattCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGgcg                             >  1:3804533/1‑100 (MQ=255)
                                                                         attCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGgcgc                            <  2:1635409/100‑1 (MQ=255)
                                                                         attCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGgcgc                            >  2:2857729/1‑100 (MQ=255)
                                                                           tcttCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTg                          <  1:4192029/100‑1 (MQ=255)
                                                                           tctcCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTg                          <  2:1247306/100‑1 (MQ=255)
                                                                           tctcCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGGCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTg                          >  2:1309560/1‑100 (MQ=255)
                                                                            ctcCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGc                         >  1:1865821/1‑100 (MQ=255)
                                                                             tcCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCa                        <  2:2774754/100‑1 (MQ=255)
                                                                             agaaaGCCCTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCa                        >  2:937377/5‑100 (MQ=21)
                                                                              gaaaGCCCTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCaa                       >  2:117668/4‑100 (MQ=21)
                                                                              gaaaGCCCTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCaa                       >  2:597294/4‑100 (MQ=21)
                                                                               cGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCaaa                      >  2:4277086/1‑100 (MQ=255)
                                                                                gAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAAt                     <  2:3197573/100‑1 (MQ=255)
                                                                                aaGCCCTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAAt                     >  1:1919785/2‑100 (MQ=21)
                                                                                aaGCCCTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATAAGGAGGAGTCTGCTTATGACCCATCTGACTTCGCGCTTTGGCGCtcc                         >  1:948884/2‑94 (MQ=18)
                                                                                  gCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCt                   >  2:346597/1‑100 (MQ=255)
                                                                                  gCCCTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCt                   >  1:142726/1‑100 (MQ=25)
                                                                                  gCCCTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCt                   >  2:863261/1‑100 (MQ=21)
                                                                                  gCCCTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCt                   >  1:3689134/1‑100 (MQ=21)
                                                                                  gCCCTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGGCTGCTTATGACCCGTCTCTCTTCACGCTTTGGCGCTGCAAATCt                   >  1:3228756/1‑100 (MQ=21)
                                                                                   cccTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCtt                  >  1:2547865/1‑100 (MQ=255)
                                                                                   cccTGCCGCACTTGCGACAGGGCTTTCTTTATTTGTTTTACCCATCCTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCacctct                   >  1:756549/1‑94 (MQ=18)
                                                                                    ccTGTCGCACGTGCGGCAGGGCTTTCTTTCGTTGTTGCACCCATACTGAGGAGTCTGCTCATGCCCCGTCTGGCTTCGCGCTTTGGCGCCGCAAATCttc                 <  1:3505772/100‑2 (MQ=21)
                                                                                    ccTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTa                 >  2:1436023/1‑100 (MQ=255)
                                                                                    ccTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTa                 >  1:1593153/1‑100 (MQ=255)
                                                                                     ctgccgcacttGCGACAGGGCTCTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTAt                >  2:3877154/11‑100 (MQ=21)
                                                                                     cTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTAt                <  2:2285013/100‑1 (MQ=255)
                                                                                       gTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTGACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTAt                <  1:1797178/98‑1 (MQ=255)
                                                                                       gTCGCAAGTGCGGCAGGACTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTAt                >  1:2804981/1‑98 (MQ=255)
                                                                                            aaGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGt         <  2:3588680/100‑1 (MQ=255)
                                                                                             aGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTg        >  1:4326221/1‑100 (MQ=255)
                                                                                               tGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAc      <  1:2422567/100‑1 (MQ=255)
                                                                                                gCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAcc     >  1:3880878/1‑100 (MQ=255)
                                                                                                gCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAcc     <  2:3028549/100‑1 (MQ=255)
                                                                                                  ggCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAccgt   >  2:4013017/1‑100 (MQ=255)
                                                                                                   gCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCCGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGAccgtc  >  1:1270101/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TGCGTTCCTCGGTGGTGTCAGTGGCCCGCGCTGTCTGGGATTTCTTCTTCTGAGTGTTTAATTCTGCCCTATTATTCTCCGAGCCCTGTCGCAAGTGCGGCAGGGCTTTCTTTATTTGTTTTACCCATACTGAGGAGTCTGCTTATGACCCGTCTGGCTTCGCGCTTTGGCGCTGCAAATCTTATTCGCCGTGACCGTC  >  NC_000913/2768521‑2768719

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