Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,185,729 T→G 5.3% S310A (TCC→GCC)  purL → phosphoribosylformylglycinamidine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,185,7290TG 5.3% 312.4 / 9.1 115S310A (TCC→GCC) purLphosphoribosylformylglycinamidine synthase
Reads supporting (aligned to +/- strand):  ref base T (54/54);  new base G (3/3);  total (58/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-02

Only 100 of 115 total aligned reads displayed.

CGTCGGCAACGTTGCCGGCCGTTTCTTCCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTATGGCAAGCCCGAGCGTATC  >  NC_002947/1185584‑1185872
                                                                                                                                                 |                                                                                                                                                 
cgtcgGCAACGTTGCCGGCCGTTTCTTCCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCgg                                                                                                                                               <  2:616651/150‑1 (MQ=255)
 gtcgGCAACGTTGCCGGCCGTTTCTTCCCGAACCCTGAAACTCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGATCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGt                                                                                                                                              >  1:2082113/1‑150 (MQ=255)
       aaCGTTGCCGGCCGTTTCTTCCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGaa                                                                                                                                        >  1:128097/1‑150 (MQ=255)
         cGTTGCCGGCCGTTTCTTCCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGAGCCGGTggaggcgc                                                                                                                                      <  2:851377/150‑7 (MQ=255)
          gTTGCCGGCCGTTTCTTCCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATc                                                                                                                                     >  2:1821791/1‑150 (MQ=255)
             gccggccgTTTCTTCCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCcgc                                                                                                                                  >  1:514606/1‑150 (MQ=255)
              ccggccgTTTCTTCCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCcgcg                                                                                                                                 >  2:1033929/1‑150 (MQ=255)
                 gccgTTTCTTCCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCACCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGcgacg                                                                                                                              >  2:1611188/1‑150 (MQ=255)
                    gTTTCTTCCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAg                                                                                                                           <  2:2151032/150‑1 (MQ=255)
                       tcttcCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTg                                                                                                                        >  1:1244865/1‑150 (MQ=255)
                        cttcCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGc                                                                                                                       >  2:892357/1‑150 (MQ=255)
                         ttcCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCt                                                                                                                      >  2:312884/1‑150 (MQ=255)
                          tcCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTa                                                                                                                     <  2:2044467/150‑1 (MQ=255)
                           cccGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTAc                                                                                                                   >  1:626715/1‑150 (MQ=255)
                            ccGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTTCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTAcc                                                                                                                 >  1:1230997/1‑150 (MQ=255)
                            ccGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTAcc                                                                                                                 <  2:1127488/150‑1 (MQ=255)
                             cGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTAccg                                                                                                                >  2:504941/1‑150 (MQ=255)
                                aCCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTAccggcc                                                                                                             >  2:1227225/1‑150 (MQ=255)
                                  ccTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTAccggccgg                                                                                                           >  1:1143603/1‑150 (MQ=255)
                                  ccTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCAGCGGTGGCGAAATCCGCGACGAAGGTGCTAccggccgg                                                                                                           >  2:2371025/1‑150 (MQ=255)
                                  ccTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGAGCCGGTg                                                                                                                                             <  2:840771/118‑1 (MQ=255)
                                  ccTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGAGCCGGTg                                                                                                                                             >  1:840771/1‑118 (MQ=255)
                                   cTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCTCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCggg                                                                                                          >  2:742880/1‑150 (MQ=255)
                                    tGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCgggg                                                                                                         <  2:1850173/150‑1 (MQ=255)
                                    tGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGAAGAAGGTGCTACCGGCCgggg                                                                                                         >  2:231284/1‑150 (MQ=255)
                                     gAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGgc                                                                                                        <  2:514606/150‑1 (MQ=255)
                                     gAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAAGCCGCGACGAAGGTGCTACCGGCCGGGgc                                                                                                        <  2:128097/150‑1 (MQ=255)
                                           cGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAg                                                                                                  >  1:653627/1‑150 (MQ=255)
                                            gCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGATCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGc                                                                                                 <  2:2082113/150‑1 (MQ=255)
                                             ccAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGcc                                                                                                >  1:279474/1‑150 (MQ=255)
                                              cAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGccc                                                                                               <  1:312884/150‑1 (MQ=255)
                                                    ggCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCTAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCGAggcc                                                                                         >  1:798411/1‑150 (MQ=255)
                                                           tGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGa                                                                                  <  1:1147306/150‑1 (MQ=255)
                                                            gCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGAc                                                                                 <  2:717955/150‑1 (MQ=255)
                                                              aGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCg                                                                               >  1:2233080/1‑150 (MQ=255)
                                                                   ccGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGAc                                                                                 <  2:707860/143‑1 (MQ=255)
                                                                    cGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCCGCTCCGGTGGCGAACTCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCa                                                                         <  1:1145334/150‑1 (MQ=255)
                                                                          acaTCCTGATGAAGGTGGAAACGCACAACCCCCCGACCGCCACCGGCCCCTGCTCCGCCTCCTCGATCGGGTGCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCgtgt                                                                   <  1:892357/150‑1 (MQ=255)
                                                                          acaTCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGcga                                                                                                                                >  2:2085137/1‑91 (MQ=255)
                                                                          acaTCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGcga                                                                                                                                <  1:2085137/91‑1 (MQ=255)
                                                                          acaTCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCgtgt                                                                   <  1:1821791/150‑1 (MQ=255)
                                                                          acaGCCTGAGCAAGGTGGACACGCACCACCACCCGACCGCCCCCGCCACGGTCTCCGGCGCCTACACCGGATCCGGTGGCGAAATCCGCGACGAAGGTGCGACCGGCCGGGGCGCCAAGCCCAAGGCCGGCGTGACCGGCTTCACCgtgt                                                                   <  2:907083/150‑1 (MQ=255)
                                                                           caTCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCgggg                                                                                                         <  1:2275229/111‑1 (MQ=255)
                                                                            aTCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTcc                                                                 <  2:1216177/150‑1 (MQ=255)
                                                                             tCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCgggg                                                                                                         <  1:946073/109‑1 (MQ=255)
                                                                             tCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCgggg                                                                                                         >  2:946073/1‑109 (MQ=255)
                                                                             tCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCa                                                                         <  2:1404117/141‑1 (MQ=255)
                                                                             tCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCa                                                                >  1:100632/1‑150 (MQ=255)
                                                                             tCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGAGCCGGTggag                                                                                                                                          >  1:321518/1‑76 (MQ=255)
                                                                             tCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGAGCCGGTggag                                                                                                                                          <  2:321518/78‑3 (MQ=255)
                                                                               cTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAAc                                                              <  1:1033929/150‑1 (MQ=255)
                                                                               cTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAAc                                                              <  2:100632/150‑1 (MQ=255)
                                                                                   tgaAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAggccggc                                                                                      <  1:528550/122‑1 (MQ=255)
                                                                                   tgaAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAggccggc                                                                                      >  2:528550/1‑122 (MQ=255)
                                                                                    gaAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTAccggccg                                                                                                            >  2:2388730/1‑99 (MQ=255)
                                                                                    gaAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTAccggccg                                                                                                            <  1:2388730/99‑1 (MQ=255)
                                                                                     aaGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCaa                                                                                             >  1:988908/1‑113 (MQ=255)
                                                                                            aaaCGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCgg                                                 <  2:1009466/150‑1 (MQ=255)
                                                                                                 cacaACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGTCTCCGATGGCTAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTc                                            >  2:399708/1‑150 (MQ=255)
                                                                                                 cacaACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTc                                            <  1:1051473/150‑1 (MQ=255)
                                                                                                  acaACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTc                                                                          <  1:374767/119‑1 (MQ=255)
                                                                                                  acaACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTc                                                                          >  2:374767/1‑119 (MQ=255)
                                                                                                  acaACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCg                                           >  1:1896054/1‑150 (MQ=255)
                                                                                                       caccCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAg                                      >  1:1313558/1‑150 (MQ=255)
                                                                                                        accCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGc                                     >  1:2305231/1‑150 (MQ=255)
                                                                                                        accCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGc                                     <  1:1745095/150‑1 (MQ=255)
                                                                                                        accCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGc                                     >  1:377658/1‑150 (MQ=255)
                                                                                                            gACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTg                                 >  1:568873/1‑150 (MQ=255)
                                                                                                            gACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTg                                 <  2:653627/150‑1 (MQ=255)
                                                                                                                gCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCCGCGCATTCCGGGCTTCGAACAGCCTTGGGAg                             >  1:1263300/1‑150 (MQ=255)
                                                                                                                 ccATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGaaa                                                                                                                                       >  2:1705206/1‑45 (MQ=255)
                                                                                                                 ccATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGaaa                                                                                                                                       <  1:1705206/45‑1 (MQ=255)
                                                                                                                  cATCTCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCa                           <  1:742880/150‑1 (MQ=255)
                                                                                                                  cATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCGTTGGGAGCa                           >  1:1334199/1‑150 (MQ=255)
                                                                                                                    tCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAgg                         >  2:1443374/1‑150 (MQ=255)
                                                                                                                       ccccGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTgc                                                          >  2:2082211/1‑114 (MQ=255)
                                                                                                                       ccccGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTgc                                                          <  1:2082211/114‑1 (MQ=255)
                                                                                                                       ccccGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCt                      >  1:208443/1‑150 (MQ=255)
                                                                                                                       ccccGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCt                      >  2:1147436/1‑150 (MQ=255)
                                                                                                                         ccGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTAt                    >  1:1791733/1‑150 (MQ=255)
                                                                                                                         ccGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTAt                    <  2:2349200/150‑1 (MQ=255)
                                                                                                                         ccGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCCGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTAt                    <  2:1263300/150‑1 (MQ=255)
                                                                                                                           gTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCTACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTATgg                  <  1:343141/150‑1 (MQ=255)
                                                                                                                           gTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTATgg                  <  2:1244865/150‑1 (MQ=255)
                                                                                                                            ttCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCAGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCAcc                                                                       <  2:740259/96‑1 (MQ=255)
                                                                                                                            ttCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCAGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCAcc                                                                       >  1:740259/1‑96 (MQ=255)
                                                                                                                               tcCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTATGGCAAg              >  2:1825982/1‑150 (MQ=255)
                                                                                                                                  ggCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTCCAATCACGGCAGTGGGCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTATGGCAAGccc           <  2:63887/150‑1 (MQ=255)
                                                                                                                                     gcCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCt                                                                            <  2:1148910/82‑1 (MQ=255)
                                                                                                                                     gcCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCt                                                                            >  1:1148910/1‑82 (MQ=255)
                                                                                                                                      ccTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTATGGCAAGCCCGAGc       <  2:1896054/150‑1 (MQ=255)
                                                                                                                                      ccTCCACCGGCTCCGGGGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTATGGCAAGCCCGAGc       <  2:568873/150‑1 (MQ=255)
                                                                                                                                       cTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTAt                    <  1:1303683/136‑1 (MQ=255)
                                                                                                                                       cTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTATGGc                 <  1:732629/139‑1 (MQ=255)
                                                                                                                                       cTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTATGGc                 >  2:732629/1‑139 (MQ=255)
                                                                                                                                         ccACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTATGGCAAGCCCGAGCGTa    <  2:208443/150‑1 (MQ=255)
                                                                                                                                         ccACCGGAGCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTATGGCAAGCCCGAGCGTa    >  1:612624/1‑150 (MQ=255)
                                                                                                                                          cACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTATGGCAAGCCCGAGCGTAt   >  2:1213692/1‑150 (MQ=255)
                                                                                                                                           aCCGGGTCGGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTATGGCAAGCCCGAGCGTAtc  <  1:1611188/150‑1 (MQ=255)
                                                                                                                                           aCCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTATGGCAAGCCCGAGCGTAtc  <  2:983858/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CGTCGGCAACGTTGCCGGCCGTTTCTTCCCGAACCCTGAAACCCGCCAGTACGGCGCGGTGCAGGAGCCGGTGCACATCCTGATGAAGGTGGAAACGCACAACCACCCGACCGCCATCGCCCCGTTCTCCGGCGCCTCCACCGGCTCCGGTGGCGAAATCCGCGACGAAGGTGCTACCGGCCGGGGCGCCAAGCCCAAGGCCGGCCTGACCGGCTTCACCGTGTCCAACCTGCGCATTCCGGGCTTCGAACAGCCTTGGGAGCAGGCCTATGGCAAGCCCGAGCGTATC  >  NC_002947/1185584‑1185872

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