Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,633,785 C→A 21.9% intergenic (‑40/‑12) yjjX ← / → ytjC non‑canonical purine NTP phosphatase, ITPase/XTPase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,633,7850CA21.9% 112.1 / 25.3 75intergenic (‑40/‑12)yjjX/ytjCnon‑canonical purine NTP phosphatase, ITPase/XTPase/phosphatase
Reads supporting (aligned to +/- strand):  ref base C (34/24);  new base A (8/9);  total (42/33)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAACCCGTGCCAAAGAG  >  NC_000913/4633643‑4633931
                                                                                                                                              |                                                                                                                                                  
cGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGaaaa                                                                                                                                              >  1:222146/1‑149 (MQ=255)
 ggATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATCGCCTGAATTATAGCGGGATTGGTGGTCGCACAGACAACTTCGTCTATAATCAGCAATACTCATAAAATTAACGTTACAGCATTCAACGaaaaaa                                                                                                                                             >  2:665590/1‑149 (MQ=255)
  gATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCGGAATTTTAGCGGGAGTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGaaaaaa                                                                                                                                            >  2:153221/1‑149 (MQ=255)
   aTGCAATATGGCAGGATCAGTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACATCTTGGTGCATAATCAGCATTACTCAGAATATTAACGTTACAGCAGAATACGGaaaaaaa                                                                                                                                           <  2:354643/149‑1 (MQ=255)
     gCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGc                                                                                                                                         <  1:808439/149‑1 (MQ=255)
     gCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGc                                                                                                                                         <  1:217832/149‑1 (MQ=255)
       aaTATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCTGAATCGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGAATACTTCGTGCAGAATCAGCATTACTAAGAAAATTAACGTTACAGCATTATACGGAAAAAAAGCag                                                                                                                                       >  1:596189/1‑148 (MQ=255)
          atGGCAGGATCCTTCGCCGAAGATCTCGAGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATGTCTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGtt                                                                                                                                    >  2:560498/1‑149 (MQ=255)
            ggCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTCCTCAGAAAATTAACGTTACAGCAGTATACGCAAAAAAAGCATGTTAc                                                                                                                                  >  2:98866/1‑149 (MQ=255)
             gCAGGATCCTTCGCCGAAGAGCTCGTGAAATGCCTGCACAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTTCAGCAGTATACGGAAAAAAAGCAAGTTACa                                                                                                                                 >  2:721666/1‑149 (MQ=255)
               aGGTTCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAgg                                                                                                                               <  1:405764/149‑1 (MQ=255)
                ggATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGAGTGGGGGTCGCACAGACAACTTGGTGCATAATCAGCATAACTCACAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGt                                                                                                                              <  2:245055/149‑1 (MQ=255)
                  aTCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGtat                                                                                                                            <  2:813145/149‑1 (MQ=255)
                    ccTTAGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATAc                                                                                                                          <  1:237172/149‑1 (MQ=255)
                     cTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATAcc                                                                                                                         >  1:48185/1‑149 (MQ=255)
                      ttCGCCGAAGATCTCGTGTAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCt                                                                                                                        >  1:554363/1‑149 (MQ=255)
                      ttCGCCGAAGATCTCGTGAAATGCCTGCTGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTATCGTTACAGCAGTAGACGGAAAAAAAGCATGTTTCAGGTATACCt                                                                                                                        >  2:871401/1‑149 (MQ=255)
                       tCGCCGAAGATCTCGTTAAATTCCAGCAGAATTGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTTCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTa                                                                                                                       <  2:1003442/149‑1 (MQ=255)
                         gCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGCCAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAACAAGCATGTTACAGGTATACCTAGt                                                                                                                     >  2:742020/1‑149 (MQ=255)
                          ccGAACATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCTGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACCGCAGTATACGGAAAAAAAGCATGTTACAGGGATACCTAGTc                                                                                                                    <  1:1023148/149‑1 (MQ=255)
                           ctaaGATCTGGTGAAATCCCTGCAGAATTGCCTGAATTTTAGCGTGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTcc                                                                                                                   <  2:288912/147‑1 (MQ=255)
                           cGAAGATCTAGTGAAATGCCTGCAGAATGGCCTGAATTATAGCGGGATTGGTGGTCGCACTGACAACTTGGTGCATAATCAGCATTACTCAGAAATTTAACGTTACAGCAGTTTAATGAAAAAAAGCATGTTACAGGTATACCTAGTcc                                                                                                                   >  2:863166/1‑149 (MQ=255)
                               gATCTCGGGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCa                                                                                                               <  2:345853/149‑1 (MQ=255)
                                      tcaaaTGCCTGCATAATGGCAGGAATTTAAGCGTGATTGGTGGTCGCAAAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACCGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGaa                                                                                                        <  2:234337/147‑1 (MQ=255)
                                      tGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGaa                                                                                                        <  1:420796/149‑1 (MQ=255)
                                       gAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGaaa                                                                                                       <  1:272442/149‑1 (MQ=255)
                                        aaaTGCCTGCAGAATGGCCTGAATTTTAGCCTGATGGGTGGTCGCACAGACAACTTGGTGCATACTCAGCATTACTCAGAAAATGAACGTTACAGCAGGAGACGGAACAAAAGCATGTTACATGTATACCTAGTCCGCAACGGAGCAAc                                                                                                      >  2:662306/1‑149 (MQ=255)
                                              cTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTg                                                                                                <  1:22848/149‑1 (MQ=255)
                                               tGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGt                                                                                                           >  1:17444/1‑137 (MQ=255)
                                               tGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGt                                                                                                           <  2:17444/137‑1 (MQ=255)
                                                gCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCGGCATAACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATCTTACAGGTATACCTAGTCCGCCACGCTGAAACGCAGTggc                                                                                              >  1:805751/1‑148 (MQ=255)
                                                gCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGa                                                                                              >  2:523765/1‑149 (MQ=255)
                                                     aTGGCCTGAATTTTATCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCAATACTCAGAAAAATAACGTTACAGCAGTATACGGAAAAAAAGCCTGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGcc                                                                                         <  1:728637/149‑1 (MQ=255)
                                                          cTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATCCCTAGTCCGCCACGGTGAACCGCAGTGGACCGCcgcgcg                                                                                    >  2:645117/1‑149 (MQ=255)
                                                          cTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCg                                                                                    >  2:549375/1‑149 (MQ=255)
                                                           tGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGa                                                                                   >  1:870180/1‑149 (MQ=255)
                                                              aTTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATAAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGt                                                                                <  1:530364/149‑1 (MQ=255)
                                                                     cGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAATGAAAAAAAGCATGTTTCAGTTATACCTTGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCCgg                                                                         >  2:462161/1‑149 (MQ=255)
                                                                     cGGGAGTGGTGGGAGCACAGACAACTTGGTGCATAATCAGCATTACTCAGACAATTCACGTTACAGCAGTATACGGAAAAAACGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAgg                                                                         <  1:98866/149‑1 (MQ=255)
                                                                      tggtttGGTGGTCGCACAGAAAACTTGGTGTATTTTGAGGATTACTCAGAAAATTAACGTTAGAGCAGTATACGGTAATAAAGCATGTTACATGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTAAggg                                                                        <  1:665590/145‑1 (MQ=255)
                                                                       ggATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTACGCGACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGc                                                                       >  2:56621/1‑149 (MQ=255)
                                                                       ggATTGGTGGGCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGc                                                                       <  1:698366/149‑1 (MQ=255)
                                                                                gtCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAAGGTTACAGCAGTATACGGAAAAAACGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTAGTCAGGGCCAGTCTGAc                                                              >  2:203228/1‑149 (MQ=255)
                                                                                 tCGGACAGACAACTTGGTGCATAATCAGCATTAATCAGAAATTTAACGGTACAGCTGTATACGGAATAAATGCATGTTACAGGTATTCCTAGTCCGCCACGGGGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACa                                                             <  2:77491/149‑1 (MQ=255)
                                                                                 tCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAATCATGTTACAGGTATACCAAGTCCGCCACGGTGAAACGCAGTGTAACGCCGAGCGAAGTATTCAGGGCCAGTCTGACa                                                             <  2:718250/149‑1 (MQ=255)
                                                                                     acaGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGCAAAAAAAGCATGTTACAGGTATACCAAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCGGACAGccc                                                         >  2:693331/1‑149 (MQ=255)
                                                                                       aGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGc                                                       >  2:804149/1‑149 (MQ=255)
                                                                                           aaCTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGAc                                                   >  2:817461/1‑149 (MQ=255)
                                                                                           aaCTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATCCCTAGTCCGCCACGGGGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGACAGTCTGACAGCCCGCTGAc                                                   >  2:469242/1‑149 (MQ=255)
                                                                                                 gTGCATAATCAGCATTACTCAGAAAATTAACGGTACAGCAGGTTAAGGAAAAAAAGCATGTTACTGGATTACCTAGTCGGCCACGGTGAAGCGCAGGGGATCGCCGAGCGGCGGGGTCCGGGGCGGGGGGCCGGCCCGCTGACCGccga                                             >  2:1002945/1‑147 (MQ=255)
                                                                                                   gaaTAATCAGCATTACTCATAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAg                                           <  2:468090/147‑1 (MQ=255)
                                                                                                   gCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGGATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAg                                           >  2:9853/1‑149 (MQ=255)
                                                                                                    cATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAgg                                          >  1:458519/1‑149 (MQ=255)
                                                                                                       aaTCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGa                                       <  1:153221/149‑1 (MQ=255)
                                                                                                       aaTCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGa                                       <  1:877568/149‑1 (MQ=255)
                                                                                                         tCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGc                                     >  2:132022/1‑149 (MQ=255)
                                                                                                            gCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGAAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAAc                                  <  2:222146/149‑1 (MQ=255)
                                                                                                              aTTACTCAGAAAATTAACGTTACAGCAGTATACGGACAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGGGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAg                                <  2:422998/149‑1 (MQ=255)
                                                                                                              aTTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAg                                >  1:557677/1‑149 (MQ=255)
                                                                                                               ttACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGAATGTTACAGGTATAACTAGTCCGCCCCGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAgg                               >  1:68011/1‑149 (MQ=255)
                                                                                                                tACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGc                              >  1:111708/1‑149 (MQ=255)
                                                                                                                 aCTCAGAAAACTAACGTTACAGCAGTATACGGAAAAAAAGCATGATACAGGTATACCTAGTCCGCCACGGGGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCg                             <  1:645117/149‑1 (MQ=255)
                                                                                                                   tCAGAAAATTAACGGTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTCGTCCGCCACGGGGAAACGCAGTTGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGAt                           <  1:560498/149‑1 (MQ=255)
                                                                                                                     aGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGCCAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGc                         >  2:334074/1‑149 (MQ=255)
                                                                                                                     aGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGc                         >  2:339124/1‑149 (MQ=255)
                                                                                                                      gAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCa                        >  1:298541/1‑149 (MQ=255)
                                                                                                                       aaaaTTTACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAg                       >  1:258022/1‑149 (MQ=255)
                                                                                                                       aaaaTTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGccg                       >  1:364812/1‑147 (MQ=255)
                                                                                                                         aatacaCGATACAGCAGTATACGGAAACAAAGCATGTTACAGGTATACCTAGTCCGCCCCGGGGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGt                     <  1:721666/144‑1 (MQ=255)
                                                                                                                          aTTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGAAGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTg                    >  1:423469/1‑149 (MQ=255)
                                                                                                                           tttaagTTACTCCATTATAAGGAAAAATATCATGTAACATGTATTCCTAGTCCTCCACGGTGAAAAGCAGTGTAACGAGGAGCGACGTTATCAGGGCAAGTCTGACAGCCCGCTGACCGCCATAGTTGAGCAACAGGCGATGCAGGTgg                   <  1:983551/144‑1 (MQ=255)
                                                                                                                                   aCAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGTAGGTGGCAACCCGt           <  2:870180/149‑1 (MQ=255)
                                                                                                                                     agcagTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAACCCGTGc         >  1:178900/1‑149 (MQ=255)
                                                                                                                                         gTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGAGAAGTATTCAGGGCCAATCTGACAGCCCGCTGACCGACAAAGGTGAGCAACAGGCGATGAATGTGGCAACACGTGCCaaa     >  1:644455/1‑149 (MQ=255)
                                                                                                                                            taCGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCCGGTGGCCACCCGTGCCCAcgag  >  1:127429/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
CGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAACCCGTGCCAAAGAG  >  NC_000913/4633643‑4633931

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