Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 352,642 T→A 6.8% V127D (GTT→GAT)  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,6420TA 6.8% 205.1 / 2.5 69V127D (GTT→GAT) hisFimidazole glycerol phosphate synthase subunit HisF
Reads supporting (aligned to +/- strand):  ref base T (33/31);  new base A (2/3);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGGTTGGTGGTGGCGTGCGGACCGTGCAGGACATCCGCAACCTGCTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTGGGCGCTGGCGAGATCCT  >  NC_002947/352497‑352786
                                                                                                                                                 |                                                                                                                                                
aCGGTTGGTGGTGGCGTGCGGACCGTGCAGGACATCCGCAACCTGCTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGtt                                                                                                                                                 <  7:50690/147‑1 (MQ=255)
 cGGTTGGTGGTGGCGTGCGGACCGTGCAGGACATCCGCAACCTGCTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTc                                                                                                                                              <  5:10812/149‑1 (MQ=255)
    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                                                                                                            <  6:140605/109‑1 (MQ=255)
                                                                           aaCACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTg                                                                                                            >  5:140605/1‑109 (MQ=255)
                                                                           aaCACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGgcgc                                                                      >  3:217489/1‑147 (MQ=255)
                                                                           aaCACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCaa                                                                    <  8:73203/149‑1 (MQ=255)
                                                                            acacCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTc                                                                                                                       <  2:170501/97‑1 (MQ=255)
                                                                                cgccgcGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGa                                                               >  3:364769/1‑149 (MQ=255)
                                                                                cgccgcGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGa                                                               >  5:328926/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                                                   <  3:145423/148‑1 (MQ=255)
                                                                                             aaCCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATg                                                   <  8:168092/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                                              >  7:137602/1‑137 (MQ=255)
                                                                                                             gTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTc                                              <  8:137602/137‑1 (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)
                                                                                                                                                 |                                                                                                                                                
ACGGTTGGTGGTGGCGTGCGGACCGTGCAGGACATCCGCAACCTGCTCAACGCCGGTGCCGACAAGGTCTCGATCAACACCGCCGCGGTGTTCAACCCAGAGTTTGTCGGTGAGGCGGCGGACCGTTTCGGTTCGCAGTGCATCGTTGTCGCCATCGATGCCAAGAAAGTGTCCGGGCCAGGTGAAGCGCCGCGCTGGGAGATTTTCACTCACGGCGGGCGCAAGCCGACCGGGCTGGATGCAGTCGAGTGGGCGAAGAAGATGGAAGGCCTGGGCGCTGGCGAGATCCT  >  NC_002947/352497‑352786

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: ‑