Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,175,593 T→C 6.3% I284T (ATC→ACC)  PP_2787 → transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,175,5930TC 6.3% 181.3 / 2.0 63I284T (ATC→ACC) PP_2787transporter
Reads supporting (aligned to +/- strand):  ref base T (23/36);  new base C (2/2);  total (25/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-01

ACCAAACAGGACGACGCTCAGGCCCAGTTGCGTGCCTATCACCAGGTCCGCGAGATCATCGCGAAGTACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCGTGGTGCTGCCGCTGTTCAATATCCTGTTCTGCACCATCTGGACCCTTGGCCTGATGGCGCTGC  >  NC_002947/3175445‑3175742
                                                                                                                                                    |                                                                                                                                                      
aCCAAACAGGACGACGCTCAGGCCCAGTTGCGTGCCTATCACCAGGTCCGCGAGATCATCGCGAAGTACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAatcatc                                                                                                                                                       >  1:457989/1‑150 (MQ=255)
aCCAAACAGGACGACGCTCAGGCCCAGTTGCGTGCCTATCACCAGGTCCGCGAGATCATCGCGAAGTACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAatcatc                                                                                                                                                       <  2:457989/150‑1 (MQ=255)
aCCAAACAGGACGACGCTCAGGCCCAGTTGCGTGCCTATCACCAGGTCCGCGAGATCATCGCGAAGTACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAatcatc                                                                                                                                                       >  2:1748276/1‑150 (MQ=255)
                  cAGGCCCAGTTGCGTGCCTATCACCAGGTCCGCGAGATCATCGCGAAGTACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATg                                                                                                                                     <  2:1139480/150‑1 (MQ=255)
                  cAGGCCCAGTTGCGTGCCTATCACCAGGTCCGCGAGATCATCGCGAAGCACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATg                                                                                                                                     <  2:578328/150‑1 (MQ=255)
                    ggCCCAGTTGCGTGCCTATCACCAGGTCCGCGAGATCATCGCGAAGTACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCt                                                                                                                                   <  1:1128895/150‑1 (MQ=255)
                     gCCCAGTTGCGTGCCTATCACCAGGTCCGCGAGATCATCGCGAAGTACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTg                                                                                                                                  >  2:1090144/1‑150 (MQ=255)
                         aGTTGCGTGCCTATCACCAGGTCCGCGAGATCATCGCGAAGTACCAGGCCTCCCACCCGGAGTGCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGTCGATGCTGTTcc                                                                                                                              >  2:1054698/1‑150 (MQ=255)
                         aGTTGCGTGCCTATCACCAGGTCCGCGAGATCATCGCGAAGTACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTcc                                                                                                                              >  1:249648/1‑150 (MQ=255)
                         aGTTGCGTGCCTATCACCAGGTCCGCGAGATCATCGCGAAGTACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTcc                                                                                                                              <  2:13254/150‑1 (MQ=255)
                              cGTGCCTATCACCAGGTCCGCGAGATCATCGCGAAGTACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATc                                                                                                                         <  1:305505/150‑1 (MQ=255)
                                     aTCACCAGGTCCGCGAGATCATCGCGAAGTACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTcc                                                                                                                 <  1:1090144/150‑1 (MQ=255)
                                        aCCAGGTCCGCGAGATCATCGCGAAGTACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCAc                                                                                                                                                 <  2:224851/116‑1 (MQ=255)
                                        aCCAGGTCCGCGAGATCATCGCGAAGTACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCAGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCAc                                                                                                                                                 >  1:224851/1‑116 (MQ=255)
                                         ccAGGTCCGCGAGATCATCGCGAAGTACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGtt                                                                                                             <  2:752875/150‑1 (MQ=255)
                                          cAGGTCCGCGAGATCACCGCGAAGTACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTg                                                                                                            <  2:262349/150‑1 (MQ=255)
                                            ggTCCGCGAGATCATCGCGAAGTACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGt                                                                                                          <  1:1896285/150‑1 (MQ=255)
                                             gTCCGCGAGATCATCGCGAAGTACCAGGCCGGCCACCCGGAGGTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCGTTGCGGCCGAGCAGGAAGGCATCGATCACGTCCTGGCGGTGCTGTTCGCGATCGTGGTTCttgttgtg                                                                                                          >  2:395233/1‑148 (MQ=255)
                                              tCCGCGAGATCATCGCGAAGTACCAGGCCGCCCACCCGGAGTTCACCGATGAGCTGTTCATCGCCGGTATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGa                                                                                                        <  2:1032497/150‑1 (MQ=255)
                                                 gcgAGATCATCGCGAAGTACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACct                                                                                                     <  2:1335257/150‑1 (MQ=255)
                                                 gcgAGATCATCGCGAAGTACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACct                                                                                                     >  2:389183/1‑150 (MQ=255)
                                                      atcatcGCGAAGTACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGAc                                                                                                                                              <  1:1739972/105‑1 (MQ=255)
                                                      atcatcGCGAAGTACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGAc                                                                                                                                              >  2:1739972/1‑105 (MQ=255)
                                                           cgcgAAGTACCAGGCCGCCCACCCGGAGTTCACCCCGTATCTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGt                                                                                           <  1:1040521/150‑1 (MQ=255)
                                                                  tACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCACCGATCa                                                                                                                                                  >  2:1405240/1‑89 (MQ=255)
                                                                  tACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCACCGATCa                                                                                                                                                  >  2:1405495/1‑89 (MQ=255)
                                                                  tACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCACCGATCa                                                                                                                                                  <  1:1405495/89‑1 (MQ=255)
                                                                  tACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCACCGATCa                                                                                                                                                  <  1:1405240/89‑1 (MQ=255)
                                                                     cAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCttt                                                                                 >  2:114828/1‑150 (MQ=255)
                                                                     cAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCttt                                                                                 <  1:549027/150‑1 (MQ=255)
                                                                     cAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGATGTTTTTCTTCttt                                                                                 >  2:1080209/1‑150 (MQ=255)
                                                                       ggCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGCCCTCTCTGCTGGTTTTCTTCTTTct                                                                               >  2:1885394/1‑149 (MQ=255)
                                                                               aCCCGGAGTGCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCAc                                                                       >  1:1547465/1‑150 (MQ=255)
                                                                                cccGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACt                                                                      >  2:455588/1‑150 (MQ=255)
                                                                                cccGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACt                                                                      >  1:867476/1‑150 (MQ=255)
                                                                                 ccGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTcc                                                                                                                              >  1:727691/1‑94 (MQ=255)
                                                                                 ccGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTcc                                                                                                                              <  2:727691/94‑1 (MQ=255)
                                                                                         cACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCgtggt                                                             <  1:389183/150‑1 (MQ=255)
                                                                                          aCCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCgtggtg                                                            <  1:841921/150‑1 (MQ=255)
                                                                                            cGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCgtg                                                               >  2:504370/1‑145 (MQ=255)
                                                                                            cGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCgtg                                                               <  1:504370/145‑1 (MQ=255)
                                                                                               tGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGAt                                                                                                                                      <  2:625939/72‑1 (MQ=255)
                                                                                               tGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGAt                                                                                                                                      >  1:625939/1‑72 (MQ=255)
                                                                                               tGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTACTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCGTGGTGCTGcc                                                       <  1:1054698/150‑1 (MQ=255)
                                                                                                 aGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCGTGGTGCTGCCGc                                                     >  1:1705114/1‑150 (MQ=255)
                                                                                                  ggTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCGTGGTGCTGCCGCt                                                    <  1:1306887/150‑1 (MQ=255)
                                                                                                   gtgtTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCGTGGTGCTGCCGCTg                                                   <  1:1748276/150‑1 (MQ=255)
                                                                                                       tCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCGTGGTGCTGCCGCTGTTCa                                               <  2:1378984/150‑1 (MQ=255)
                                                                                                          tCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCGTGGTGCTGCCGCTGTTCAata                                            <  1:1689301/150‑1 (MQ=255)
                                                                                                             ccGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCGTGGTGCTGCCGCTGTTCAATATcc                                         <  2:249648/150‑1 (MQ=255)
                                                                                                              cGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCGTGGTGCTGCCGCTGTTCAATATCCt                                        >  2:1635921/1‑150 (MQ=255)
                                                                                                                 cATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCGTGGTGCTGCCGCTGTTCAATATCCTGtt                                     <  1:1902851/150‑1 (MQ=255)
                                                                                                                 cATGCCGATCTTCATTCGGGCCCAACAGGAAATCATCGATCATGACCTGGCGATGCTGTGCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCGGGTTTTCTTCTTTCGCAAACCACTTGGCGTGGTGCTGCCGCTGTTCAATATCCTGtt                                     <  1:292602/150‑1 (MQ=255)
                                                                                                                   tGCCGATCTTCATTGCGGCCCAGGAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTCTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCCCTTGGCGTGGTACTGCCGCTGTTCAATAGCCTGTTCt                                   >  2:946177/1‑150 (MQ=255)
                                                                                                                   tGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCGTGGTGCTGCCGCTGTTCAATATCCTGTTCt                                   <  1:822613/150‑1 (MQ=255)
                                                                                                                     ccGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTttcttc                                                                                    <  2:1017320/99‑1 (MQ=255)
                                                                                                                     ccGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTttcttc                                                                                    >  1:1017320/1‑99 (MQ=255)
                                                                                                                      cGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTTGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTATCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCGTGGTGCTGCCGCTGTTCAATATCCTGTTCTGCa                                <  2:1095719/150‑1 (MQ=255)
                                                                                                                       gATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCGTGGTGCTGCCGCTGTTCAATATCCTGt                                      >  1:1054105/1‑143 (MQ=255)
                                                                                                                       gATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCGTGGTGCTGCCGCTGTTCAATATCCTGt                                      <  2:1054105/143‑1 (MQ=255)
                                                                                                                              aTTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCGTGGTGCTGCCGCTGTTCAATCTCCTGTTCTGCACCATCTgg                        >  1:143404/1‑150 (MQ=255)
                                                                                                                                    gCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCGTGGTGCTGCCGCTGTTCAATATCCTGTTCTGCACCATCTGGACCCtt                  <  2:547482/150‑1 (MQ=255)
                                                                                                                                     cccAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCGTGGTGCTGCCGCTGTTCAATATCCTGTTCTGCACCATCTGGACCCTTg                 <  1:220047/150‑1 (MQ=255)
                                                                                                                                     cccAACAGGAAATCATCGATCACTCCGTGGCGATGCTCTTCCCGATCGTGTTTCTGTTCGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCTCTTGGCGTGCCGCTCCCGCTGCTCAATATCCTGTTCTGCACCATCTGGACCCTTg                 <  2:53456/150‑1 (MQ=255)
                                                                                                                                           agGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCGTGGTGCTGCCGCTGTTCAATATCCTGTTCTGCACCATCTGGACCCTTGGCCTGa           >  2:619593/1‑150 (MQ=255)
                                                                                                                                            ggAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCGTGGTCCGGCCGCTGTTCAATATCCTGTTCTGCACCATCTGGACCCTTGGCCTGAt          <  1:987007/150‑1 (MQ=255)
                                                                                                                                                   atcGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCGTGGTGCTGCCGCTGTTCAATATCCTGTTCTGCACCATCTGGACCCTTGGCCTGATGGCgctg   >  2:420111/1‑150 (MQ=255)
                                                                                                                                                   atcGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTATTCTTTCGCAAACCAGTTTGCGGGGTGCTGCCGCTGTTCAATATCCGTTTCTGCACCATCTGGACCCTTGTTCTGATGGCgctg   >  1:390658/1‑150 (MQ=255)
                                                                                                                                                    tcGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCGTGGTGCTGCCGCTGTTCAATATCCTGTTCTGCACCATCTGGACCCTTGGCCTGATGGCgctgc  >  2:798468/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
ACCAAACAGGACGACGCTCAGGCCCAGTTGCGTGCCTATCACCAGGTCCGCGAGATCATCGCGAAGTACCAGGCCGCCCACCCGGAGTTCACCGATGAGGTGTTCATCGCCGGCATGCCGATCTTCATTGCGGCCCAGCAGGAAATCATCGATCACGACCTGGCGATGCTGTTCCCGATCGTGTTCCTGTTGGTGACCTCTCTGCTGGTTTTCTTCTTTCGCAAACCACTTGGCGTGGTGCTGCCGCTGTTCAATATCCTGTTCTGCACCATCTGGACCCTTGGCCTGATGGCGCTGC  >  NC_002947/3175445‑3175742

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