Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 352,648 C→A 7.3% A129D (GCC→GAC)  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,6480CA 7.3% 215.4 / 5.0 68A129D (GCC→GAC) hisFimidazole glycerol phosphate synthase subunit HisF
Reads supporting (aligned to +/- strand):  ref base C (33/30);  new base A (2/3);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

TTGGTGGTGGCGTGCGGACCGTGCAGGACATCCGCAACCTGCTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTGGGCGCTGGCGAGATCCTGCTGACCAGC  >  NC_002947/352501‑352796
                                                                                                                                                   |                                                                                                                                                    
ttGGTGGTGGCGTGCGGACCGTGCAGGACATCCGCAACCTGCTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGAAGTGCATCGTTGTCGc                                                                                                                                                      >  3:375997/1‑148 (MQ=255)
 tggtggtggCGTGCGGACCGTGCAGGACATCCGCAACCTGCTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCa                                                                                                                                                    <  5:362891/149‑1 (MQ=255)
   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                                                                                                                         >  3:92129/1‑138 (MQ=255)
                                       tGCTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAg                                                                                                                         <  4:92129/138‑1 (MQ=255)
                                         cTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCgg                                                                                                                              <  5:62754/131‑1 (MQ=255)
                                         cTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCgg                                                                                                                              >  6:62754/1‑131 (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:39226/1‑147 (MQ=255)
                                                                                          aCCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATg                                                             >  4:44098/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:44098/149‑1 (MQ=255)
                                                                                                                         tttCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTgg                            <  3:39226/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                  >  6:132123/1‑145 (MQ=255)
                                                                                                                                       gCATCGATGGCGACAACGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTGGGCGCTGGCga                  <  5:132123/145‑1 (MQ=255)
                                                                                                                                         aTCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTGGGCGCTGGCGAGATCct            >  6:316073/1‑149 (MQ=255)
                                                                                                                                                  gCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGAGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTGGGCGCTGGCGAGATCCTGCTGACCAg   >  1:181579/1‑149 (MQ=255)
                                                                                                                                                   ccATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTGGGCGCTGGCGAGATCCTGCTGACCAg   >  6:163600/1‑148 (MQ=255)
                                                                                                                                                   ccATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTGGGCGCTGGCGAGATCCTGCTGACCAg   >  4:30131/1‑148 (MQ=255)
                                                                                                                                                   ccATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTGGGCGCTGGCGAGATCCTGCTGACCAg   >  1:162365/1‑148 (MQ=255)
                                                                                                                                                   ccATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTGGGCGCTGGCGAGATCCTGCTGACCAGc  <  3:160189/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTGGTGGTGGCGTGCGGACCGTGCAGGACATCCGCAACCTGCTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTGGGCGCTGGCGAGATCCTGCTGACCAGC  >  NC_002947/352501‑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: