Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 781,672 A→G 5.2% intergenic (+20/‑413) lysQ → / → nadA tRNA‑Lys/quinolinate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913781,6720AG 5.2% 207.6 / 2.5 77intergenic (+20/‑413)lysQ/nadAtRNA‑Lys/quinolinate synthase
Reads supporting (aligned to +/- strand):  ref base A (40/33);  new base G (2/2);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

CCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAA  >  NC_000913/781528‑781817
                                                                                                                                                |                                                                                                                                                 
ccTGCACGACCCACCATACTGAATGATAAAGGAAACATAAACAAGCAAGGGGTCGTTAGCTAAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAAACCTGAACGACACACCAATGTAAAAAAGCGCCCTAAAGgcgc                                                                                                                                              >  2:1002738/1‑150 (MQ=17)
ccTGCACGACCCACCAACCTGAATGATTAAGCCAGCATAAACCCGAAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGCTCGCAGGTTCGAATCATCAACGACCCACCAATGTAAAAAAGCGCCCTAAAGgcgc                                                                                                                                              >  2:486910/1‑150 (MQ=18)
   gcgcGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCttt                                                                                                                                           >  2:1798814/4‑150 (MQ=255)
   gCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCttt                                                                                                                                           <  2:605653/150‑1 (MQ=255)
    cACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttt                                                                                                                                          >  1:740504/1‑150 (MQ=255)
       gACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCtttttta                                                                                                                                       <  2:406068/150‑2 (MQ=255)
         cccACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttttt                                                                                                                                        >  2:1364194/1‑147 (MQ=255)
         cccACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttttt                                                                                                                                        <  1:1364194/147‑1 (MQ=255)
            accaTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATc                                                                                                                                  <  1:167027/150‑1 (MQ=255)
              caTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTg                                                                                                                                <  1:1506539/150‑1 (MQ=255)
               aTCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGc                                                                                                                               <  2:1908335/150‑1 (MQ=255)
                  cTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGAt                                                                                                                            >  1:264481/1‑150 (MQ=255)
                   tGAATTATTAAGGCACCCTAATCCCGCAAGCGGTCGATATCTCTTCTGTTATAGCATTTTACATTTTATCAATTTGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATa                                                                                                                           <  2:1518820/150‑1 (MQ=14)
                          ttAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAg                                                                                                                    <  1:365352/150‑1 (MQ=255)
                              ggCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTc                                                                                                                <  1:1759982/150‑1 (MQ=255)
                              ggCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTc                                                                                                                <  2:1900705/150‑1 (MQ=255)
                                cagcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGa                                                                                                              <  1:1885578/150‑1 (MQ=255)
                                   caTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGACGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAAcc                                                                                                           >  2:1298127/1‑150 (MQ=255)
                                     tAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTg                                                                                                         >  1:1996352/1‑150 (MQ=255)
                                      aaTCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTgc                                                                                                        <  2:982891/150‑1 (MQ=40)
                                      aaTCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTgc                                                                                                        <  2:975739/150‑1 (MQ=40)
                                       aTCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTgca                                                                                                       <  1:1323769/150‑1 (MQ=40)
                                        tCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTgcag                                                                                                      >  1:1528909/1‑150 (MQ=39)
                                        tCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTgcag                                                                                                      >  1:120025/1‑150 (MQ=39)
                                           cGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGcagcag                                                                                                   >  2:95991/1‑150 (MQ=37)
                                            gCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAgg                                                                                                  >  1:1471208/1‑150 (MQ=35)
                                                  ggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgag                                                                                            <  2:1418340/150‑1 (MQ=32)
                                                   gTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagt                                                                                           >  1:561326/1‑150 (MQ=32)
                                                     cGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCg                                                                                                               <  1:120885/128‑1 (MQ=32)
                                                     cGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCg                                                                                                               >  2:120885/1‑128 (MQ=32)
                                                      gTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga                                                                                        >  1:457285/1‑150 (MQ=32)
                                                        tAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCAAACCTGCTGCAGGTtcgaatctttg                                                                                      <  1:787447/150‑8 (MQ=11)
                                                          gCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgcgc                                                                                    >  1:943377/1‑150 (MQ=32)
                                                          gCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAgcgc                                                                                    >  1:950489/1‑150 (MQ=32)
                                                           cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATTTAAAAAAGCGCACTAAAGGCGCTTTTTTGCTATCACCGATACTAAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                                                                   >  2:1132181/1‑150 (MQ=11)
                                                             cAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGGAGCAGGTTTGAGTTGAGCGCAGc                                                                                 >  2:1859818/1‑150 (MQ=17)
                                                             cAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGc                                                                                 >  2:769668/1‑150 (MQ=32)
                                                             cAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGc                                                                                 >  1:984661/1‑150 (MQ=32)
                                                              aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCg                                                                                >  2:424219/1‑150 (MQ=32)
                                                                ttGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGaa                                                                              >  1:1640180/1‑150 (MQ=32)
                                                                  ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaac                                                                            <  2:238595/150‑1 (MQ=32)
                                                                  ggTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaac                                                                            >  1:1795726/1‑150 (MQ=32)
                                                                   gTAGAGCAGTTGACTTTTAATCAATTGGGCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAaaca                                                                           <  1:871439/150‑1 (MQ=17)
                                                                          aGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGc                                                                    <  2:120025/150‑1 (MQ=32)
                                                                           gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGTTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGcc                                                                   <  1:1370918/150‑1 (MQ=17)
                                                                           gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGcc                                                                   >  1:1078223/1‑150 (MQ=32)
                                                                                cTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTcg                                                              >  2:724399/1‑150 (MQ=32)
                                                                                    tAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGCCGCATTTTTGCTATCTACCATACACAAAGATTCCAAACTGCACCAGATTTCAAAACAACACAACGAAACAACCCAACCGCTCccggc                                                          >  2:890353/1‑150 (MQ=2)
                                                                                       tCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGAc                                                       <  2:1471208/150‑1 (MQ=32)
                                                                                                cGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAg                                              >  2:802470/1‑150 (MQ=32)
                                                                                                    ggTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcg                                          >  2:1202865/1‑150 (MQ=32)
                                                                                                    ggTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTATAGCAATCTCCTATACTACAAAATTCGAAACTGCACCAGGTTTGAGATGAGCGCAGCGAAACAAAGAAAAAAAACGCGGCGACGGCCCGAAGGgcg                                          >  2:679544/1‑150 (MQ=11)
                                                                                                       tCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagc                                       <  1:802470/150‑1 (MQ=32)
                                                                                                       tCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagc                                       <  1:52289/150‑1 (MQ=32)
                                                                                                       tCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagc                                       <  1:1132181/150‑1 (MQ=32)
                                                                                                            tCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAAACTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagc                                  >  2:180491/1‑150 (MQ=17)
                                                                                                             ccTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagcg                                 <  1:724399/150‑1 (MQ=32)
                                                                                                             ccTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATATACGATACTCAAAGATTCGAACGAGGCGAAGGAATGAGTTGAGAAAAGAGAAACAAAGGAGCCGCTCGCGGAGAAGGAACGAAGGGCgagcggagcg                                 >  2:1934486/1‑150 (MQ=11)
                                                                                                             ccTGCACGACCCACCAATGTAAAAAACCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGCACCTGCAGCAGGTTTGAGTCGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagcg                                 <  2:358807/150‑1 (MQ=12)
                                                                                                                gCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGt                              <  1:1202865/150‑1 (MQ=32)
                                                                                                                gCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGATCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGt                              >  2:1547464/1‑150 (MQ=17)
                                                                                                                     aCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATcc                         >  2:214421/1‑150 (MQ=32)
                                                                                                                     aCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATcc                         <  2:264481/150‑1 (MQ=32)
                                                                                                                     aCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATcc                         <  2:110968/150‑1 (MQ=32)
                                                                                                                         accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCa                     <  2:1528909/150‑1 (MQ=32)
                                                                                                                           caATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACg                   <  1:769668/150‑1 (MQ=32)
                                                                                                                              tGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGAcc                <  2:740504/150‑1 (MQ=32)
                                                                                                                               gTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGTAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGAccc               <  1:1859818/150‑1 (MQ=17)
                                                                                                                                 aaaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac             >  1:401142/1‑150 (MQ=32)
                                                                                                                                  aaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc            >  1:1065427/1‑150 (MQ=12)
                                                                                                                                  aaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc            >  1:1055277/1‑150 (MQ=12)
                                                                                                                                  aaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc            >  1:1325241/1‑150 (MQ=32)
                                                                                                                                  aaaaaGCGACCTAAAGGCGCTTTTTTGCTATCTTCGATACTCAAAGATTCCAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc            <  1:1495484/150‑1 (MQ=12)
                                                                                                                                   aaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca           >  2:359280/1‑150 (MQ=32)
                                                                                                                                     aaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAAt         >  2:1748925/1‑150 (MQ=32)
                                                                                                                                           ccTAAAGGCGCTTTTTTGCCATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaa   >  1:1497681/1‑150 (MQ=17)
                                                                                                                                            cTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaaa  >  1:1474955/1‑150 (MQ=32)
                                                                                                                                                |                                                                                                                                                 
CCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAA  >  NC_000913/781528‑781817

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