Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 546,893 G→T 8.5% G81C (GGC→TGC)  rpsG → 30S ribosomal protein S7

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947546,8930GT 8.5% 237.5 / 6.8 79G81C (GGC→TGC) rpsG30S ribosomal protein S7
Reads supporting (aligned to +/- strand):  ref base G (38/34);  new base T (3/4);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGTCGATGGCTCTGCGTCTGGCCGGCGAACTGCTGGATGCTGCCGAAG  >  NC_002947/546748‑547040
                                                                                                                                                 |                                                                                                                                                   
aaGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTcggcgg                                                                                                                                                 <  1:223909/150‑1 (MQ=255)
 aGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGTAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCTGCGGt                                                                                                                                                >  2:655993/1‑150 (MQ=255)
     gCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCa                                                                                                                                            <  1:519965/150‑1 (MQ=255)
          aagGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTAc                                                                                                                                       <  1:1367959/150‑1 (MQ=255)
          aagGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTAc                                                                                                                                       <  1:632126/150‑1 (MQ=255)
            ggCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCTCTGAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCa                                                                                                                                     <  1:919765/150‑1 (MQ=255)
            ggCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAGGTAACGTAAGCGCGTGTCGGCGGGGCCACGTACCa                                                                                                                                     >  2:1586569/1‑150 (MQ=255)
              ccGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGTTAATGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTAGGCGGGGCCACTGACCAgg                                                                                                                                   >  1:920345/1‑150 (MQ=255)
                        gcgcATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGa                                                                                                                         >  1:714645/1‑150 (MQ=255)
                           cATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGTAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCTGCGGTGCCACTTACCAGGTCCCGGTTGAAGt                                                                                                                      <  1:655993/150‑1 (MQ=255)
                              cGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCg                                                                                                                   <  1:181260/150‑1 (MQ=255)
                                tttACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGCGACAGCTCTCGACGCCATCGCTCCGCTGGTCGCAGTAAAGTCCCGTCGTGTCTGCGGCGCCCCTTACCAGGTCACGGTTGAAGATCGTc                                                                                                                 >  2:89986/1‑150 (MQ=255)
                                    cGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATc                                                                                                             <  1:596348/150‑1 (MQ=255)
                                     ggTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATcc                                                                                                            >  1:695501/1‑150 (MQ=255)
                                        gCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCcgt                                                                                                         <  2:714645/150‑1 (MQ=255)
                                           cTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCcgtcgt                                                                                                      >  1:134519/1‑150 (MQ=255)
                                                 aCCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGct                                                                                                >  1:1500701/1‑150 (MQ=255)
                                                 aCCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACAAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGct                                                                                                >  1:310749/1‑150 (MQ=255)
                                                    gTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGATAAAGCTATCAACGTGATGGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTg                                                                                             <  2:793099/150‑1 (MQ=255)
                                                    gTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTg                                                                                             <  1:1554763/150‑1 (MQ=255)
                                                    gTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTg                                                                                             >  2:1554763/1‑150 (MQ=255)
                                                         aGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCAc                                                                                                                                           <  2:911402/99‑1 (MQ=255)
                                                         aGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCAc                                                                                                                                           >  1:911402/1‑99 (MQ=255)
                                                         aGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGt                                                                                                                  <  2:819320/124‑1 (MQ=255)
                                                         aGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGt                                                                                                                  >  1:819320/1‑124 (MQ=255)
                                                         aGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACTCTCTGGCAAt                                                                                        <  2:920345/150‑1 (MQ=255)
                                                                aaGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTgg                                                                                 >  2:1567659/1‑150 (MQ=255)
                                                                  gAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCt                                                                               >  1:1076195/1‑150 (MQ=255)
                                                                      aGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTa                                                                           <  1:1567659/150‑1 (MQ=255)
                                                                        cGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCTGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGa                                                                         <  2:891415/150‑1 (MQ=255)
                                                                           cccccTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCTGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTa                                                                      >  1:509470/1‑150 (MQ=255)
                                                                              ccTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATc                                                                                                             <  1:1481018/108‑1 (MQ=255)
                                                                              ccTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATc                                                                                                             >  2:1481018/1‑108 (MQ=255)
                                                                              ccTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGCCCATCCCGTCGTAACGCCCTGGCAACGCGCTGGCTCGTACACTCCGc                                                                   >  2:817590/1‑150 (MQ=255)
                                                                              ccTGAAAACCTTCAAGAAAGCTCTCAACGCCATCGCTCCGCTGCTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGc                                                                   >  2:881467/1‑150 (MQ=255)
                                                                                   aaaTCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCa                                                              >  2:864895/1‑150 (MQ=255)
                                                                                     aTCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAg                                                            >  1:452712/1‑150 (MQ=255)
                                                                                           gagaAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCgcggc                                                      >  2:930789/1‑150 (MQ=255)
                                                                                             gaAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCTGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCga                                                    <  1:89986/150‑1 (MQ=255)
                                                                                             gaAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTgg                                                                                            <  1:79428/110‑1 (MQ=255)
                                                                                             gaAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTgg                                                                                            >  2:79428/1‑110 (MQ=255)
                                                                                             gaAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCga                                                    >  2:523462/1‑150 (MQ=255)
                                                                                                aGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGaa                                                 <  2:1076195/150‑1 (MQ=255)
                                                                                                  ctctCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCACGTCGTAAAGCTGTGGCAAGGCGCTGGCTTGTAGACTACGCCCGCAAGCGCGGCGAGAAGt                                               >  1:773969/1‑150 (MQ=255)
                                                                                                   tctcGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGTc                                              <  1:1586569/150‑1 (MQ=255)
                                                                                                   tctcGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGTc                                              >  2:1321037/1‑150 (MQ=255)
                                                                                                       gACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCTGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGTCGATg                                          <  2:509470/150‑1 (MQ=255)
                                                                                                       gACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTc                                                                                                                 >  1:733412/1‑79 (MQ=255)
                                                                                                       gACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTc                                                                                                                 <  2:733412/79‑1 (MQ=255)
                                                                                                          gCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGACGATAAGTCGaatgct                                       >  2:913666/1‑145 (MQ=255)
                                                                                                           ccATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTcc                                                                                                                <  1:334708/76‑1 (MQ=255)
                                                                                                           ccATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTcc                                                                                                                >  2:334708/1‑76 (MQ=255)
                                                                                                           ccATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGTCGATGGctc                                      >  1:423740/1‑150 (MQ=255)
                                                                                                             aTCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGt                                               >  2:254609/1‑139 (MQ=255)
                                                                                                             aTCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGt                                               <  1:254609/139‑1 (MQ=255)
                                                                                                                 cTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGTCGa                                            >  2:498723/1‑138 (MQ=255)
                                                                                                                 cTCCGATGGTGGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGTCGa                                            <  1:498723/138‑1 (MQ=255)
                                                                                                                  tCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAAt                                                                                        <  2:366746/93‑1 (MQ=255)
                                                                                                                  tCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAAt                                                                                        >  1:366746/1‑93 (MQ=255)
                                                                                                                  tCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGTCGATGGCTCTGCGTCt                               >  1:235957/1‑150 (MQ=255)
                                                                                                                   cctttGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGTCGATGGCTCTGCGTCTg                              <  1:881467/146‑1 (MQ=255)
                                                                                                                     gCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCcgtc                                                                                                        >  1:1186944/1‑74 (MQ=255)
                                                                                                                     gCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCcgtc                                                                                                        <  2:1186944/74‑1 (MQ=255)
                                                                                                                         gTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGTCGATGGCTCTGCGTCTGGCCGGc                        >  2:935398/1‑150 (MQ=255)
                                                                                                                            gAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCt                                                                                   <  2:341545/88‑1 (MQ=255)
                                                                                                                            gAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCt                                                                                   >  1:341545/1‑88 (MQ=255)
                                                                                                                                tAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGTCGATGGCTCTGCGTCTGGCCGGCGAActgc                 >  2:640678/1‑150 (MQ=255)
                                                                                                                                 aaaGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGTCGATGGCTCTGCGTCTGGCCGGCGAActgct                <  2:235957/150‑1 (MQ=255)
                                                                                                                                  aaGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGCAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGTCGCTGGCTCTGCGTCTGGTCGGCGAActgctg               <  1:817590/150‑1 (MQ=255)
                                                                                                                                  aaGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATc                                                                                                             >  2:1643057/1‑56 (MQ=255)
                                                                                                                                  aaGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATc                                                                                                             <  1:1643057/56‑1 (MQ=255)
                                                                                                                                  aaGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGTCGATGGCTCTGCGTCTGGCCGGCGAActgctg               <  2:423740/150‑1 (MQ=255)
                                                                                                                                   aGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGTCGATGGCTCTGCGTCTGGCCGGCGAACTGCTgg              <  2:134519/150‑1 (MQ=255)
                                                                                                                                    gTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGCTCTGGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGTCGATGGCTCTGCGTCTGGCCGGCGAACTGCTGGa             <  1:1251279/150‑1 (MQ=255)
                                                                                                                                        cgtcgtGTCGGCGGTGCCACTTACCAGGTCCCTGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGTCGATGGCTCTGCGTCTGGCCGGCGAACTGCTGGAtgct         >  2:703054/1‑150 (MQ=255)
                                                                                                                                         gtcgtGTCGGCGGTGCCACTTACCAGGTCCCTGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGTCGATGGCTCTGCGTCTGGCCGGCGAACTGCTGGAtgctg        <  1:703054/150‑1 (MQ=255)
                                                                                                                                              gtCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGTCGATGGCTCTGCGTCTGGCCGGCGAACTACTGGATGCTGCCGaa   >  1:1272040/1‑150 (MQ=255)
                                                                                                                                               tCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGTCGATGGCTCTGc                                   >  2:375964/1‑117 (MQ=255)
                                                                                                                                               tCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGTCGATGGCTCTGc                                   <  1:375964/117‑1 (MQ=255)
                                                                                                                                               tCGGCGGTGCCACTTACCAGGTCCAGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGTCGATGGCTCTGCGTCTGGCCGGCGAACTGCTGGATGCTGCCGAAg  <  2:695501/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
AAGCGGCAAGAAGGCCGTAGCCGAGCGCATCGTTTACGGTGCCCTGGATACCGTCAAAGCACGCAAGAACAGCGACCCCCTGGAAATCTTCGAGAAAGCTCTCGACGCCATCGCTCCGCTGGTCGAAGTAAAGTCCCGTCGTGTCGGCGGTGCCACTTACCAGGTCCCGGTTGAAGTTCGTCCATCCCGTCGTAACGCTCTGGCAATGCGCTGGCTCGTAGACTACGCCCGCAAGCGCGGCGAGAAGTCGATGGCTCTGCGTCTGGCCGGCGAACTGCTGGATGCTGCCGAAG  >  NC_002947/546748‑547040

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