Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,500,464 T→A 5.2% F789I (TTC→ATC)  PP_2193 → ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,500,4640TA 5.2% 243.0 / 2.3 77F789I (TTC→ATC) PP_2193ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (38/35);  new base A (2/2);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

AAGGACGATGCCGATCACGAGATTGACGGCTACACCACCTTCGACCTGCTGGGCAGCCAGGACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCG  >  NC_002947/2500317‑2500606
                                                                                                                                                   |                                                                                                                                              
aaGGACGATGCCGATCACGAGATTGACGGCTACACCACCTTCGACCTGCTGGGCAGCCAGGACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACgttgt                                                                                                                                                >  6:16479/1‑148 (MQ=255)
 aGGACGATGCCGATCACGAGATTGACGGCTACACCACCTTCGACCTGCTGGGCAGCCAGGACACCGGTTTCGGCACCTTCAGTGCCGGTACCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTc                                                                                                                                              <  2:349876/149‑1 (MQ=255)
           cGATCACGAGATTGACGGCTACGCCACCTTCGACCTGCTGGGCAGCCAGGACGCCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTCCAGCACGGTATGGGGCCAGCGTGCGAGGTTGTGCTATTCGCCgc                                                                                                                                    >  6:317524/1‑148 (MQ=255)
                aCGAGATTGACGGCTACACCACCTTCGACCTGCTGGGCAGCCAGGACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACTTTGTTCTATTCGCCGACCTAc                                                                                                                               <  5:317524/149‑1 (MQ=255)
                 cGAGATTGACGGCTACACCACCTTCGACCTGCTGGGCAGCCAGGACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACg                                                                                                                              >  5:143112/1‑149 (MQ=255)
                         aCGGCTACACCACCTTCGACCTGCTGGGCAGCCAGGACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGc                                                                                                                       <  3:327532/148‑1 (MQ=255)
                         aCGGCTACACCACCTTCGACCTGCTGGGCAGCCAGGACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGc                                                                                                                       <  6:143112/148‑1 (MQ=255)
                           ggCTACACCACCTTCGACCTGCTGGGCAGCCAGTACACCGGTTTCGGCACCTTCAGTGCCCGTATACAGAACCTGCTCGACAACCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGt                                                                                                                     >  8:106851/1‑148 (MQ=255)
                               acacCACCTTCGACCTGCTGGGCAGCCAGGACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTg                                                                                                                >  2:95305/1‑149 (MQ=255)
                                 accaccTTCGACCTGCTGGGCAGCCAGGACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGATCTATTCGCCGACCTACGGGCCGGCGTATCTGt                                                                                                               >  7:263654/1‑148 (MQ=255)
                                 accaccTTCGACCTGCTGGGCAGCCAGGACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGATCTATTCGCCGACCTACGGGCCGGCGTATCTGt                                                                                                               >  7:263657/1‑148 (MQ=255)
                                 accaccTTCGACCTGCTGGGCAGCCAGGACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGATCTATTCGCCGACCTACGGGCCGGCGTATCTGt                                                                                                               <  8:263654/148‑1 (MQ=255)
                                 accaccTTCGACCTGCTGGGCAGCCAGGACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGATCTATTCGCCGACCTACGGGCCGGCGTATCTGt                                                                                                               <  8:263657/148‑1 (MQ=255)
                                   caccTTCGACCTGCTGGGCAGCCAGGACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTAt                                                                                                             >  8:143412/1‑148 (MQ=255)
                                      cTTCGACCTGCTGGGCAGCCAGGACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACt                                                                                                         <  7:26124/149‑1 (MQ=255)
                                      cTTCGACCTGCTGGGCAGCCAGGACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACt                                                                                                         <  2:190542/149‑1 (MQ=255)
                                       ttCGACCTGCTGGGCAGCCAGGACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTa                                                                                                        >  7:354457/1‑149 (MQ=255)
                                        tCGACCTGCTGGGCAGCCAGGACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTAc                                                                                                       >  7:96672/1‑149 (MQ=255)
                                             ctgctgGGCAGCCAGGACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAggg                                                                                                  <  6:273129/149‑1 (MQ=255)
                                              tgctgGGCAGCCAGGACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAgggg                                                                                                 >  4:299910/1‑149 (MQ=255)
                                                    gCAGCCAGGACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTggg                                                                                           >  8:78942/1‑149 (MQ=255)
                                                    gCAGCCAGGACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTggg                                                                                           <  8:96672/149‑1 (MQ=255)
                                                         cAGGACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTAt                                                                                                             >  5:207759/1‑126 (MQ=255)
                                                         cAGGACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTAt                                                                                                             <  6:207759/126‑1 (MQ=255)
                                                          aGGACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGAcc                                                                                     >  1:267083/1‑149 (MQ=255)
                                                          aGGACACCCGTTTCGGCACCTTCAGAGCCCCTATCCAGAACCTGCTGCACAAGCAGTACAGCACGGTATGGGGACACCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATATCTAAGACAAACACGGGCGAGGGCGGAc                                                                                      >  1:300245/1‑148 (MQ=255)
                                                            gACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCta                                                                                   <  5:16479/149‑1 (MQ=255)
                                                               acCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGc                                                                                                                       <  5:238815/110‑1 (MQ=255)
                                                               acCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGc                                                                                                                       >  6:238815/1‑110 (MQ=255)
                                                                 cGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTAt                                                                                                                                           >  4:281465/1‑88 (MQ=255)
                                                                 cGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTAt                                                                                                                                           <  3:281465/88‑1 (MQ=255)
                                                                   gTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACAGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCt                                                                             >  8:130989/1‑148 (MQ=255)
                                                                        ggCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTg                                                                       <  2:267083/149‑1 (MQ=255)
                                                                         gCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCTTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTgg                                                                      <  1:95305/149‑1 (MQ=255)
                                                                         gCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTg                                                                       >  2:357573/1‑148 (MQ=255)
                                                                            ccTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAg                                                                    >  4:86005/1‑148 (MQ=255)
                                                                             cTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACATGGAGCTGGAGCa                                                                  >  1:141083/1‑149 (MQ=255)
                                                                              ttCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGAc                                                                                      >  1:350576/1‑128 (MQ=255)
                                                                              ttCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGAc                                                                                      <  2:350576/128‑1 (MQ=255)
                                                                               tCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATg                                                                >  2:68323/1‑149 (MQ=255)
                                                                                 aGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGc                                                              <  3:86005/149‑1 (MQ=255)
                                                                                  gTGCCGGGATCCAGAACCTGCTGGACAAGCAGTACAGGACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATAc                                                                                <  2:373956/130‑1 (MQ=255)
                                                                                      cGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACt                                                         >  3:76832/1‑149 (MQ=255)
                                                                                        gTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGCCTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGa                                                       >  3:239664/1‑149 (MQ=255)
                                                                                           tCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGAtt                                                     <  6:300213/148‑1 (MQ=255)
                                                                                            cgaGAACCTGCTGGACTAGCAGTAGCGCAGGGTATGGGGTCAGCGTTGGACGGTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATtg                                                    <  4:239664/146‑1 (MQ=255)
                                                                                            ccAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTggg                                                                                           <  1:142801/109‑1 (MQ=255)
                                                                                            ccAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTggg                                                                                           >  2:142801/1‑109 (MQ=255)
                                                                                             cAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATtgtg                                                  <  7:78942/149‑1 (MQ=255)
                                                                                             cAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATtgtg                                                  >  1:170710/1‑149 (MQ=255)
                                                                                             cAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATtgtg                                                  <  6:363918/149‑1 (MQ=255)
                                                                                                  ccTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGc                                             >  2:258861/1‑149 (MQ=255)
                                                                                                      ctgGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTg                                          >  5:64443/1‑148 (MQ=255)
                                                                                                      ctgGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCCCCGCCCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTg                                          >  2:255493/1‑148 (MQ=255)
                                                                                                         gACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTgg                                                                      <  8:357356/117‑1 (MQ=255)
                                                                                                         gACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTgg                                                                      >  7:357356/1‑117 (MQ=255)
                                                                                                            aaGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGc                                   >  3:249850/1‑149 (MQ=255)
                                                                                                              gCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCt                                 >  7:31171/1‑149 (MQ=255)
                                                                                                                aGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGt                               >  1:368514/1‑149 (MQ=255)
                                                                                                                 gTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGtt                              <  1:68323/149‑1 (MQ=255)
                                                                                                                 gTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGtt                              <  8:354457/149‑1 (MQ=255)
                                                                                                                      gCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAgg                          >  2:339692/1‑148 (MQ=255)
                                                                                                                      gCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGc                         <  1:339692/149‑1 (MQ=255)
                                                                                                                       cACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGc                         >  3:176292/1‑148 (MQ=255)
                                                                                                                       cACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTACACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGc                         >  2:344209/1‑148 (MQ=255)
                                                                                                                            tATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAAcgc                   <  7:100073/149‑1 (MQ=255)
                                                                                                                                gggACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAAcgccgc                <  1:258861/148‑1 (MQ=255)
                                                                                                                                gggACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAAcgccgc                >  7:229759/1‑148 (MQ=255)
                                                                                                                                gggACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAAcgccgc                <  6:64443/148‑1 (MQ=255)
                                                                                                                                gggACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAAcgccgc                >  5:15775/1‑148 (MQ=255)
                                                                                                                                 ggACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAAcgccgc                >  4:12236/1‑147 (MQ=255)
                                                                                                                                 ggACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCa              <  3:299910/149‑1 (MQ=255)
                                                                                                                                 ggACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCa              <  4:249850/149‑1 (MQ=255)
                                                                                                                                       cGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTAGTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAAt        >  3:235296/1‑149 (MQ=255)
                                                                                                                                           cGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCt     <  1:357573/148‑1 (MQ=255)
                                                                                                                                            gACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTACACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTgc   <  1:344209/149‑1 (MQ=255)
                                                                                                                                            gACGTTGTTCTATTCGCCGACCTACAGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTg    <  7:130989/148‑1 (MQ=255)
                                                                                                                                             aCGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTgcg  <  4:76832/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
AAGGACGATGCCGATCACGAGATTGACGGCTACACCACCTTCGACCTGCTGGGCAGCCAGGACACCGGTTTCGGCACCTTCAGTGCCGGTATCCAGAACCTGCTGGACAAGCAGTACAGCACGGTATGGGGACAGCGTGCGACGTTGTTCTATTCGCCGACCTACGGGCCGGCGTATCTGTATGACTACCAGGGGCGTGGGCGGACCTATACGCTGAGCTGGAGCATGGCTTACTGATTGTGGTGGCCTGAACGGGCCTGTTCAGGCAAACGCCGCCAGCAAATCCTGCG  >  NC_002947/2500317‑2500606

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