Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,820,248 Δ1 bp 5.5% 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. 5.5% 304.4 / 9.6 91coding (923/1683 nt)ligBDNA ligase B
Reads supporting (aligned to +/- strand):  ref base A (45/41);  new base . (3/2);  total (48/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

CCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCAACTTCA  >  NC_000913/3820099‑3820397
                                                                                                                                                      |                                                                                                                                                     
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:240484/150‑5 (MQ=255)
           cAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTGTAGCT‑CAGGTTgccga                                                                                                                                            <  1:378371/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:424724/150‑1 (MQ=255)
                           cAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGcc                                                                                                                            >  2:6638/1‑150 (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                                                                                                     <  2:782936/150‑1 (MQ=255)
                                                  ttATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCaa                                                                                                     >  1:518496/1‑150 (MQ=255)
                                                        tcCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCgg                                                                                               >  1:1566228/1‑150 (MQ=255)
                                                        tcCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCgg                                                                                               >  1:1560414/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                                                                      >  2:285522/1‑150 (MQ=255)
                                                                                 cAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTa                                                                      >  1:758371/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:130877/1‑150 (MQ=255)
                                                                                     caccGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAAc                                                                  >  1:136933/1‑150 (MQ=255)
                                                                                     caccGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAAc                                                                  >  1:481921/1‑150 (MQ=255)
                                                                                      accGATATTTTACCGCTCTTACCCACCGTAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTcg                                                                               <  1:887254/136‑1 (MQ=255)
                                                                                      accGATATTTTACCGCTCTTACCCACCGTAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTcg                                                                               >  2:887254/1‑136 (MQ=255)
                                                                                          atatTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGc                                                             <  2:321056/150‑1 (MQ=255)
                                                                                          atatTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGc                                                             <  1:101001/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:244717/150‑1 (MQ=255)
                                                                                                          aCCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAggg                                             >  1:27099/1‑150 (MQ=255)
                                                                                                           cccACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGt                                            >  2:495465/1‑150 (MQ=255)
                                                                                                           cccACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGt                                            >  1:1134391/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                 >  2:695394/5‑150 (MQ=255)
                                                                                                                                      atatcAACCTGTAGCT‑CAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGt                 >  1:254781/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)
                                                                                                                                                  cgCTACAGGTTGATATTTCCCGGCCACCCGCCACTCTGCCTGGCCCGGTAGCCCATGGCGGGATTCGGGCTCTTTCGCCGCTCGTACAACTCCGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCAAc      <  2:277952/149‑1 (MQ=255)
                                                                                                                                                    cTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCAACtt    >  2:950975/1‑150 (MQ=255)
                                                                                                                                                      aCAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCAACTtca  >  1:1182552/1‑150 (MQ=255)
                                                                                                                                                      aCAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCAACTtca  >  2:20223/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CCCACTCCTGCCAGCGCCTGACGGAACCAATATTCACCCGCTGGACTTTTTTATCATCCAGCATGACAGGTGCGAGTGACGCAACCACCGATATTTTACCGCTCTTACCCACCGCAAACTGAATTGCCTTCACTTCGGCAACCTG‑AGCTACAGGTTGATATTTCCAGGCCACCAGCCACTCTGCCTGGCCCGGTAGCCAATGGCGGGATTCTGGCTCTTTCGCCGCTCGTACAACTACGCCATCGGTGACGAAGGGTAATTCCGCTTTCCACCACTCATTGCGTACGCGCGCAACTTCA  >  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: ‑