Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 352,651 T→A 6.6% I130N (ATC→AAC)  hisF → imidazole glycerol phosphate synthase subunit HisF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947352,6510TA 6.6% 216.1 / 2.8 71I130N (ATC→AAC) hisFimidazole glycerol phosphate synthase subunit HisF
Reads supporting (aligned to +/- strand):  ref base T (36/30);  new base A (2/3);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGTGGCGTGCGGACCGTGCAGGACATCCGCAACCTGCTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTGGGCGCTGGCGAGATCCTGCTGACCAGC  >  NC_002947/352504‑352796
                                                                                                                                                   |                                                                                                                                                 
gtggtggCGTGCGGACCGTGCAGGACATCCGCAACCTGCTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATc                                                                                                                                                  <  8:249068/149‑1 (MQ=255)
 tggtggCGTGCGGACCGTGCAGGACATCCGCAACCTGCTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATc                                                                                                                                                  >  5:146526/1‑148 (MQ=255)
         tGCGGACCGTGCAGGACATCCGCAACCTGCTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGATGGCGACAACGATGcactg                                                                                                                                         <  8:122790/149‑5 (MQ=255)
         tGCGGACCCTGCAGGACATCCGCAACCTGCTCAAAGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGATGGCGACAACGATGcactg                                                                                                                                         <  4:5308/149‑5 (MQ=255)
            ggACCGTGCAGGACATCCGCAACCTGCTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGaaa                                                                                                                                      <  2:15963/149‑1 (MQ=255)
              aCCGTGCAGGACATCCGCAACCTGCTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAgt                                                                                                                                    >  2:211021/1‑149 (MQ=255)
                cGTGCAGGACATCCGCAACCTGCTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAgtgt                                                                                                                                  >  7:73203/1‑149 (MQ=255)
                    cAGGACATCCGCAACCTGCTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCg                                                                                                                               <  3:358695/148‑1 (MQ=255)
                     aGGACATCCGCAACCTGCTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCAGCGGTGTTAAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTAGGTTCGCAGTGCATCGTTGTCGCCATAGATGACAAGAAAGTGTCCggg                                                                                                                             >  6:27287/1‑149 (MQ=255)
                                aaCCTGCTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAgc                                                                                                                  <  1:211021/149‑1 (MQ=255)
                                  ccTGCTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAgcgc                                                                                                                >  7:168092/1‑149 (MQ=255)
                                    tGCTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAg                                                                                                                         <  4:92129/138‑1 (MQ=255)
                                    tGCTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAg                                                                                                                         >  3:92129/1‑138 (MQ=255)
                                      cTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCgg                                                                                                                              >  6:62754/1‑131 (MQ=255)
                                      cTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCgg                                                                                                                              <  5:62754/131‑1 (MQ=255)
                                      cTCAACGCCGGTGCCGACAAGGTCTCGAGCAACACCGCCGCGGTGTTCAACCCCGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGcgccg                                                                                                              <  5:207315/147‑1 (MQ=255)
                                         aaCGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGcgcc                                                                                                               >  7:20010/1‑143 (MQ=255)
                                         aaCGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGcgcc                                                                                                               <  8:20010/143‑1 (MQ=255)
                                         aaCGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTg                                                                                                         >  5:163729/1‑149 (MQ=255)
                                           cGCAGGTTCCGACAAGTTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTggg                                                                                                       <  4:375997/149‑1 (MQ=255)
                                            gCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTgg                                                                                                        >  7:143565/1‑147 (MQ=255)
                                              cGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTgg                                                                                                        >  4:145423/1‑145 (MQ=255)
                                                     gACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTtcac                                                                                             >  2:308875/1‑149 (MQ=255)
                                                                    aaCACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTg                                                                                                                      >  5:140605/1‑109 (MQ=255)
                                                                    aaCACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTg                                                                                                                      <  6:140605/109‑1 (MQ=255)
                                                                    aaCACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGgcgc                                                                                >  3:217489/1‑147 (MQ=255)
                                                                    aaCACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCaa                                                                              <  8:73203/149‑1 (MQ=255)
                                                                     acacCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTc                                                                                                                                 <  2:170501/97‑1 (MQ=255)
                                                                         cgccgcGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGa                                                                         >  5:328926/1‑149 (MQ=255)
                                                                         cgccgcGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGa                                                                         >  3:364769/1‑149 (MQ=255)
                                                                           ccgcGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGAcc                                                                       >  6:356538/1‑149 (MQ=255)
                                                                                gtgtTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCt                                                                  <  6:146526/149‑1 (MQ=255)
                                                                                   ttCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCAcggcgg                                                                                    >  1:114358/1‑128 (MQ=255)
                                                                                   ttCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCAcggcgg                                                                                    <  2:114358/128‑1 (MQ=255)
                                                                                   ttCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTgg                                                                <  1:308875/148‑1 (MQ=255)
                                                                                   ttCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTgg                                                                <  6:350828/148‑1 (MQ=255)
                                                                                      aaCCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATg                                                             <  8:168092/148‑1 (MQ=255)
                                                                                      aaCCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATg                                                             <  3:145423/148‑1 (MQ=255)
                                                                                       aCCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATg                                                             >  4:44098/1‑147 (MQ=255)
                                                                                       aCCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATg                                                             >  4:39226/1‑147 (MQ=255)
                                                                                        cccATAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGATCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTCTCCGGGCCAGGTGAAGCGCCGCGCTCGGAGATTTTCACTGACTGCGGGCGCAAGCCGACCGGGCTGGATGCAg                                                          >  6:317718/1‑149 (MQ=255)
                                                                                         ccAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGt                                                         >  8:354178/1‑149 (MQ=255)
                                                                                          cAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTAACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTc                                                        >  2:104515/1‑149 (MQ=255)
                                                                                            gagTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGa                                                      <  6:147968/149‑1 (MQ=255)
                                                                                                ttGTCGTTTAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGTATGCAGTCGAGTg                                                   <  5:27287/148‑1 (MQ=255)
                                                                                                ttGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTg                                                   <  1:104515/148‑1 (MQ=255)
                                                                                                     ggTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCga                                              <  7:354178/148‑1 (MQ=255)
                                                                                                      gTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTc                                                        <  8:137602/137‑1 (MQ=255)
                                                                                                      gTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTc                                                        >  7:137602/1‑137 (MQ=255)
                                                                                                      gTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCg                                               >  7:122087/1‑146 (MQ=255)
                                                                                                        gAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCgaag                                            >  8:341712/1‑147 (MQ=255)
                                                                                                         aGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGaagaag                                         <  6:163729/149‑1 (MQ=255)
                                                                                                           gcggcgGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGAt                                       >  2:364502/1‑149 (MQ=255)
                                                                                                              gcggACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATgg                                     >  4:160189/1‑148 (MQ=255)
                                                                                                              gcggACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGa                                    <  8:143565/149‑1 (MQ=255)
                                                                                                               cggACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGaa                                   <  4:364769/149‑1 (MQ=255)
                                                                                                                  aCCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGc                                >  5:189018/1‑149 (MQ=255)
                                                                                                                    cGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCt                              >  7:94722/1‑149 (MQ=255)
                                                                                                                      tttCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTgg                            <  3:39226/149‑1 (MQ=255)
                                                                                                                      tttCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTgg                            <  3:44098/149‑1 (MQ=255)
                                                                                                                              cGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTGGGCGCTgg                     <  4:217489/148‑1 (MQ=255)
                                                                                                                                 ttggCATCGATGGCGACAACGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTGGGCGCTGGc                    >  7:99558/4‑146 (MQ=255)
                                                                                                                                   tGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTGGGCGCTGGc                    >  8:314607/1‑144 (MQ=255)
                                                                                                                                    gCATCGATGGCGACAACGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTGGGCGCTGGCga                  <  5:132123/145‑1 (MQ=255)
                                                                                                                                    gCATCGATGGCGACAACGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTGGGCGCTGGCga                  >  6:132123/1‑145 (MQ=255)
                                                                                                                                      aTCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTGGGCGCTGGCGAGATCct            >  6:316073/1‑149 (MQ=255)
                                                                                                                                               gCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGAGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTGGGCGCTGGCGAGATCCTGCTGACCAg   >  1:181579/1‑149 (MQ=255)
                                                                                                                                                ccATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTGGGCGCTGGCGAGATCCTGCTGACCAg   >  1:162365/1‑148 (MQ=255)
                                                                                                                                                ccATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTGGGCGCTGGCGAGATCCTGCTGACCAg   >  4:30131/1‑148 (MQ=255)
                                                                                                                                                ccATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTGGGCGCTGGCGAGATCCTGCTGACCAg   >  6:163600/1‑148 (MQ=255)
                                                                                                                                                ccATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTGGGCGCTGGCGAGATCCTGCTGACCAGc  <  3:160189/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GTGGTGGCGTGCGGACCGTGCAGGACATCCGCAACCTGCTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTGGGCGCTGGCGAGATCCTGCTGACCAGC  >  NC_002947/352504‑352796

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: