Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,821,244 A→T 5.0% intergenic (‑30/+55) chbC ← / ← chbB N,N'‑diacetylchitobiose‑specific PTS enzyme IIC component/N,N'‑diacetylchitobiose‑specific PTS enzyme IIB component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,821,2440AT 5.0% 327.7 / 4.9 99intergenic (‑30/+55)chbC/chbBN,N'‑diacetylchitobiose‑specific PTS enzyme IIC component/N,N'‑diacetylchitobiose‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base A (45/49);  new base T (3/2);  total (48/51)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

AACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCATCGACTTTGCCATAAAGCAGCGAGTCAATTACTT  >  NC_000913/1821102‑1821390
                                                                                                                                              |                                                                                                                                                  
aaCGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGc                                                                                                                                             <  1:1197948/150‑1 (MQ=255)
 aCGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCg                                                                                                                                            <  2:693952/150‑1 (MQ=255)
  cGGCATTAAGCGAATAAAGCAATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACACTTAAAAACCGCAATTTAAATATTGCgg                                                                                                                                           >  2:21907/1‑150 (MQ=255)
   ggCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGt                                                                                                                                          >  2:85773/1‑150 (MQ=255)
   ggCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGt                                                                                                                                          <  2:250407/150‑1 (MQ=255)
     cATTAAGCGAATAAAGCCATTTTTGATTGCAAAAACGTGGGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTAt                                                                                                                                        <  2:842042/150‑1 (MQ=255)
      aTTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTAtt                                                                                                                                       >  2:1141228/1‑150 (MQ=255)
         aaGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGAt                                                                                                                                    <  1:228226/150‑1 (MQ=255)
          aGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATATTTAAATTGCGGtttttaag                                                                                                                                   <  1:1276980/150‑5 (MQ=255)
           gCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGAttt                                                                                                                                  >  1:1176175/1‑150 (MQ=255)
            cGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAAAAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTa                                                                                                                                 <  2:808641/150‑1 (MQ=255)
              aaTAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATg                                                                                                                               <  2:882641/150‑1 (MQ=255)
               aTAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCGATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGa                                                                                                                              <  2:590765/150‑1 (MQ=255)
               aTAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGa                                                                                                                              <  2:961756/150‑1 (MQ=255)
                tAAATCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGaa                                                                                                                             <  2:367439/150‑1 (MQ=255)
                tAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGaa                                                                                                                             <  1:492806/150‑1 (MQ=255)
                    gCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATaa                                                                                                                         >  2:1428493/1‑150 (MQ=255)
                       aTTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAActc                                                                                                                      <  2:93565/150‑1 (MQ=255)
                         ttttGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCtt                                                                                                                    <  2:992490/150‑1 (MQ=255)
                            tGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGa                                                                                                                 <  1:787295/150‑1 (MQ=255)
                              aTTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACg                                                                                                               >  1:1172196/1‑150 (MQ=255)
                               ttGCATTAACTTGTGGCTGCTTTCCTATTTTAACTGGAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACgg                                                                                                              <  2:1021275/150‑1 (MQ=255)
                               ttGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACgg                                                                                                              >  2:802595/1‑150 (MQ=255)
                                tGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTATTTGACggg                                                                                                             >  2:697440/1‑150 (MQ=255)
                                 gCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGa                                                                                                            >  1:1006906/1‑150 (MQ=255)
                                 gCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGa                                                                                                            >  1:853039/1‑150 (MQ=255)
                                  cATTAACGTGCGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAAATCTTTGACGGGaa                                                                                                           >  2:1600250/1‑150 (MQ=255)
                                   aTTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGaaa                                                                                                          <  1:1509576/150‑1 (MQ=255)
                                    ttAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGaaaa                                                                                                         <  2:592892/150‑1 (MQ=255)
                                    ttAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGaaaa                                                                                                         <  2:775427/150‑1 (MQ=255)
                                      aaCGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAAtt                                                                                                       <  1:342833/150‑1 (MQ=255)
                                         gtgtGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTaa                                                                                                    <  2:681732/150‑1 (MQ=255)
                                          tgtgGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTaaa                                                                                                   <  2:7054/150‑1 (MQ=255)
                                          tgtgGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTaaa                                                                                                   <  1:1176704/150‑1 (MQ=255)
                                            tgGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAaaat                                                                                                 <  2:879796/150‑1 (MQ=255)
                                               ctgctTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATaa                                                                                                                         >  2:911321/1‑123 (MQ=255)
                                               ctgctTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATaa                                                                                                                         <  1:911321/123‑1 (MQ=255)
                                               ctgctTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAaaataaa                                                                                              >  2:767813/1‑150 (MQ=255)
                                               ctgctTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAaaataaa                                                                                              <  1:1332275/150‑1 (MQ=255)
                                                  ctTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAaatta                                                                                           <  1:1268001/150‑1 (MQ=255)
                                                  ctTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAATAaattaa                                                                                          >  1:902231/1‑150 (MQ=255)
                                                   tttCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAaattaa                                                                                          >  2:1437883/1‑150 (MQ=255)
                                                    ttCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAaattaat                                                                                         >  1:354397/1‑150 (MQ=255)
                                                       cTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTtg                                                                                      <  1:761247/150‑1 (MQ=255)
                                                                cTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCtt                                                                             >  1:1561573/1‑150 (MQ=255)
                                                                cTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCtt                                                                             >  1:84105/1‑150 (MQ=255)
                                                                cTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCtt                                                                             >  1:139287/1‑150 (MQ=255)
                                                                   cAAAAGGGAGGACTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAATGGCAATTTAAATATTGGGGTATTTATTTATGAAATAACTCTTTTACGTGAAAATTTAAAATAATTTAATTTGCTGCGGCttttt                                                                          <  2:498469/150‑1 (MQ=255)
                                                                       aGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAAt                                                                      >  2:945064/1‑150 (MQ=255)
                                                                         ggaggaGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCg                                                                    <  1:85773/150‑1 (MQ=255)
                                                                          gaggagTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGc                                                                   >  1:532776/1‑150 (MQ=255)
                                                                           aggagTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCt                                                                  <  1:210259/150‑1 (MQ=255)
                                                                             gagTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGc                                                                <  2:912233/150‑1 (MQ=255)
                                                                              agTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCg                                                                                 <  1:61765/132‑1 (MQ=255)
                                                                              agTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCg                                                                                 >  2:61765/1‑132 (MQ=255)
                                                                              agTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCa                                                               >  1:1600565/1‑150 (MQ=255)
                                                                              agTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCa                                                               <  2:1077681/150‑1 (MQ=255)
                                                                              agTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTCCa                                                               >  2:1045765/1‑150 (MQ=255)
                                                                                tACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAAc                                                             >  1:645040/1‑150 (MQ=255)
                                                                                tACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAAc                                                             >  1:554195/1‑150 (MQ=255)
                                                                                 aCCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAAcc                                                            >  2:677800/1‑150 (MQ=255)
                                                                                 aCCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAAcc                                                            <  1:130565/150‑1 (MQ=255)
                                                                                   cTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGc                                                          <  1:884380/150‑1 (MQ=255)
                                                                                      ttCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTATAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGc                                                       <  2:1480423/150‑1 (MQ=255)
                                                                                        cAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAAttt                                                                                       >  1:440510/1‑116 (MQ=255)
                                                                                        cAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAAttt                                                                                       <  2:440510/116‑1 (MQ=255)
                                                                                        cAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCt                                                     >  2:459238/1‑150 (MQ=255)
                                                                                         aaGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCtt                                                    >  1:7498/1‑150 (MQ=255)
                                                                                           gCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTaa                                                  >  2:1576473/1‑150 (MQ=255)
                                                                                           gCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTaa                                                  <  2:307772/150‑1 (MQ=255)
                                                                                           gCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTaa                                                  >  1:1222565/1‑150 (MQ=255)
                                                                                           gCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTaa                                                  >  1:857466/1‑150 (MQ=255)
                                                                                            cGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAg                                                 <  1:886740/150‑1 (MQ=255)
                                                                                             gATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGc                                                <  1:1428493/150‑1 (MQ=255)
                                                                                             gATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGc                                                <  2:354397/150‑1 (MQ=255)
                                                                                                 cAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAATTTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGc                                            >  1:767832/1‑150 (MQ=255)
                                                                                                 cAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGc                                            >  2:81929/1‑150 (MQ=255)
                                                                                                    tAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTa                                         <  2:139287/150‑1 (MQ=255)
                                                                                                      aCATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTaaa                                       >  1:394645/1‑150 (MQ=255)
                                                                                                        aTTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAAcc                                     <  1:802595/150‑1 (MQ=255)
                                                                                                         ttACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCa                                    >  2:516110/1‑150 (MQ=255)
                                                                                                             tCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCATCGa                                >  2:374472/1‑150 (MQ=255)
                                                                                                             aaaaataaaaaTACCCTAAAAAACCGAAATTAAAATATTCCGTTATTTATTAATCAAAAAACCCTTTCACGGCAAAATTAAAAAAAAATAAATTTGCTGCGCCTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCATCGAc                               <  1:1282290/144‑1 (MQ=255)
                                                                                                              cATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCATCGAc                               >  2:271861/1‑150 (MQ=255)
                                                                                                              cATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCATCGAc                               <  1:942365/150‑1 (MQ=255)
                                                                                                               aTAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCATCGACt                              >  2:843825/1‑150 (MQ=255)
                                                                                                               aTAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCATCGACt                              >  2:1513563/1‑150 (MQ=255)
                                                                                                                   aaaaaTACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCATCGACTTTGc                          >  2:358027/1‑150 (MQ=255)
                                                                                                                        tACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCATCGACTTTGCCATaa                     >  2:264403/1‑150 (MQ=255)
                                                                                                                          cccTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCATCGACTTTGCCATAAag                   <  1:677800/150‑1 (MQ=255)
                                                                                                                           ccTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCATCGACTTTGCCATAAagc                  >  1:717853/1‑150 (MQ=255)
                                                                                                                                 ataCCGCAATATTTAAATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCATCGACTTTGc                          >  2:298007/3‑136 (MQ=255)
                                                                                                                                 ataCCGCAATATTTAAATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCATCGACTTTGc                          <  1:298007/134‑1 (MQ=255)
                                                                                                                                  aaCCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCATCGACTTTGCCATAAAGCAGCGAGt           <  1:1447359/150‑1 (MQ=255)
                                                                                                                                   aCCGCAATATTTAAATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCATCGACTTTGCCATAAAGCAGCGAGTc          >  1:67580/1‑150 (MQ=255)
                                                                                                                                    ccGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCATCGACTTTGCCATAAAGCAGCGAGTCa         >  2:830740/1‑150 (MQ=255)
                                                                                                                                    ccGCAATATTTAAATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCATCGACTTTGCCATAAAGCAGCGAGTCa         >  2:1089662/1‑150 (MQ=255)
                                                                                                                                         aTTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCATCGACTTTGCCATAAAGCAGCGAGTCAATTAc    <  1:459238/150‑1 (MQ=255)
                                                                                                                                           ttAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCATCGACTTTGCCATAAAGCAGCGAGTCAATTACtt  <  1:21907/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
AACGGCATTAAGCGAATAAAGCCATTTTTGATTGCATTAACGTGTGGCTGCTTTCCTATTTTAACTGCAAAAGGGAGGAGTACCTTTTCAAGCGATGCAATAACATTACTCATAGAAAAATACCCTTAAAAACCGCAATTTAAATATTGCGGTATTGATTTATGAAATAACTCTTTGACGGGAAAATTTAAAATAAATTAATTTGCTGCGGCTTTTTTAATCGCTGCAACCGCAGCCTTAAGCACGCCTAAACCATCGACTTTGCCATAAAGCAGCGAGTCAATTACTT  >  NC_000913/1821102‑1821390

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