Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,366,393 A→T 5.4% F103I (TTT→ATT) ‡ PP_1191 ← S4 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,366,3930AT 5.4% 288.8 / 4.7 93F103I (TTT→ATT) ‡PP_1191S4 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (47/41);  new base T (2/3);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

CGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGG  >  NC_002947/1366252‑1366540
                                                                                                                                              |                                                                                                                                                   
cgccgcGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGccc                                                                                                                                               <  6:236244/148‑1 (MQ=255)
 gccgcGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCATTTTGCCCgg                                                                                                                                             >  4:261854/1‑148 (MQ=255)
 gccgcGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGt                                                                                                                                             >  8:314559/1‑149 (MQ=255)
   cgcgCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCACGGGCATTTTGCCCggc                                                                                                                                            <  6:12013/149‑3 (MQ=255)
   cgcgCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCACGGGCATTTTGCCCgg                                                                                                                                             >  5:12013/1‑147 (MQ=255)
   cgcgCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCATTTTGCCCggca                                                                                                                                           <  3:261854/149‑4 (MQ=255)
    gcgcACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCa                                                                                                                                          <  2:114825/149‑1 (MQ=255)
     cgcACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCACGGGCATTTTGCCCggcag                                                                                                                                          <  6:152594/149‑5 (MQ=255)
     cgcACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAg                                                                                                                                         >  2:327729/1‑149 (MQ=255)
      gcACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAgg                                                                                                                                        >  7:344505/1‑149 (MQ=255)
        aCTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGTCCGTGCAGGgc                                                                                                                                      >  3:175878/1‑149 (MQ=255)
            cAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCaaa                                                                                                                                                     >  8:46652/1‑130 (MQ=255)
             aaTGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCa                                                                                                                                 >  8:273375/1‑149 (MQ=255)
                  tgaACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCt                                                                                                                            >  3:41834/1‑149 (MQ=255)
                   gaACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCtt                                                                                                                           >  1:111438/1‑149 (MQ=255)
                    aaCTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTc                                                                                                                          <  8:344505/149‑1 (MQ=255)
                      cTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCa                                                                                                                                          >  4:79777/1‑131 (MQ=255)
                      cTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCa                                                                                                                                          <  3:79777/131‑1 (MQ=255)
                      cTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTcgc                                                                                                                        >  1:279377/1‑149 (MQ=255)
                         cTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGt                                                                                                                     <  2:154991/149‑1 (MQ=255)
                            cTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTgg                                                                                                                   <  1:218531/148‑1 (MQ=255)
                             tCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGc                                                                                                                 >  6:239514/1‑149 (MQ=255)
                             tCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGc                                                                                                                 >  3:270959/1‑149 (MQ=255)
                              ccAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGcc                                                                                                                >  6:321739/1‑149 (MQ=255)
                               cAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCCCGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCt                                                                                                               >  8:225136/1‑149 (MQ=255)
                               cAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTgg                                                                                                                   >  8:81307/1‑145 (MQ=255)
                               cAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTgg                                                                                                                   <  7:81307/145‑1 (MQ=255)
                               cAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGc                                                                                                                 <  1:253613/147‑1 (MQ=255)
                               cAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCt                                                                                                               >  8:219456/1‑149 (MQ=255)
                               cAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCt                                                                                                               <  7:314559/149‑1 (MQ=255)
                                aGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGTAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTc                                                                                                              >  1:63349/1‑149 (MQ=255)
                                            tCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAgg                                                                                                       <  5:80139/144‑1 (MQ=255)
                                            tCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAgg                                                                                                       >  6:80139/1‑144 (MQ=255)
                                            tCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGcc                                                                                                  <  1:327729/149‑1 (MQ=255)
                                                ccTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATg                                                                                              <  2:111438/149‑1 (MQ=255)
                                                 cTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGTCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGt                                                                                             <  4:175878/149‑1 (MQ=255)
                                                 cTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGt                                                                                             >  1:351415/1‑149 (MQ=255)
                                                    gATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTc                                                                                           <  1:158249/148‑1 (MQ=255)
                                                    gATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTc                                                                                           <  4:41834/148‑1 (MQ=255)
                                                            gggATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGgcg                                                                                  >  5:75784/1‑149 (MQ=255)
                                                            gggATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGgcg                                                                                  >  4:344272/1‑149 (MQ=255)
                                                               aTACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCTTTTCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCgg                                                                               <  1:307444/149‑1 (MQ=255)
                                                               aTACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGt                                                                                             >  1:48600/1‑135 (MQ=255)
                                                               aTACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGt                                                                                             <  2:48600/135‑1 (MQ=255)
                                                                   ccGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCGCATGTTCATGCTGGCGCGGGCAt                                                                           >  4:239687/1‑149 (MQ=255)
                                                                   ccGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCAt                                                                           >  7:24008/1‑149 (MQ=255)
                                                                     gCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATcg                                                                         >  8:304894/1‑148 (MQ=255)
                                                                     gCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCt                                                                         >  2:149255/1‑149 (MQ=255)
                                                                       cAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCtt                                                                        >  7:252195/1‑148 (MQ=255)
                                                                       cAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCtt                                                                        >  3:107351/1‑148 (MQ=255)
                                                                       cAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGAAGGGCGCGCACGCCATCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCtt                                                                        >  8:231502/1‑148 (MQ=255)
                                                                        aGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCa                                                                      <  7:219456/149‑1 (MQ=255)
                                                                        aGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCa                                                                      <  7:225136/149‑1 (MQ=255)
                                                                        aGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCa                                                                      <  2:279377/149‑1 (MQ=255)
                                                                          tCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACt                                                                    >  4:7449/1‑149 (MQ=255)
                                                                                 tAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATg                                                             >  4:64571/1‑149 (MQ=255)
                                                                                    acacCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCg                                                          >  7:345182/1‑149 (MQ=255)
                                                                                     cacCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCgg                                                         <  6:75784/149‑1 (MQ=255)
                                                                                      acCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCggc                                                        >  1:49166/1‑149 (MQ=255)
                                                                                          gCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCggctggc                                                    >  7:3530/1‑149 (MQ=255)
                                                                                            aGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGgctggctg                                                  >  5:158878/1‑149 (MQ=255)
                                                                                                 gCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGtt                                               >  3:165756/1‑147 (MQ=255)
                                                                                                 gCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTca                                             <  2:49166/149‑1 (MQ=255)
                                                                                                 gCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTca                                             <  7:273375/149‑1 (MQ=255)
                                                                                                       gcCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTcagcagca                                       <  5:239514/149‑1 (MQ=255)
                                                                                                       gcCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTcagcag                                         >  8:188835/1‑147 (MQ=255)
                                                                                                         cACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTcagcagca                                       >  5:372414/1‑147 (MQ=255)
                                                                                                         cACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTcagcagca                                       <  6:372414/147‑1 (MQ=255)
                                                                                                         cACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAgg                                     >  5:257545/1‑149 (MQ=255)
                                                                                                             cTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCa                                 >  3:216460/1‑149 (MQ=255)
                                                                                                              tCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCAc                                >  5:153486/1‑149 (MQ=255)
                                                                                                              tCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCAc                                >  6:97429/1‑149 (MQ=255)
                                                                                                                cAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGc                                                                                       >  3:56291/1‑92 (MQ=255)
                                                                                                                cAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGc                                                                                       <  4:56291/92‑1 (MQ=255)
                                                                                                                 aGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCg                                                          >  2:324340/1‑120 (MQ=255)
                                                                                                                 aGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCg                                                          <  1:324340/120‑1 (MQ=255)
                                                                                                                 aGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCgg                             <  4:54437/149‑1 (MQ=255)
                                                                                                                        ccACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAgg                      <  8:252195/149‑1 (MQ=255)
                                                                                                                        ccACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAgg                      <  3:344272/149‑1 (MQ=255)
                                                                                                                        ccACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAg                       <  2:63349/148‑1 (MQ=255)
                                                                                                                         cACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAggg                     >  5:164707/1‑149 (MQ=255)
                                                                                                                          aCGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGt                    <  2:351415/149‑1 (MQ=255)
                                                                                                                          aCGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGt                    <  6:153486/149‑1 (MQ=255)
                                                                                                                           cGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCAATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGtt                   >  2:359216/1‑149 (MQ=255)
                                                                                                                            gCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGttt                  <  7:188835/149‑1 (MQ=255)
                                                                                                                            gCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGttt                  >  2:135684/1‑149 (MQ=255)
                                                                                                                             cAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTc                 <  4:216460/149‑1 (MQ=255)
                                                                                                                             cAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTc                 >  3:284954/1‑149 (MQ=255)
                                                                                                                             cAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTc                 <  3:64571/149‑1 (MQ=255)
                                                                                                                              aGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTGGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCg                <  3:46912/149‑1 (MQ=255)
                                                                                                                              aGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCg                <  8:345182/149‑1 (MQ=255)
                                                                                                                                     ccGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACgg          <  1:135684/148‑1 (MQ=255)
                                                                                                                                     ccGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACgg          <  6:158878/148‑1 (MQ=255)
                                                                                                                                           cAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGTTCCAGGGTTTCGGCACGGGCGccc    <  7:231502/149‑1 (MQ=255)
                                                                                                                                            aaaaTGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTAGGCACGGGCGCCCg   >  5:214506/1‑149 (MQ=255)
                                                                                                                                             aaaTGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCgg  <  6:214506/149‑1 (MQ=255)
                                                                                                                                              aaTGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGc                                                                                       >  2:298908/1‑63 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
CGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCAGGCCTGCCGGGCAAAATGCCCGTGCAGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGG  >  NC_002947/1366252‑1366540

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: