Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,039,701 A→C 6.3% I284S (ATC→AGC) ‡ PP_4446 ← group II intron‑encoding maturase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,039,7010AC 6.3% 265.3 / 8.1 95I284S (ATC→AGC) ‡PP_4446group II intron‑encoding maturase
Reads supporting (aligned to +/- strand):  ref base A (39/50);  new base C (2/4);  total (41/54)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

GATTACGCTGGCTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGCCCAGGAAGTCAAACGTCT  >  NC_002947/5039552‑5039848
                                                                                                                                                     |                                                                                                                                                   
gATTACGCTGGCTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGa                                                                                                                                                     <  1:383508/150‑1 (MQ=255)
   tACGCTGGCTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATgg                                                                                                                                                  <  1:1907313/150‑1 (MQ=255)
   tACGCTGGCTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATgg                                                                                                                                                  <  1:1849178/150‑1 (MQ=255)
      gCTGGCTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGCTTCCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGtt                                                                                                                                               <  2:1662459/150‑1 (MQ=255)
        tGGCTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCATCTGGGTTCa                                                                                                                                             >  1:1402434/1‑150 (MQ=255)
         ggCTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCAt                                                                                                                                            <  1:241717/150‑1 (MQ=255)
         ggCTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTGTCCCGGCACCGCTTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCAt                                                                                                                                            <  2:1662675/150‑1 (MQ=255)
         ggCCACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCATCTGGGTTCaa                                                                                                                                            <  2:1402434/150‑2 (MQ=255)
           ttACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCATCTGGGTTCaagg                                                                                                                                          <  1:402515/149‑4 (MQ=255)
           cTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCATCTGGGTTCaagg                                                                                                                                          <  1:1529889/150‑4 (MQ=255)
           cTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGa                                                                                                                                          <  2:252647/150‑1 (MQ=255)
           cTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGa                                                                                                                                          <  2:443871/150‑1 (MQ=255)
               gacgTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGt                                                                                                                                                <  2:1801923/140‑1 (MQ=255)
               gacgTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGt                                                                                                                                                >  1:1801923/1‑140 (MQ=255)
               gacgTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTc                                                                                                                                      >  1:824915/1‑150 (MQ=255)
                acgTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCATTTAGACGTGTCAGTTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCg                                                                                                                                     <  1:1346750/150‑1 (MQ=255)
                acgTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCg                                                                                                                                     <  2:994114/150‑1 (MQ=255)
                 cgTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCgg                                                                                                                                    >  2:54614/1‑150 (MQ=255)
                  gTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCa                                                                                                                                             >  2:755238/1‑140 (MQ=255)
                  gTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCa                                                                                                                                             <  1:755238/140‑1 (MQ=255)
                  gTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGgt                                                                                                                                   <  2:2281240/150‑1 (MQ=255)
                   tttGAGGACTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGgtg                                                                                                                                  >  1:1415293/1‑150 (MQ=255)
                         ggCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGgtgt                                                                                                                                 >  1:988702/1‑145 (MQ=255)
                         ggCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGgtgt                                                                                                                                 <  2:988702/145‑1 (MQ=255)
                          gCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGc                                                                                                                           <  2:793901/150‑1 (MQ=255)
                          gCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGc                                                                                                                           <  2:517830/150‑1 (MQ=255)
                               ccgccATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATc                                                                                                                      >  1:638024/1‑150 (MQ=255)
                                cgccATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATcc                                                                                                                     <  1:54614/150‑1 (MQ=255)
                                   cATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCATCTCATCCCGc                                                                                                                  <  2:2279362/150‑1 (MQ=255)
                                   cATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGc                                                                                                                  <  2:1951344/150‑1 (MQ=255)
                                    aTAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCa                                                                                                                 >  2:1372450/1‑150 (MQ=255)
                                    aTAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCa                                                                                                                 >  1:253813/1‑150 (MQ=255)
                                     tAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCAt                                                                                                                <  1:714829/150‑1 (MQ=255)
                                      aaaCGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATg                                                                                                               >  1:985973/1‑150 (MQ=255)
                                      aaaCGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTACCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATg                                                                                                               <  1:270995/150‑1 (MQ=255)
                                          ggCAGACCGCCGAACGAAAACCACCCAGACGTGACAGGTTTCACGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGCCGTCGACGGTTCTGGTCATCCCGCATGGccc                                                                                                           >  2:552251/1‑149 (MQ=255)
                                             aGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTTAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCcagca                                                                                                        <  1:458456/150‑1 (MQ=255)
                                             aGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCcagca                                                                                                        <  1:736394/150‑1 (MQ=255)
                                                ccgccgAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCACGGCCAgcagca                                                                                                     <  2:1055040/150‑1 (MQ=255)
                                                 cgccgAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCGGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCGGCGCATCCCGCATGGCCAgcggcac                                                                                                    >  2:1192480/1‑149 (MQ=255)
                                                 cgccgAACGAAAACCACCCAGACATGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCAGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCGGCACATCCCGCATGGCCAGCTGCAt                                                                                                    >  2:1192376/1‑150 (MQ=255)
                                                           aaaCCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGc                                                                                                                           <  2:1090703/117‑1 (MQ=255)
                                                           aaaCCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGc                                                                                                                           >  1:1090703/1‑117 (MQ=255)
                                                            aaCCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTc                                                                                         >  1:2205124/1‑150 (MQ=255)
                                                              ctaactaGACATGTCAGGTTTACCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAACCATTGGCCTTGAACCCAGATGGGTTCATGACATCGGTGTTTCAGCTCATCCCGAATGGCCAGCAGCATCGCACGCATTCGc                                                                                       >  2:1140911/7‑150 (MQ=255)
                                                               caccCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGct                                                                                      >  1:1349727/1‑150 (MQ=255)
                                                               caccCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGct                                                                                      >  2:1750185/1‑150 (MQ=255)
                                                                 cccAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGcttc                                                                                    <  1:1521938/150‑1 (MQ=255)
                                                                 cccAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGcttc                                                                                    <  2:1349727/150‑1 (MQ=255)
                                                                  ccAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGc                                                                                                                           >  2:2151297/1‑110 (MQ=255)
                                                                  ccAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGc                                                                                                                           <  1:2151297/110‑1 (MQ=255)
                                                                       cGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGACCCCTCATGGGTTCATGACGTCGGTGTTTCCGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGAcc                                                                              <  1:1118936/150‑1 (MQ=255)
                                                                         tgtCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGt                                                                            <  1:2290937/150‑1 (MQ=255)
                                                                           tCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCATCTGGGTTCaa                                                                                                                                            >  2:2340878/1‑83 (MQ=255)
                                                                           tCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCATCTGGGTTCaa                                                                                                                                            <  1:2340878/84‑2 (MQ=255)
                                                                           tCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGtt                                                                                                                                               >  2:2335686/1‑81 (MQ=255)
                                                                           tCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGtt                                                                                                                                               <  1:2335686/81‑1 (MQ=255)
                                                                           tCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCa                                                                          >  2:558383/1‑150 (MQ=255)
                                                                               gTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCg                                                                      <  2:2453308/150‑1 (MQ=255)
                                                                               gTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCg                                                                      >  2:2109215/1‑150 (MQ=255)
                                                                                    ccGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTa                                                                 >  1:1858842/1‑150 (MQ=255)
                                                                                    ccGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCAAGCATTCGCTTCTTGACCGTCAGTCGCAGTa                                                                 >  1:1183088/1‑150 (MQ=255)
                                                                                          ccGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCTCTTCTTGACCGTCAGTCGCAGTATTTGaa                                                           <  2:1415293/150‑1 (MQ=255)
                                                                                           cgcgTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGaaa                                                          <  2:1183088/150‑1 (MQ=255)
                                                                                           cgcgTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGaaa                                                          <  1:1372450/150‑1 (MQ=255)
                                                                                               tGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCa                                                      <  1:568066/150‑1 (MQ=255)
                                                                                                  tAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATccc                                                                                                                    >  2:1671939/1‑85 (MQ=255)
                                                                                                  tAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATccc                                                                                                                    <  1:1671939/85‑1 (MQ=255)
                                                                                                  tAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCg                                                   >  1:716042/1‑150 (MQ=255)
                                                                                                    aTTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCt                                                 >  1:2310923/1‑150 (MQ=255)
                                                                                                        aaGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTc                                                                       <  2:781668/124‑1 (MQ=255)
                                                                                                        aaGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTc                                                                       >  1:781668/1‑124 (MQ=255)
                                                                                                            aaCCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTgg                                         <  1:558383/150‑1 (MQ=255)
                                                                                                            aaCCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTgg                                         <  2:2205124/150‑1 (MQ=255)
                                                                                                            aaCCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCACTTCGTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTgg                                         <  1:552251/150‑1 (MQ=255)
                                                                                                                  tCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGtactac                                   <  2:253813/150‑1 (MQ=255)
                                                                                                                  tCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGtactac                                   <  2:2305694/150‑1 (MQ=255)
                                                                                                                    accaccCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTacaa                                 >  1:2134109/1‑150 (MQ=255)
                                                                                                                       accCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAg                              >  1:933973/1‑150 (MQ=255)
                                                                                                                         ccGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTg                            >  1:760420/1‑150 (MQ=255)
                                                                                                                         ccGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTg                            <  2:1858842/150‑1 (MQ=255)
                                                                                                                          cGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGc                           <  1:1325699/150‑1 (MQ=255)
                                                                                                                               gAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGaa                      <  2:1931782/150‑1 (MQ=255)
                                                                                                                                 gCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCa                                                                                                                             >  1:801139/1‑45 (MQ=255)
                                                                                                                                 gCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCa                                                                                                                             <  2:801139/45‑1 (MQ=255)
                                                                                                                                  ccATTGCCCTTGAACCCAGATGGGTCCATGACGTTGGTTGTTCATCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACTGTCAGTCGCAGTATTTGAAACCTACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGcc                   >  1:655178/1‑150 (MQ=255)
                                                                                                                                    aTTGCCCTTGAACCCAGATGGGTTCATGACGTCGTTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGc                                                                                       >  2:706866/1‑80 (MQ=255)
                                                                                                                                    aTTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGc                                                                                       <  1:706866/80‑1 (MQ=255)
                                                                                                                                    aTTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTacaa                                 >  1:1832771/1‑134 (MQ=255)
                                                                                                                                    aTTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTacaa                                 <  2:1832771/134‑1 (MQ=255)
                                                                                                                                    aTTGCCCTTGAACCCAGATGGATTCATAACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAACATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCTCTTCTGTTGGTACTACAACAGTGCGTGAAGCCCa                 <  1:1140911/150‑1 (MQ=255)
                                                                                                                                        cccTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGCCCAGGaa             <  2:1940047/150‑1 (MQ=255)
                                                                                                                                        cccTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGCCCAGGaa             >  1:241590/1‑150 (MQ=255)
                                                                                                                                        cccTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGCCCAGGaa             >  1:2151198/1‑150 (MQ=255)
                                                                                                                                          cTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGCCCAGGAAGt           >  1:544879/1‑150 (MQ=255)
                                                                                                                                                   aGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGCCCAGGAAGTCAAACGtct  >  2:1161193/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GATTACGCTGGCTACGACGTTTGAGGGCTTGCCGCCATAAACGGCAGACCGCCGAACGAAAACCACCCAGACGTGTCAGGTTTCCCGGCACCGCGTGGTAATTGAAGTAACCGCTCACCACCCGTGTGAGCCATTGCCCTTGAACCCAGATGGGTTCATGACGTCGGTGTTTCAGCTCATCCCGCATGGCCAGCAGCATCGCACGCATTCGCTTCTTGACCGTCAGTCGCAGTATTTGAAACCCACCGCTTCTGTTGGTACTACAACAGTGCGTGAAGCCCAGGAAGTCAAACGTCT  >  NC_002947/5039552‑5039848

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