Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 669,244 T→C 26.2% intergenic (‑208/+52) rsfS ← / ← cobC ribosomal silencing factor/putative alpha‑ribazole‑5'‑P phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913669,2440TC26.2% 182.2 / 76.1 120intergenic (‑208/+52)rsfS/cobCribosomal silencing factor/putative alpha‑ribazole‑5'‑P phosphatase
Reads supporting (aligned to +/- strand):  ref base T (38/48);  new base C (18/13);  total (57/63)
Fisher's exact test for biased strand distribution p-value = 2.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.16e-01

Only 100 of 124 total aligned reads displayed.

CATCCCACCCCGCGCCGCAGATTGCATCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGT  >  NC_000913/669145‑669343
                                                                                                   |                                                                                                   
catcCCACCCCGCGCCGCAGATTGCATCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATAGCAGCCt                                                                                                     >  1:99485/1‑100 (MQ=255)
 atcCCACCCCGCGCCGCAGATTGCATCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTg                                                                                                    >  2:4211873/1‑100 (MQ=255)
  tcCCACCCCGCGCCGCAGATTGCATCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGc                                                                                                   <  2:2628314/100‑1 (MQ=255)
  tcCCACCCCGCGCCGCAGATTGCATCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGc                                                                                                   <  2:818216/100‑1 (MQ=255)
     cACCCCGCGCCGCAGATTGCATCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCg                                                                                                >  1:3125624/1‑100 (MQ=255)
      aCCCCGCGCCGCAGATTGCATCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCgg                                                                                               >  2:1538467/1‑100 (MQ=255)
       ccccGCGCCGCGGATTGCATCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCggg                                                                                              <  1:1116152/100‑1 (MQ=255)
       ccccGCGCCGCAGATTGCATCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCggg                                                                                              >  2:3204289/1‑100 (MQ=255)
        cccGCGCCGCAGATTGCATCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGt                                                                                             >  2:121753/1‑100 (MQ=255)
         ccGCGCCGCAGATTGCATCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTc                                                                                            <  2:4045691/100‑1 (MQ=255)
         ccGCGCCGCAGATCGCATCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTc                                                                                            <  2:4074738/100‑1 (MQ=255)
           gcgcCGCAGATTGCATCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTg                                                                                          <  1:2283455/100‑1 (MQ=255)
            cgccgcAGATTGCATCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGt                                                                                         <  2:610817/100‑1 (MQ=255)
             gccgcAGATTGCATCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTc                                                                                        >  2:4347983/1‑100 (MQ=255)
             gccgcAGATTGCATCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTc                                                                                        <  2:1120904/100‑1 (MQ=255)
              ccgcAGATTGCATCGCCATTTTTGTTAAACACTTGCTACAAAGTCGTGTCTGGTGGAAAAAGTCTGGCTGCGTAGAATATCAGCCTGCCCGGGTCTGTc                                                                                        <  1:718307/99‑1 (MQ=255)
               cgcAGATTGCATCGGCGTTTTTGGAAAGCAATTTGTACAATGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTc                                                                                        <  1:4017209/98‑1 (MQ=255)
                gcAGATTGCATCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGGGGCAAAAGTCTGGCTGCGGAGAATATCAGCTCGCCCGGGTCTGTCAAt                                                                                     >  1:1803307/1‑100 (MQ=255)
                gcAGATTGCATCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTc                                                                                        <  1:3534102/97‑1 (MQ=255)
                   gATTGCATCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGt                                                                                  >  1:1042496/1‑100 (MQ=255)
                     ttGCATCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCg                                                                                <  2:371641/100‑1 (MQ=255)
                      tgaagcGCCATTTTTGTAAAACAATTTTTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCCGCCCGGGTCTGTCAATGGTCg                                                                                <  1:1093502/94‑1 (MQ=255)
                      tGCATCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCCGCCCGGGTCTGTCAATGGTCGt                                                                               >  1:3301858/1‑100 (MQ=255)
                        cATCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCCGCCCGGGTCTGTCAATGGTCGTTg                                                                             >  2:3516873/1‑100 (MQ=255)
                         aTCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGCGAATATCAGCCTGCCCCGGTCTGTCAATGGTCGTTGc                                                                            >  1:3742389/1‑100 (MQ=255)
                          tCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCg                                                                           <  1:3111068/100‑1 (MQ=255)
                           cGCCATTTTTGTGAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCgg                                                                          <  1:626799/100‑1 (MQ=255)
                            gCCATTTTTGTAAAACAATTTCTACACAGTCGTGTCTGGCGGAAAAAGTCTGCCTGCGGAGAATACCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGc                                                                         >  1:1713729/1‑100 (MQ=255)
                            gCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGc                                                                         <  1:433065/100‑1 (MQ=255)
                                 ttttGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCCGCCCGGGTCTGTCAATGGTCGTTGCGGCTTAcc                                                                    >  2:3220728/1‑100 (MQ=255)
                                 ttttGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCCGCCCGGGTCTGTCAATGGTCGTTGCGGCTTAca                                                                    <  1:2588363/100‑2 (MQ=255)
                                 ttttGTAAAACAATTTCTACAAAGTCAAAAATAAAAGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTAcc                                                                    <  1:3591709/100‑1 (MQ=255)
                                   ttGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCg                                                                  <  2:3273179/100‑1 (MQ=255)
                                   ttGTAAAACAATTTCTACAAAGTCGTGTCTGGCGAAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCg                                                                  <  2:2648815/100‑1 (MQ=255)
                                     gTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCCGCCCGGGTCTGTCAATGGTCGTTGCGGCTTAccagat                                                                >  2:1742920/1‑96 (MQ=255)
                                       aaaaCAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCCGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTaaa                                                              >  1:2152270/1‑100 (MQ=255)
                                         aaCAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTaaaaa                                                            >  2:3866217/1‑100 (MQ=255)
                                          aCAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTaaaaaa                                                           <  2:2776671/100‑1 (MQ=255)
                                          aCAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTaaaaaa                                                           >  1:2991989/1‑100 (MQ=255)
                                            aaTTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAAc                                                         <  2:803631/100‑1 (MQ=255)
                                             aTTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCCGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACa                                                        >  1:1870118/1‑100 (MQ=255)
                                              tttCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACTCGTAAAAAAACAg                                                       <  1:2727592/100‑1 (MQ=255)
                                              ttccTACAAAGTAGTGTCTGCCGGAAAACGTCTGGCTGCGGAGAATCCCAGCGTTACCGGGTCGGCCAATGGTCGTTGCGGCATACCCGTCAACAAACAg                                                       >  1:1525602/4‑100 (MQ=255)
                                               ttCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGa                                                      <  2:1350253/100‑1 (MQ=255)
                                                 cTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGaaa                                                    >  1:1407317/1‑100 (MQ=255)
                                                  tACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCCGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGaaaa                                                   >  2:2438390/1‑100 (MQ=255)
                                                   aCAAAGTCGTGTCTGGTGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAg                                                  >  2:3900778/1‑100 (MQ=255)
                                                        gTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATg                                             <  2:2274506/100‑1 (MQ=255)
                                                        gTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATg                                             >  2:4314125/1‑100 (MQ=255)
                                                        gTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATg                                             <  1:3204271/100‑1 (MQ=255)
                                                         tCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCCGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGc                                            <  2:3516851/100‑1 (MQ=255)
                                                          cGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCa                                           >  2:766101/1‑100 (MQ=255)
                                                            tgtCTGGCGGAAAATGTCTGGCTGGGGAGAATATCAGCCGGCGCGGTTCTGTCAGTGGTCGTTTCTGCTTAACCGTAAAAAAACAGAAAAGTCa                                               <  1:1888638/94‑1 (MQ=255)
                                                            tgtCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCAtt                                         >  1:1650285/1‑100 (MQ=255)
                                                             gtCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCCGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATtc                                        <  2:794186/100‑1 (MQ=255)
                                                               cTGGCGGAAAAAGTCTGGATGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATtctc                                      <  1:2026472/100‑1 (MQ=255)
                                                                 ggCGGAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGAc                                   <  2:2683676/100‑1 (MQ=255)
                                                                  gCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGAc                                   >  1:3161212/1‑100 (MQ=255)
                                                                    ggAAAAAGTCTGGCTTCGGACAATATCAACCTGCACAGGTCTGTCAATTGTCGTTGCTTAATAAAAAAAAAAAAACAGAAAAGTCATGCATTCTCGAccc                                 <  1:787223/100‑1 (MQ=255)
                                                                    ggAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGAccc                                 >  1:680357/1‑100 (MQ=255)
                                                                    ggAAAAAGTCTGGCTGCGGAGAATATCAGCCCGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGAccc                                 <  2:2152277/100‑1 (MQ=255)
                                                                      aaaaaGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCg                               >  2:3631448/1‑100 (MQ=255)
                                                                      aaaaaGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCg                               >  2:194115/1‑100 (MQ=255)
                                                                        aaaGTCTGGCTGCGGAGAATATCAGCCCGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGAt                             <  1:19291/100‑1 (MQ=255)
                                                                         aaGTCTGGCTGCGGAGAATATCAGCCCGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATg                            <  1:2438379/100‑1 (MQ=255)
                                                                          aGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATgg                           >  2:1799433/1‑100 (MQ=255)
                                                                          aGTCTGGCTGCGGAGAATATCAGCCGGCCCTGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATgg                           <  1:1413506/100‑1 (MQ=255)
                                                                          aGTCTGGCTGCGGAGAATATCAGCCCGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATgg                           >  1:3234391/1‑100 (MQ=255)
                                                                           gTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGc                          >  2:4413433/1‑100 (MQ=255)
                                                                           gTCTGGCTGCGGAGAATATCAGCCCGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGc                          <  2:1472453/100‑1 (MQ=255)
                                                                               ggCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACgg                      <  2:2025907/100‑1 (MQ=255)
                                                                               ggCTGCGGAGAATATCAGCCCGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACgg                      >  2:2370815/1‑100 (MQ=255)
                                                                                gCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGc                     <  2:458326/100‑1 (MQ=255)
                                                                                 cTGCGGAGAATATCAGCCCGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCt                    >  2:3864444/1‑100 (MQ=255)
                                                                                  tGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTa                   >  2:1106331/1‑100 (MQ=255)
                                                                                  tGCGGAGAATATCAGCCCGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCt                    <  1:4084756/99‑1 (MQ=255)
                                                                                    cGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTAtt                 >  2:68561/1‑100 (MQ=255)
                                                                                     ggAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTg                <  1:1325510/100‑1 (MQ=255)
                                                                                     ggAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTg                <  2:3639235/100‑1 (MQ=255)
                                                                                        gaATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTgtac             >  1:2035734/1‑97 (MQ=255)
                                                                                        gaATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAgg             <  1:1737624/100‑1 (MQ=255)
                                                                                        gaATATCAGCCCGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGatg             <  2:2980356/100‑3 (MQ=255)
                                                                                         aaTATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGa            <  1:431514/100‑1 (MQ=255)
                                                                                          atatCAGCCCGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGAc           >  2:3391211/1‑100 (MQ=255)
                                                                                           tatCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGAcg          <  2:159319/100‑1 (MQ=255)
                                                                                           tatCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGAcg          >  1:4192988/1‑100 (MQ=255)
                                                                                            atCAGCCCGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGAcgc         >  2:3504897/1‑100 (MQ=255)
                                                                                            atCAGCCCGCCCGGGTCTGTCAATGGTCGTTGCGGCGGACACGTAAAAACGCAGAAAAGTCGTGCATTCGCGACTCCTAGGGCACGGCTATTg                >  1:1528935/1‑93 (MQ=255)
                                                                                             tcggCCGGCCTGGGTCTGTCATTGGTCGTTGCGGCTTAGCCGTAAAAAAACAGAAAAGTCATGCATTGTCGACCTTGGTGGCACGGCTATTGAGGAcgc         >  1:3512772/4‑99 (MQ=255)
                                                                                             tCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGAcgcg        <  2:1758950/100‑1 (MQ=255)
                                                                                             tCAGCCTGCCCGGGTCTGTCAAGGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGAcgcg        <  1:1149731/100‑1 (MQ=255)
                                                                                             tCAGCCCGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGAcgcg        >  2:1224010/1‑100 (MQ=255)
                                                                                             tCAGCCCGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGAcgcg        >  2:819082/1‑100 (MQ=255)
                                                                                                gcccgcccGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGt       >  1:1726816/1‑98 (MQ=255)
                                                                                                gcccgcccGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAg     >  1:3682453/1‑100 (MQ=255)
                                                                                                gCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTa      >  2:1666426/1‑99 (MQ=255)
                                                                                                gCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAg     >  2:729479/1‑100 (MQ=255)
                                                                                                 ccTCCCCGGGTCTGTCAAGGGTCGTTGCGGCTTACCAATAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAg     <  1:4313338/99‑1 (MQ=255)
                                                                                                  cTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCAGGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCg   <  1:4327748/100‑1 (MQ=255)
                                                                                                   tGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGt  >  2:4255602/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CATCCCACCCCGCGCCGCAGATTGCATCGCCATTTTTGTAAAACAATTTCTACAAAGTCGTGTCTGGCGGAAAAAGTCTGGCTGCGGAGAATATCAGCCTGCCCGGGTCTGTCAATGGTCGTTGCGGCTTACCCGTAAAAAAACAGAAAAGTCATGCATTCTCGACCCCGATGGCACGGCTATTGAGGACGCGTAGCGT  >  NC_000913/669145‑669343

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: