Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,820,248 Δ1 bp 6.8% coding (923/1683 nt) ligB ← DNA ligase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,820,2480A. 6.8% 531.8 / 28.8 161coding (923/1683 nt)ligBDNA ligase B
Reads supporting (aligned to +/- strand):  ref base A (85/65);  new base . (4/7);  total (89/72)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

Only 100 of 165 total aligned reads displayed.

CCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCAACTTCA  >  NC_000913/3820099‑3820397
                                                                                                                                                      |                                                                                                                                                     
cccACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTa                                                                                                                                                       >  2:28647/1‑150 (MQ=255)
  cACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACa                                                                                                                                                     >  2:491754/1‑150 (MQ=255)
     tCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGt                                                                                                                                                  >  2:1262171/1‑150 (MQ=255)
         gCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGTAGCTCAGGTTgcc                                                                                                                                              <  1:1195181/150‑3 (MQ=255)
          ccAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGTAGCTCAGGTTgccg                                                                                                                                             <  1:974546/150‑4 (MQ=255)
            aGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGTAGCTCAGGTTgccgaa                                                                                                                                           <  1:1310244/150‑6 (MQ=255)
               gcCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGGGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTcc                                                                                                                                        <  1:574749/150‑1 (MQ=255)
                ccTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCa                                                                                                                                       <  1:276702/150‑1 (MQ=255)
                 cTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAg                                                                                                                                      <  1:491754/150‑1 (MQ=255)
                   gACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGc                                                                                                                                    >  1:236195/1‑150 (MQ=255)
                     cGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGGGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGcca                                                                                                                                  <  1:90783/150‑1 (MQ=255)
                          ccAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGTAGCTCAGGTTgcc                                                                                                                                              >  2:1074661/1‑131 (MQ=255)
                           cAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGcc                                                                                                                            >  2:27934/1‑150 (MQ=255)
                           cAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGcc                                                                                                                            <  2:1201692/150‑1 (MQ=255)
                            aaTATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCa                                                                                                                           <  2:1068230/150‑1 (MQ=255)
                            aaTATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCa                                                                                                                           >  2:438918/1‑150 (MQ=255)
                            aaTATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCa                                                                                                                           >  2:781799/1‑150 (MQ=255)
                             atatTCACCCGCTGGGCTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCAc                                                                                                                          >  2:728050/1‑150 (MQ=255)
                             atatTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCAc                                                                                                                          <  1:10187/150‑1 (MQ=255)
                                  cACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCCCCAGCCACTCTGc                                                                                                                     >  2:1367174/1‑150 (MQ=255)
                                   aCCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGcc                                                                                                                    <  2:1187011/150‑1 (MQ=255)
                                   aCCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGcc                                                                                                                    <  1:1426528/150‑1 (MQ=255)
                                     ccGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTg                                                                                                                  <  1:1173369/150‑1 (MQ=255)
                                      cGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGc                                                                                                                             >  1:486370/1‑138 (MQ=255)
                                            aCTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCActc                                                                                                                        <  2:17477/137‑1 (MQ=255)
                                            aCTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGt                                                                                                           >  2:414566/1‑150 (MQ=255)
                                             cTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTa                                                                                                          >  2:382692/1‑150 (MQ=255)
                                             cTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTa                                                                                                          >  1:742708/1‑150 (MQ=255)
                                             cTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTa                                                                                                          <  2:1030431/150‑1 (MQ=255)
                                              ttttttATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAg                                                                                                         <  1:763227/150‑1 (MQ=255)
                                                ttttATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGcc                                                                                                       >  1:960302/1‑150 (MQ=255)
                                                ttttATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGcc                                                                                                       >  1:855142/1‑150 (MQ=255)
                                                      catcCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGc                                                                                                 >  2:532825/1‑150 (MQ=255)
                                                        tcCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCgg                                                                                               >  1:946965/1‑150 (MQ=255)
                                                           aGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGAt                                                                                            >  1:448887/1‑150 (MQ=255)
                                                           aGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGAt                                                                                            >  2:1211345/1‑150 (MQ=255)
                                                               tGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTACCTGGCCCGGTAGCCAATGGCGGGATTCTg                                                                                        >  2:1109399/1‑150 (MQ=255)
                                                                 aCAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGc                                                                                      >  2:323542/1‑150 (MQ=255)
                                                                 aCAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGc                                                                                      <  2:258218/150‑1 (MQ=255)
                                                                  cAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGct                                                                                     >  2:977833/1‑150 (MQ=255)
                                                                       gCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTc                                                                                >  2:85808/1‑150 (MQ=255)
                                                                          aGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTcgcc                                                                             >  1:1112101/1‑150 (MQ=255)
                                                                           gTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTcgccg                                                                            <  2:589103/150‑1 (MQ=255)
                                                                           gTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTcgccg                                                                            >  1:366370/1‑150 (MQ=255)
                                                                              aCGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCAGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTc                                                                         >  2:1238343/1‑150 (MQ=255)
                                                                               cGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCg                                                                        >  1:1259296/1‑150 (MQ=255)
                                                                                    ccaccGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTATCGCCGCTCGTACaa                                                                   >  2:856141/1‑150 (MQ=255)
                                                                                      accGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACt                                                                 >  1:971378/1‑150 (MQ=255)
                                                                                       ccGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTAc                                                                     >  2:1281443/1‑145 (MQ=255)
                                                                                       ccGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTAc                                                                     <  1:1281443/145‑1 (MQ=255)
                                                                                       ccGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTa                                                                >  2:336932/1‑150 (MQ=255)
                                                                                       ccGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTCCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTa                                                                >  1:1362941/1‑150 (MQ=255)
                                                                                               ttACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCg                                                        >  1:821345/1‑150 (MQ=255)
                                                                                                tACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCgg                                                       <  1:590413/150‑1 (MQ=255)
                                                                                                   cGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGa                                                    <  1:1405188/150‑1 (MQ=255)
                                                                                                    gCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGAc                                                   <  2:414192/150‑1 (MQ=255)
                                                                                                     ctctTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTTCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACg                                                  >  1:114957/1‑150 (MQ=255)
                                                                                                      tctTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGa                                                 >  1:68741/1‑150 (MQ=255)
                                                                                                         tACCCACCGCACACTGATTTGCCTTCACTTCGGCAAGCTTATCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTTCCTGGCCCGTTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAgg                                              <  2:494200/150‑1 (MQ=255)
                                                                                                         tACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAgg                                              <  2:660127/150‑1 (MQ=255)
                                                                                                          aCCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAggg                                             >  2:1299928/1‑150 (MQ=255)
                                                                                                          aCCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAggg                                             <  1:85808/150‑1 (MQ=255)
                                                                                                           cccACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGt                                            <  1:1262171/150‑1 (MQ=255)
                                                                                                           cccACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGt                                            >  1:584392/1‑150 (MQ=255)
                                                                                                              aCCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAAt                                         >  2:179778/1‑150 (MQ=255)
                                                                                                              aCCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAAt                                         >  2:356601/1‑150 (MQ=255)
                                                                                                              aCCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAAt                                         <  1:831133/150‑1 (MQ=255)
                                                                                                                 gCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTTCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTcc                                      <  2:114957/150‑1 (MQ=255)
                                                                                                                 gCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTcc                                      <  1:532825/150‑1 (MQ=255)
                                                                                                                   aaaCTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGc                                    <  1:1199365/150‑1 (MQ=255)
                                                                                                                   aaaCTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGc                                    >  2:924312/1‑150 (MQ=255)
                                                                                                                    aaCTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCt                                   <  1:728050/150‑1 (MQ=255)
                                                                                                                     aCTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCtt                                  >  1:1388010/1‑150 (MQ=255)
                                                                                                                     aCTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCtt                                  <  1:438918/150‑1 (MQ=255)
                                                                                                                      cTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCttt                                 >  2:929578/1‑150 (MQ=255)
                                                                                                                        gAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTcc                               >  1:968248/1‑150 (MQ=255)
                                                                                                                        gAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTcc                               <  2:14595/150‑1 (MQ=255)
                                                                                                                             gCCTTCCCTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCCCCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCccccac                          >  2:206987/1‑150 (MQ=255)
                                                                                                                             gCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCcaccac                          >  1:658772/1‑150 (MQ=255)
                                                                                                                             gCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCcaccac                          <  2:1414460/150‑1 (MQ=255)
                                                                                                                             gCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGAGGGATTCTGGCTCATTCGCCGCTCGTACAACTAAGACATCTGTGACGCAGGGTAATTCCGCTTTccaccaa                          >  2:212819/1‑149 (MQ=255)
                                                                                                                               cTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTc                        >  1:708808/1‑150 (MQ=255)
                                                                                                                                ttCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCa                       >  2:464265/1‑150 (MQ=255)
                                                                                                                                ttCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCa                       >  1:213046/1‑150 (MQ=255)
                                                                                                                                 tCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCAt                      <  2:575198/150‑1 (MQ=255)
                                                                                                                                  cACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCAtt                     >  2:851913/1‑150 (MQ=255)
                                                                                                                                   aCTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTg                    <  1:28647/150‑1 (MQ=255)
                                                                                                                                   aCTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTg                    <  1:1422575/150‑1 (MQ=255)
                                                                                                                                   aCTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTg                    >  1:299426/1‑150 (MQ=255)
                                                                                                                                     ttCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGTCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCg                  >  2:747307/1‑150 (MQ=255)
                                                                                                                                      tCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCAGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGt                 >  2:694232/1‑150 (MQ=255)
                                                                                                                                        ggCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTAc               >  2:590661/1‑150 (MQ=255)
                                                                                                                                        ggCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTAc               >  2:472791/1‑150 (MQ=255)
                                                                                                                                         gCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTAcg              >  2:883984/1‑150 (MQ=255)
                                                                                                                                         gCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTAcg              >  2:1286548/1‑150 (MQ=255)
                                                                                                                                              cTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACgcgcgc         <  1:10646/150‑1 (MQ=255)
                                                                                                                                               tGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCa        <  1:464265/150‑1 (MQ=255)
                                                                                                                                               tGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCa        <  1:132051/150‑1 (MQ=255)
                                                                                                                                                  aGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCTGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCAAc      <  1:1238343/150‑1 (MQ=255)
                                                                                                                                                      aCAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCAACt     <  1:1296358/147‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGAGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCAACTTCA  >  NC_000913/3820099‑3820397

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: