Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,471,656 A→T 6.0% intergenic (‑467/‑427) PP_3961 ← / → PP_3962 transporter/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,471,6560AT 6.0% 204.7 / 2.6 66intergenic (‑467/‑427)PP_3961/PP_3962transporter/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/34);  new base T (2/2);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

GAGGACTTGCCGCATGCCGCGAAGGTTGAGGATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGA  >  NC_002947/4471511‑4471798
                                                                                                                                                  |                                                                                                                                              
gaggaCTTGCCGCATGCCGCGAAGGTTGAGGATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGaataa                                                                                                                                             >  3:333366/1‑149 (MQ=255)
    aCTTGCCGCATGCCGCGAAGGTTGAGGATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAg                                                                                                                                         <  5:139058/149‑1 (MQ=255)
      ttGCCGCATGCCGCGAAGGTTGAGGATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAggg                                                                                                                                       >  3:131535/1‑149 (MQ=255)
            cATGCCGCGAAGGTTGAGGATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCAGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCCTGATTATTCAAGGCttttt                                                                                                                                 <  4:300053/149‑3 (MQ=255)
                ccGCGAAGGTTGAGGATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCAt                                                                                                                             >  7:321820/1‑149 (MQ=255)
                  gcgAAGGTTGAGGATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGc                                                                                                                           >  3:271953/1‑149 (MQ=255)
                        gTTGAGGATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTgg                                                                                                                     <  2:182259/149‑1 (MQ=255)
                        gTTGACGATTACGTATCGCTATTGCCGTCGAACGCAAAGTCGGGGAAGCAACGCTAAACCAAGCCCTGACCTTTCGGGAGGAATGGAACCTGGATCGGGTAAGAAACGCAAAATTCCTTGAATAATCAGGGCTTTCCTCTTGCctcttg                                                                                                                     <  2:134742/149‑3 (MQ=255)
                           gAGGATTACGAATCGCTATTGCGTTGGAAGTTAAAGTCGGTGAAGCCACGTTAAACCAAGCCATGAGCTTGCGTCAGGAATTGAAGGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCAt                                                                                                                             <  5:302643/138‑1 (MQ=255)
                           gAGGATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTAGGATCGGTTAAGAAACGCAAAAAGCCTTTAATAATCAGGGCTTTCCTCATGCCTCTGGc                                                                                                                    <  1:342206/148‑1 (MQ=255)
                              gATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCAtt                                                                                                               >  3:406186/1‑149 (MQ=255)
                                   cGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGc                                                                                                           <  4:333366/148‑1 (MQ=255)
                                   cGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCACTGGCCCATTGGGc                                                                                                           <  6:19847/148‑1 (MQ=255)
                                      aTCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGgc                                                                                                       >  8:395313/1‑149 (MQ=255)
                                         gCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGgcgc                                                                                                     <  2:287492/148‑1 (MQ=255)
                                          cTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGgcgcgc                                                                                                   >  5:82463/1‑149 (MQ=255)
                                           tATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCt                                                                                                  >  3:20962/1‑149 (MQ=255)
                                             ttGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGAGCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCt                                                                                                  <  7:296533/147‑1 (MQ=255)
                                                  gTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTa                                                                                           <  6:124928/149‑1 (MQ=255)
                                                        gTGCAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGTAACGTGGATCGTTTAAGAACCGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCAttt                                                                                     <  6:142997/149‑1 (MQ=255)
                                                        gTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCAttt                                                                                     <  6:82463/149‑1 (MQ=255)
                                                        gTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCACATTGGGCCGGCGCGCTCCGCTTACCAttt                                                                                     <  2:118614/149‑1 (MQ=255)
                                                          gAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCTGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGa                                                                                   >  6:207995/1‑149 (MQ=255)
                                                               tCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACAcc                                                                               >  8:97537/1‑148 (MQ=255)
                                                                  gggAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGa                                                                           <  4:20962/149‑1 (MQ=255)
                                                                     aaGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTc                                                                        <  7:97537/149‑1 (MQ=255)
                                                                       gCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCaa                                                                      <  4:131535/149‑1 (MQ=255)
                                                                       gCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCaa                                                                      <  7:395313/149‑1 (MQ=255)
                                                                         cACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAAt                                                                     <  6:114450/148‑1 (MQ=255)
                                                                         cACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAAt                                                                     >  5:114450/1‑148 (MQ=255)
                                                                         cACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAAt                                                                     >  2:213082/1‑148 (MQ=255)
                                                                          aCGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGa                                                                   >  3:360977/1‑149 (MQ=255)
                                                                              tAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTg                                                               <  4:271953/149‑1 (MQ=255)
                                                                                  ccAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAAcc                                                           >  2:125554/1‑149 (MQ=255)
                                                                                   cAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGAGACCAAGAGTCAATGACGTGAACCg                                                          >  8:365876/1‑149 (MQ=255)
                                                                                    aaGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGt                                                         >  5:61387/1‑149 (MQ=255)
                                                                                    aaGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGt                                                         >  3:77324/1‑149 (MQ=255)
                                                                                       cccTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGc                                                      <  8:321820/149‑1 (MQ=255)
                                                                                        ccTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAGCCGTCGCa                                                     >  1:50452/1‑149 (MQ=255)
                                                                                         cTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCa                                                     >  2:165761/1‑148 (MQ=255)
                                                                                         cTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCa                                                     >  8:190384/1‑148 (MQ=255)
                                                                                         cTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCa                                                     >  6:143157/1‑148 (MQ=255)
                                                                                          tGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCata                                                   >  2:162132/1‑149 (MQ=255)
                                                                                              cTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTc                                               >  7:201486/1‑149 (MQ=255)
                                                                                               ttGCTGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCg                                              <  5:207995/149‑1 (MQ=255)
                                                                                                   ggCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCg                                              >  4:114548/1‑145 (MQ=255)
                                                                                                   ggCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGa                                          <  4:406186/149‑1 (MQ=255)
                                                                                                          aGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCAcgcg                                   <  7:75802/149‑1 (MQ=255)
                                                                                                           gggAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGc                                             >  3:264022/1‑138 (MQ=255)
                                                                                                           gggAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGc                                             <  4:264022/138‑1 (MQ=255)
                                                                                                            ggAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCATCATGACTCTGGCCAATTGGGCCGGCCCGAAACGCTTACCATTTGACACCAACATTCAATGACGTGAACCGACGCAAATGTCGCAGATCAcaagaa                                 >  2:83311/1‑144 (MQ=255)
                                                                                                                cGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATggg                             <  7:190384/149‑1 (MQ=255)
                                                                                                                         ggTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCt                     <  6:61387/149‑1 (MQ=255)
                                                                                                                         ggTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCt                     >  8:88478/1‑149 (MQ=255)
                                                                                                                          gTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCAt                                >  6:65887/1‑137 (MQ=255)
                                                                                                                          gTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCAt                                <  5:65887/137‑1 (MQ=255)
                                                                                                                           ttAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGc                   >  6:23148/1‑149 (MQ=255)
                                                                                                                           ttAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGc                   <  3:114548/149‑1 (MQ=255)
                                                                                                                              aGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGtt                <  4:77324/149‑1 (MQ=255)
                                                                                                                                aaaCGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTgg              <  1:213082/149‑1 (MQ=255)
                                                                                                                                 aaCGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGc             <  1:165761/149‑1 (MQ=255)
                                                                                                                                   cGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAGCCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGa           <  2:50452/149‑1 (MQ=255)
                                                                                                                                     cAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAAc         <  5:92358/149‑1 (MQ=255)
                                                                                                                                     aggaaaGCCCTGATTATTCAAGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAAc         >  6:171171/4‑149 (MQ=255)
                                                                                                                                      ggaaaGCCCTGATTATTCAAGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAAcc        >  5:392745/3‑149 (MQ=255)
                                                                                                                                      aaaaaGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAAcc        <  3:16348/149‑1 (MQ=255)
                                                                                                                                            ccctGATTATTCAAGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGa  <  6:392745/146‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GAGGACTTGCCGCATGCCGCGAAGGTTGAGGATTACGAATCGCTATTGCCGTGGAAGTGAAAGTCGGGGAAGCCACGGTAAACCAAGCCCTGACCTTGCGGCAGGAAGGGAACGTGGATCGGTTAAGAAACGCAAAAAGCCTTGAATAATCAGGGCTTTCCTCATGCCTCTGGCCCATTGGGCCGGCGCGCTCCGCTTACCATTTGACACCAAGATTCAATGACGTGAACCGTCGCATATGTCGCAGATCACGCGCATGGGCTGGCCCTGCGTTGGCGAACCCTGCGA  >  NC_002947/4471511‑4471798

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