Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,820,243:1 +T 5.6% coding (928/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,2431.T 5.6% 359.3 / 16.1 90coding (928/1683 nt)ligBDNA ligase B
Reads supporting (aligned to +/- strand):  ref base . (45/40);  new base T (3/2);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

AATATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCA  >  NC_000913/3820094‑3820391
                                                                                                                                                      |                                                                                                                                                    
aaTATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTg                                                                                                                                                       <  2:1070902/150‑1 (MQ=255)
aaTATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTg                                                                                                                                                       >  2:61664/1‑150 (MQ=255)
     cccACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTa                                                                                                                                                 >  2:970637/1‑150 (MQ=255)
         ctacTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAgg                                                                                                                                             >  1:424146/4‑150 (MQ=255)
         cTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAgg                                                                                                                                             >  1:68503/1‑150 (MQ=255)
             tGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGa                                                                                                                                         <  1:1253368/150‑1 (MQ=255)
                cAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGTAGCT‑CAGGTTgccga                                                                                                                                      <  1:378371/150‑5 (MQ=255)
                cAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGTAGCT‑CAGGTTgccga                                                                                                                                      <  1:240484/150‑5 (MQ=255)
                cAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGatat                                                                                                                                      >  2:1058316/1‑150 (MQ=255)
                       tGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAgg                                                                                                                               <  1:143726/150‑1 (MQ=255)
                         aCGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGcc                                                                                                                             <  2:687087/150‑1 (MQ=255)
                         aCGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGcc                                                                                                                             <  2:658277/150‑1 (MQ=255)
                          cGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGcca                                                                                                                            <  2:88255/150‑1 (MQ=255)
                           ggAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGccac                                                                                                                           >  2:766730/1‑150 (MQ=255)
                            gAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGccacc                                                                                                                          <  1:602317/150‑1 (MQ=255)
                                cAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGcc                                                                                                                      >  2:6638/1‑150 (MQ=255)
                                cAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGcc                                                                                                                      <  2:424724/150‑1 (MQ=255)
                                      tCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTg                                                                                                                <  1:61664/150‑1 (MQ=255)
                                       cACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCATCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGc                                                                                                               >  2:351634/1‑150 (MQ=255)
                                           cGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTgg                                                                                                           >  2:101001/1‑150 (MQ=255)
                                              tGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGccc                                                                                                        <  2:837590/150‑1 (MQ=255)
                                               ggACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCg                                                                                                       <  1:766730/150‑1 (MQ=255)
                                               ggACTTTTTTATCATCCAGCATGACAGGAGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACACACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTACAGGCCACCAGCCACTCTGCCTGGCCCg                                                                                                       >  2:1144094/1‑150 (MQ=255)
                                                gACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCgg                                                                                                      <  2:953018/150‑1 (MQ=255)
                                                  cTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTa                                                                                                    <  1:1531558/150‑1 (MQ=255)
                                                   ttttttATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAg                                                                                                   <  2:1152314/150‑1 (MQ=255)
                                                   ttttttATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAg                                                                                                   >  1:782936/1‑150 (MQ=255)
                                                   ttttttATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTAATTTCTGCAACCTG‑TGCTACAGGTTGATATTTCCAGGCCACAATCCACTCTGCCTGGCGAGGTAg                                                                                                   <  2:1263459/150‑1 (MQ=255)
                                                    tttttATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGTTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTa                                                                                                    <  2:1524789/148‑1 (MQ=255)
                                                    tttttATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTa                                                                                                    >  1:1524789/1‑148 (MQ=255)
                                                      tttATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCa                                                                                                <  2:326389/150‑1 (MQ=255)
                                                       ttATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCaa                                                                                               >  1:518496/1‑150 (MQ=255)
                                                       ttATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCaa                                                                                               <  2:782936/150‑1 (MQ=255)
                                                             tcCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCgg                                                                                         >  1:1560414/1‑150 (MQ=255)
                                                             tcCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCgg                                                                                         >  1:1566228/1‑150 (MQ=255)
                                                                         ggTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGctct                                                                             >  1:978188/1‑150 (MQ=255)
                                                                          gTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTAACCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCAACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCtt                                                                            >  1:1095667/1‑150 (MQ=255)
                                                                             ggAGTGACGCAACCACCGATATTTTACCGCTCTTACCCAGCGCAAACTGAATTGCCTTCACTTCGGCAACGTG‑AGCTACAGGTTGATATTTCGAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTcg                                                                         <  1:1421925/149‑1 (MQ=255)
                                                                                  gACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGCATTCTGGCTCTTTCGCCGCt                                                                    >  1:897053/1‑150 (MQ=255)
                                                                                   aCGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTc                                                                   >  1:1588092/1‑150 (MQ=255)
                                                                                    cGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCg                                                                  <  1:6638/150‑1 (MQ=255)
                                                                                      cAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTa                                                                >  1:758371/1‑150 (MQ=255)
                                                                                      cAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTa                                                                >  2:285522/1‑150 (MQ=255)
                                                                                       aaCCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTAc                                                               <  2:607498/150‑1 (MQ=255)
                                                                                       aaCCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTAc                                                               >  2:245104/1‑150 (MQ=255)
                                                                                        accaccGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACa                                                              >  1:1036123/1‑150 (MQ=255)
                                                                                        accaccGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACa                                                              <  1:1668461/150‑1 (MQ=255)
                                                                                         ccaccGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACaa                                                             >  2:195032/1‑150 (MQ=255)
                                                                                          caccGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAAc                                                            >  1:136933/1‑150 (MQ=255)
                                                                                          caccGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAAc                                                            >  1:481921/1‑150 (MQ=255)
                                                                                          caccGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAAc                                                            >  1:130877/1‑150 (MQ=255)
                                                                                           accGATATTTTACCGCTCTTACCCACCGTAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTcg                                                                         >  2:887254/1‑136 (MQ=255)
                                                                                           accGATATTTTACCGCTCTTACCCACCGTAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTcg                                                                         <  1:887254/136‑1 (MQ=255)
                                                                                               atatTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGc                                                       <  1:101001/150‑1 (MQ=255)
                                                                                               atatTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGc                                                       <  2:321056/150‑1 (MQ=255)
                                                                                                tatTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTAAAACTACGcc                                                      >  2:102984/1‑150 (MQ=255)
                                                                                                  ttttACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCaa                                                    >  1:1044776/1‑149 (MQ=255)
                                                                                                   tttACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTa                                                                                                    >  2:510873/1‑101 (MQ=255)
                                                                                                   tttACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTa                                                                                                    <  1:510873/101‑1 (MQ=255)
                                                                                                        cGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGa                                              <  2:518496/150‑1 (MQ=255)
                                                                                                              tACCCACCGCAATCTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGTCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCAGCTCGTACTACTACGCCATCTTTGACGAAgg                                        <  2:221750/150‑1 (MQ=255)
                                                                                                               aCCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAggg                                       >  1:27099/1‑150 (MQ=255)
                                                                                                               aCCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAggg                                       <  1:244717/150‑1 (MQ=255)
                                                                                                                cccACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGt                                      >  1:1134391/1‑150 (MQ=255)
                                                                                                                cccACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGt                                      >  2:495465/1‑150 (MQ=255)
                                                                                                                cccACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTGGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGt                                      <  2:1312348/150‑1 (MQ=255)
                                                                                                                  cACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTaa                                    >  2:1499901/1‑150 (MQ=255)
                                                                                                                   aCCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAAt                                   >  2:1263150/1‑150 (MQ=255)
                                                                                                                     cGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGGTATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTc                                 >  1:468063/1‑150 (MQ=255)
                                                                                                                      gCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTcc                                >  1:1253318/1‑150 (MQ=255)
                                                                                                                       cAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCg                               <  2:442114/150‑1 (MQ=255)
                                                                                                                        aaaCTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGc                              >  1:1231752/1‑150 (MQ=255)
                                                                                                                         aaCTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGt                                                                 <  2:870271/114‑1 (MQ=255)
                                                                                                                         aaCTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGt                                                                 >  1:870271/1‑114 (MQ=255)
                                                                                                                             gAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTcc                         <  1:1058316/150‑1 (MQ=255)
                                                                                                                                ttGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTccacc                      <  2:1036123/150‑1 (MQ=255)
                                                                                                                                 tGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTccacca                     <  2:424146/150‑1 (MQ=255)
                                                                                                                                  gCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTcc                         <  1:1339557/145‑1 (MQ=255)
                                                                                                                                  gCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTcc                         >  2:1339557/1‑145 (MQ=255)
                                                                                                                                   ccTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACt                   <  1:191254/150‑1 (MQ=255)
                                                                                                                                        aCTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTg              >  2:486043/1‑150 (MQ=255)
                                                                                                                                         cTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGc             >  1:1353095/1‑150 (MQ=255)
                                                                                                                                         cTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGc             >  1:1229105/1‑150 (MQ=255)
                                                                                                                                           tCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGt           >  2:627023/1‑150 (MQ=255)
                                                                                                                                           atatcAACCTGTAGCT‑CAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGt           >  1:254781/5‑150 (MQ=255)
                                                                                                                                           atatcAACCTGTAGCT‑CAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGt           >  2:695394/5‑150 (MQ=255)
                                                                                                                                            tatcAACCTGTAGCT‑CAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTa          >  1:550649/4‑150 (MQ=255)
                                                                                                                                              gCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTAcg        <  1:458268/150‑1 (MQ=255)
                                                                                                                                              gCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTAcg        <  2:1164819/150‑1 (MQ=255)
                                                                                                                                                  ccTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTAcgcgcg    >  2:376943/1‑150 (MQ=255)
                                                                                                                                                  ccTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTAcgcgcg    <  1:285522/150‑1 (MQ=255)
                                                                                                                                                    tG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCa  <  1:351634/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
AATATCCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCA  >  NC_000913/3820094‑3820391

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: ‑