Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,282,421 A→T 12.5% intergenic (+439/+101) narI → / ← rttR nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,282,4210AT12.5% 550.8 / 13.0 201intergenic (+439/+101)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base A (48/120);  new base T (13/8);  total (70/131)
Fisher's exact test for biased strand distribution p-value = 5.03e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 203 total aligned reads displayed.

TCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGT  >  CP009273/1282273‑1282570
                                                                                                                                                    |                                                                                                                                                     
tcACCCTTCCCCGCTTGGGCATAATTTTTGATTGCGGATTCGTTTGAGAATTCCGTGGCTTTAGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGAc                                                                                                                                                      <  2:102731/150‑1 (MQ=14)
    cgttacccGGTTGGGCAGCTTAAATGATTGCGTCTGCGTTTGAGGATTCGTGTGTTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTc                                                                                                                                                  <  2:3875465/145‑1 (MQ=14)
    ccgtCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTc                                                                                                                                                  <  2:1112390/147‑1 (MQ=255)
    ccTTCCCCGCTTGGGCAGATTATTTTATTGCGGATTAGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTc                                                                                                                                                  <  2:4930542/150‑1 (MQ=255)
        ccccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCACGCTGTCGCTCGAgt                                                                                                                                              >  2:283445/1‑149 (MQ=255)
         cccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAACGCGCAATCGCAGACGCTCgcgac                                                                                                                                             >  1:4453688/1‑146 (MQ=14)
          ccGCTTGGGCAGAATATTTGATTGCGGAGTCGTTTGAGAATTCCGGGGCTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCg                                                                                                                                            <  2:5519551/150‑1 (MQ=255)
             cTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGTGGCTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAc                                                                                                                                         <  2:4060223/150‑1 (MQ=255)
              ttGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCACAGCTGAGAGCTTCGCCCATCTGGTCGTCGCATGAGGCAACGCTCTCTCGCTGACCCACGAGCCGAAca                                                                                                                                        >  1:1311891/1‑149 (MQ=14)
                gggCAGAATATTTGCTTGCGGATTCGTTTGAGACTTCCGGGGCTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCtt                                                                                                                                      <  2:3539873/150‑1 (MQ=255)
                gggCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCtt                                                                                                                                      <  1:5119045/150‑1 (MQ=255)
                gggCAGAATATTTGATCGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCGGGGGCCACGCACTCTCCCTGTCGCGCGCGGCGAccccg                                                                                                                                      >  1:4331129/1‑148 (MQ=14)
                 ggCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTGCTGCGGGAACGCAATCACGCTGTCGCTAGAGCCGaaacgaa                                                                                                                                     >  1:2860203/1‑145 (MQ=14)
                  gCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAg                                                                                                                                    <  1:938958/150‑1 (MQ=255)
                    aGATTATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTc                                                                                                                                  <  1:4626357/150‑1 (MQ=255)
                    aGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAACCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTc                                                                                                                                  >  2:1612297/1‑150 (MQ=255)
                      aaTTTTTGATTGCGGATTCGTTTGAGAATTCCGGTGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGa                                                                                                                                <  1:4941364/150‑1 (MQ=255)
                      aaTATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGa                                                                                                                                <  2:1756154/150‑1 (MQ=255)
                       atatTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGaa                                                                                                                               <  1:2410259/150‑1 (MQ=255)
                        tatTTGATTGCGGATTCGTTTGAGTATTTCGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAg                                                                                                                              <  1:2730658/150‑1 (MQ=255)
                        tatTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAg                                                                                                                              <  1:2730593/150‑1 (MQ=255)
                           ttGATTGCGGATTCGTTTGAGAATTCCGTGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCtt                                                                                                                           <  1:5204177/150‑1 (MQ=255)
                            tGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTtc                                                                                                                          <  1:3020509/150‑1 (MQ=255)
                              aTTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTGCAGCGGCAACGCACTAAAGCGGTCGCTAGAGACGAACCATAATAGAAGcataaa                                                                                                                        >  1:162907/1‑145 (MQ=14)
                                tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAt                                                                                                                      <  1:1045602/150‑1 (MQ=255)
                                     aTTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCGGCGGCAACCCTCAAACGCAGTCGCGCGGACCGAAACCTAGTGGGAGCATCTCaaacaaa                                                                                                                 >  1:4201686/1‑144 (MQ=14)
                                     aTTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCAGCGGCAAACCACTCACGCTGGCGCTCGAGCCGAACCATAGTAGAAGCATCTCaaacata                                                                                                                 >  2:4186446/1‑144 (MQ=14)
                                      ttCGTTTGAGTTTTCGGGGGCTTTTGGAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTcc                                                                                                                <  1:2463466/150‑1 (MQ=255)
                                      ttCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTGCAGAGGCAACCCTCAAAACCGGACGCACGAGCCCAACCGTAGTCGAATCTTGTCaaaataac                                                                                                                >  1:3564429/1‑143 (MQ=14)
                                        cGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTcccc                                                                                                              <  2:3457485/150‑1 (MQ=255)
                                        cGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTcccc                                                                                                              <  1:679792/150‑1 (MQ=255)
                                        cGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTcccc                                                                                                              <  2:3457352/150‑1 (MQ=255)
                                        agTTTGAGTATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTcccc                                                                                                              <  2:1995078/149‑1 (MQ=255)
                                        agTTTGAGAATTGCGGGGCTTTTGAAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTcccc                                                                                                              <  1:4925031/149‑1 (MQ=255)
                                         gTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAACCATCCCCg                                                                                                             >  2:2261365/1‑150 (MQ=255)
                                         gTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCACGCTGACGCGCGAGCCGAACCATAGTCGCAGCGACTCATACATCCCCg                                                                                                             >  2:59187/1‑150 (MQ=14)
                                          tttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCTCGCTGACCCTCGAGCCGAACCATAGTCGAAGCTTCTCAACATTaccccc                                                                                                            >  2:4409713/1‑148 (MQ=14)
                                              agtaTTCCGGTGCTTTTGATAGTGATGGTGGTGGGTGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATgg                                                                                                        <  1:4669412/147‑1 (MQ=255)
                                                 tttCCGGGGATTATGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                                     <  2:4760794/149‑1 (MQ=255)
                                                   tCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCAGCGGCAAAGCACAAAAGCTGACCCTCGAGACGAACCATAGTCGAAACATATCAAACTTACCCCCAAGGGCaca                                                                                                   >  1:2658050/1‑148 (MQ=9)
                                                    ccGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCTCTCTCGCTTACGCTCGAGACGAACCATTGTCGAATCTTCTCATACATCCCCGCATGGGCAGaa                                                                                                  >  2:3451650/1‑150 (MQ=11)
                                                      tgtttttttGTAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAata                                                                                                <  1:3748661/145‑1 (MQ=255)
                                                        ggCTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAtt                                                                                              <  2:3124396/150‑1 (MQ=38)
                                                         tcTTTTGCAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                                             <  1:1743300/149‑1 (MQ=255)
                                                         gCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAAGCAAACTCGCAGTCGCTCGAGCCCAACCATAGTCGAAGCATCACATCCATCCCCCCAAGGGCAGAAAAttt                                                                                             >  2:1371856/1‑150 (MQ=11)
                                                          tttatgtaaGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTg                                                                                            <  1:447061/143‑1 (MQ=255)
                                                          cTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCACACGCTGACGCTCGAGCCGAACCATAGTAGAAAATAAACATCCTACCCCGCATGGGCAGAATATTTg                                                                                            >  2:4541413/1‑150 (MQ=11)
                                                           ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGa                                                                                           <  2:3244431/150‑1 (MQ=38)
                                                            tttGATAGTGATGGTGGTGGGGGTTGGATTTCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAt                                                                                          <  1:3661830/150‑1 (MQ=255)
                                                             ttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGTCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAtt                                                                                         <  1:2632342/150‑1 (MQ=255)
                                                             ttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCTCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAtt                                                                                         <  2:2324716/150‑1 (MQ=255)
                                                                aaaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCg                                                                                      <  1:5251147/150‑1 (MQ=39)
                                                                   gTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCTTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAt                                                                                   <  1:2996845/150‑1 (MQ=255)
                                                                   gTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAt                                                                                   <  1:2350150/150‑1 (MQ=39)
                                                                    tGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCATTGTCGAAGCTTCTAATACTTCCCCGCATGGGCAGAATATTTTATTTCGGaat                                                                                  >  2:2958150/1‑148 (MQ=11)
                                                                    tGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCCACGCTCTCTCGCTGGCGCTCGATCCGAACCGTAGTCGTAGCTTATAAAACATCCCCGCATGGGCCGAATATTTTATTGGGGAtt                                                                                  >  2:3042471/1‑150 (MQ=11)
                                                                    tGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGGTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAtt                                                                                  <  2:1104867/150‑1 (MQ=255)
                                                                    cgATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAtt                                                                                  <  1:675245/149‑1 (MQ=37)
                                                                     gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTc                                                                                 <  1:2923119/150‑1 (MQ=39)
                                                                     gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTc                                                                                 <  1:3905460/150‑1 (MQ=39)
                                                                     gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACACACGCTGACGCACGAGCCGAACCATATTCGAAGCTTATCATACATACCCCCATGGGCCGAATAATTTATTTAGGAtta                                                                                 >  2:5112646/1‑149 (MQ=11)
                                                                      aTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCg                                                                                <  1:1995504/150‑1 (MQ=39)
                                                                         gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCtt                                                                             <  2:3722373/150‑1 (MQ=255)
                                                                         gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCCACGCACACACGCGGACGCGCGAGCCGAACCATAGTCGAAGCGTATCCTACAGCCCCCCAAGGGCAGGAGAATTGATTTGGGATTCGcgt                                                                             >  1:5287948/1‑148 (MQ=9)
                                                                               ggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGt                                                                       <  1:4997340/150‑1 (MQ=255)
                                                                                 ggatggatTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTc                                                                     <  1:193535/150‑1 (MQ=255)
                                                                                 ggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTc                                                                     <  1:780196/150‑1 (MQ=255)
                                                                                  gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCGCTAACGCTGGCGCTCGCGCCGAACCTTAGTCGAAGCTTCTCAAACCTCCCCCCATGGGCAGAATATTTTATTTCGGATTGGGTTgatcggtag                                                                    >  2:1972541/1‑143 (MQ=11)
                                                                                   aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                                   <  1:2477218/150‑1 (MQ=255)
                                                                                   aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                                   <  2:3267364/150‑1 (MQ=255)
                                                                                     ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCGTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGACTATTTGATTGCGGTTTCGCTGGGGCGTTCAggg                                                                 >  2:4473452/1‑150 (MQ=14)
                                                                                      gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGTGCCGCACCATAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCCGAATATTTGATTGCGGATTAGCTTGAGAGGTCAgggc                                                                >  2:628612/1‑149 (MQ=17)
                                                                                      gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGa                                                                <  1:586483/150‑1 (MQ=255)
                                                                                           cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtttt                                                           <  1:3948518/150‑1 (MQ=255)
                                                                                           cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAAAGCTCTCTCGGTGACGCTCGAGGCGAACCGTAGGCGCAGCTTCTCATCCTTCCCCGCCGGGGCCGAATAGTTGTTTTGGGAGTCGCGTGTGAGTTCAGGGGCgttt                                                           >  2:5351947/1‑150 (MQ=14)
                                                                                                  cTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTg                                                    <  2:1482796/150‑1 (MQ=255)
                                                                                                   tGCGCGCTTCGCCCTTCGGGTCGTTGGCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCGTTGTCGGAGCTTCTCATCCTTCCCCGCCAGGGCCGAAGATTTTTTTTCGGATTCGGTTGGGGGTTCCGGGCCGTTTGTAAGTgt                                                   >  2:5461943/1‑149 (MQ=2)
                                                                                                   tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCCGAATATTTTATTTGGGTTTCGCTTGGGAGTTAGGGGGCTTTTGTAAGTgt                                                   >  2:5462041/1‑149 (MQ=14)
                                                                                                    gcgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAT                                                  <  2:2867108‑M1/150‑1 (MQ=255)
                                                                                                     cgcgcTTCGCCCTTCGGGTCGTTGCCGGCGGCCACGCGCTCTCGCTGACGCTCGCGCCGAACCTTAGTCGCACCTTCTCATCCATCCCCCCACGGGCAGAAGATTTGATTGCGGATTCGGTTGGGAGTTCAGGGACGTTTGTAAGGGAtg                                                 >  2:4295082/1‑150 (MQ=14)
                                                                                                      gcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGG                                                <  1:3297841‑M1/150‑1 (MQ=255)
                                                                                                             gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGG                                         <  2:2598151‑M1/150‑1 (MQ=255)
                                                                                                               ccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGG                                       <  1:2499049‑M1/150‑1 (MQ=255)
                                                                                                                cTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCTAGCCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGA                                      <  2:4594291‑M1/150‑1 (MQ=255)
                                                                                                                 ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAA                                     <  2:2119103‑M1/150‑1 (MQ=255)
                                                                                                                 gtcGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAA                                     <  2:2983396‑M1/149‑1 (MQ=255)
                                                                                                                  tCGTGTCGTTGCCTGCTGCATCGCTCTCTAGCTGACTCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAG                                    <  1:5154432‑M1/150‑1 (MQ=255)
                                                                                                                  tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAG                                    <  1:802393‑M1/150‑1 (MQ=255)
                                                                                                                    gggTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCGTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGA                                  <  1:184043‑M1/150‑1 (MQ=255)
                                                                                                                     ggTCGTTGCCAGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTGAGTCGAAGCATCTCATCCATCCCCGCCTGGGCAGAATATTTGATTGCGGATTCGCGTGAGCGTTCAGGGCCTTTTGAAAGTGATGGGGGGGGGGGAAGGAg                                 >  2:220975‑M1/1‑149 (MQ=255)
                                                                                                                        cGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAc                              >  2:732540/1‑150 (MQ=255)
                                                                                                                         gTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACt                             >  2:2278113/1‑150 (MQ=255)
                                                                                                                                  ggCAACGCTCTCTCGCTGACGCTCGAGCCGAACCGTAGTCGAAGCTTCTCAACCATCCCCCCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCGgc                    >  1:3481166/1‑150 (MQ=17)
                                                                                                                                       cGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGTGGAAGGATTAATCAGCGCTGAGCGCt               >  2:4244415/1‑150 (MQ=255)
                                                                                                                                       cGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCt               >  2:3687619/1‑150 (MQ=255)
                                                                                                                                           ctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCTTGTGCAGAATATTTCGTGGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATCGTGGTGGGGGACGGATTACTCAGCGCTGCGCGCTTCGc           >  2:3943407/1‑150 (MQ=14)
                                                                                                                                           ctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCAGGGGCAGAATAGTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCGGCGCGCTTCGc           >  2:3400794/1‑150 (MQ=255)
                                                                                                                                            tctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCg                                                                                >  2:4462321/1‑80 (MQ=255)
                                                                                                                                              tcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCCTCAGAGCGGCGCCCTTCTcccg        >  2:396080/1‑149 (MQ=17)
                                                                                                                                                  tGACGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCgg    <  1:1409158/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGT  >  CP009273/1282273‑1282570

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