Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,735,764 G→T 6.1% P83Q (CCG→CAG)  rluD ← 23S rRNA pseudouridine(1911/1915/1917) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,735,7640GT 6.1% 325.9 / 6.0 97P83Q (CCG→CAG) rluD23S rRNA pseudouridine(1911/1915/1917) synthase
Reads supporting (aligned to +/- strand):  ref base G (44/47);  new base T (4/2);  total (48/49)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

GCACGCGGTACATCGGCAATGGGTGGATAGTAATGAAGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGATCCATTCTTTTATTCGCGAA  >  NC_000913/2735617‑2735913
                                                                                                                                                   |                                                                                                                                                     
gCACGCGGTACATCGGCAATGGGTGGATAGTAATGAAGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCggg                                                                                                                                                     >  1:1776090/1‑150 (MQ=255)
 cACGCGGTACATCGGCAATGGGTGGATAGTAATGAAGCAACGCATTCAGTACCGTGCCATCCGTGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGa                                                                                                                                                    <  1:1963310/150‑1 (MQ=255)
      ggTACATCGGCAATGGGTGGATAGTAATGAAGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATcc                                                                                                                                               <  1:1993852/150‑1 (MQ=255)
         aCATCGGCAATGGGTGGATAGTAATGAAGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCCGCTGGATATcgtct                                                                                                                                            <  2:740459/150‑5 (MQ=255)
         aCATCGGCAATGGGTGGATAGTAATGAAGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGc                                                                                                                                            <  2:1274203/150‑1 (MQ=255)
         aCATCGGCAATGGGTGGATAGTAATGAAGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGc                                                                                                                                            <  1:1038307/150‑1 (MQ=255)
           aTCGGCAATGGGTGGATAGTAATGAAGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCgg                                                                                                                                          >  1:1467846/1‑150 (MQ=255)
              ggCAATGGGTGGATAGTAATGAAGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTc                                                                                                                                       <  1:1829086/150‑1 (MQ=255)
              ggCAATGGGTGGATAGTAATGAAGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTc                                                                                                                                       <  1:660039/150‑1 (MQ=255)
               gCAATGGGTGGATAGTAATGAAGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCa                                                                                                                                      >  1:928494/1‑150 (MQ=255)
                   tGGGTGGATAGTAATGAAGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAAc                                                                                                                                  <  1:1514789/150‑1 (MQ=255)
                      gTGGATAGTAATGAGGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAAcgcg                                                                                                                               <  1:385116/150‑1 (MQ=255)
                        ggATAGTAATGAAGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGct                                                                                                                             >  2:1448660/1‑150 (MQ=255)
                        ggATAGTAATGAAGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAAAGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGct                                                                                                                             <  1:591950/150‑1 (MQ=255)
                         gATAGTAATGAAGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGctt                                                                                                                            >  1:1714046/1‑150 (MQ=255)
                           tAGTAATGAAGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGcttct                                                                                                                          <  2:1575556/150‑1 (MQ=255)
                            aGTAATGAAGCAACGCATTCAGTACCGTGCCATCCCGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTAATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGcttctt                                                                                                                         >  2:1645149/1‑150 (MQ=255)
                             gTAATGAAGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGcttcttc                                                                                                                        >  1:206406/1‑150 (MQ=255)
                              tAATGAAGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGcttcttct                                                                                                                       <  2:6814/150‑1 (MQ=255)
                               aaTGAAGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGcttcttctt                                                                                                                      >  1:347340/1‑150 (MQ=255)
                               aaTGAAGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGcttcttctt                                                                                                                      >  2:703433/1‑150 (MQ=255)
                                aTGAAGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCttcttcttc                                                                                                                     <  1:648494/150‑1 (MQ=255)
                                aTGAAGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCttcttcttc                                                                                                                     <  2:1297230/150‑1 (MQ=255)
                                  gAAGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCaa                                                                                                                   >  2:1685328/1‑150 (MQ=255)
                                    aGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAAtc                                                                                                                 >  1:780473/1‑150 (MQ=255)
                                    aGCAACACATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAAtc                                                                                                                 <  2:316239/150‑1 (MQ=255)
                                     gCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAAtct                                                                                                                <  2:206406/150‑1 (MQ=255)
                                          gCAATCAGTACCGTGCCATCCGGGATACCCGCGCCAGGATGTAACACCAGGTAGCGCGGTTTATTAAAGATAAAAATGTCATCATCTTCATAGACGATATCCAGCGGCATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCg                                                                                                           >  2:1318287/1‑150 (MQ=255)
                                            aTTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGtt                                                                                                         <  2:1776090/150‑1 (MQ=255)
                                            aTTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGtt                                                                                                         <  1:363505/150‑1 (MQ=255)
                                             ttCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTg                                                                                                        <  1:1490062/150‑1 (MQ=255)
                                              tCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGa                                                                                                       <  2:37855/150‑1 (MQ=255)
                                               cAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGAt                                                                                                      >  1:438761/1‑150 (MQ=255)
                                                aGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATg                                                                                                     <  2:1778686/150‑1 (MQ=255)
                                                   aCCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCGATCTCAGCGTTGATGGCa                                                                                                  <  2:1714046/150‑1 (MQ=255)
                                                      gTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAAcc                                                                                               <  1:703433/150‑1 (MQ=255)
                                                      gTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAAcc                                                                                               <  2:547291/150‑1 (MQ=255)
                                                          cATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCt                                                                                           >  2:1396817/1‑150 (MQ=255)
                                                          cATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCt                                                                                           >  2:914577/1‑150 (MQ=255)
                                                              cGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGcc                                                                                       >  2:862223/1‑150 (MQ=255)
                                                                   tACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCa                                                                                  >  1:241756/1‑150 (MQ=255)
                                                                       cgcgCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATAc                                                                              <  1:84184/150‑1 (MQ=255)
                                                                                  gTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCg                                                                   <  2:914023/150‑1 (MQ=255)
                                                                                    accaccAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGc                                                                 >  1:864128/1‑150 (MQ=255)
                                                                                     ccaccaGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCt                                                                >  1:1514565/1‑150 (MQ=255)
                                                                                     ccaccaGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCt                                                                <  1:1685328/150‑1 (MQ=255)
                                                                                     ccaccaGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCt                                                                >  1:767474/1‑150 (MQ=255)
                                                                                           ggTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATcac                                                          >  1:251693/1‑150 (MQ=255)
                                                                                           ggTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATcac                                                          >  2:1414022/1‑150 (MQ=255)
                                                                                             tCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACaa                                                        >  1:602111/1‑150 (MQ=255)
                                                                                                   ggTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTg                                                  <  1:115948/150‑1 (MQ=255)
                                                                                                    gTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGc                                                 <  2:241756/150‑1 (MQ=255)
                                                                                                    gTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGc                                                 <  2:928494/150‑1 (MQ=255)
                                                                                                     tttATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACTCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGACAACCAATACTTTTTCTTTAGGCTTATAACAAAATTTgac                                                >  1:1691553/1‑148 (MQ=255)
                                                                                                     tttATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGcc                                                <  2:204423/150‑1 (MQ=255)
                                                                                                       tattaATGATAATAATGTCTTCATCTTCGTAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGt                                              >  1:1712511/1‑150 (MQ=255)
                                                                                                       tattaATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGt                                              >  1:1024682/1‑150 (MQ=255)
                                                                                                       tattaATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGt                                              >  2:427667/1‑150 (MQ=255)
                                                                                                       tattaATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGt                                              <  2:864128/150‑1 (MQ=255)
                                                                                                         ttaATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTa                                            >  2:1461803/1‑150 (MQ=255)
                                                                                                          taATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATc                                                            >  1:200835/1‑133 (MQ=255)
                                                                                                          taATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATc                                                            <  2:200835/133‑1 (MQ=255)
                                                                                                             tGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCa                                        <  1:29931/150‑1 (MQ=255)
                                                                                                             tGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCa                                        >  1:654553/1‑150 (MQ=255)
                                                                                                             tGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCa                                        >  2:627414/1‑150 (MQ=255)
                                                                                                                     atGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGc                                >  2:355703/1‑150 (MQ=255)
                                                                                                                     atGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGc                                <  2:49393/150‑1 (MQ=255)
                                                                                                                     atGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGc                                <  1:355703/150‑1 (MQ=255)
                                                                                                                      tGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCt                               <  1:427667/150‑1 (MQ=255)
                                                                                                                         cttcatcttcatAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGt                            >  2:1367381/1‑150 (MQ=255)
                                                                                                                         cttcatcttcatAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGt                            >  2:83808/1‑150 (MQ=255)
                                                                                                                            catcttcatAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGa                         <  1:607287/150‑1 (MQ=255)
                                                                                                                             atcttcatAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAg                        <  2:251693/150‑1 (MQ=255)
                                                                                                                               cttcatAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGa                      <  1:1703928/150‑1 (MQ=255)
                                                                                                                               cttcatAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGa                      >  1:1616196/1‑150 (MQ=255)
                                                                                                                                ttcatAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGAt                     >  1:1829174/1‑150 (MQ=255)
                                                                                                                                  catAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGATcc                   >  2:1487809/1‑150 (MQ=255)
                                                                                                                                     gcaggATATCCCGCTGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGATCCa                  >  1:1247972/5‑148 (MQ=255)
                                                                                                                                     gcaggATATCCCGCTGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGATCCa                  <  2:1247972/144‑1 (MQ=255)
                                                                                                                                     gcaggATATCCCGCTGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCACTGGTCGAGGATCCAtt                >  1:1618361/5‑150 (MQ=255)
                                                                                                                                     aGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGGCGTTAACCAGCACTCGCTGGTCGAGGATCCAtt                >  1:1055430/1‑150 (MQ=255)
                                                                                                                                     aGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTAGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGATCCAtt                >  1:1558698/1‑150 (MQ=255)
                                                                                                                                      gACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGATCCATTc               >  2:222488/1‑150 (MQ=255)
                                                                                                                                        ggATATCCCGCTGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGATCCATTCtt             >  1:198583/2‑150 (MQ=255)
                                                                                                                                        ggATATCCCGCTGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGATCCATTCtt             >  2:738762/2‑150 (MQ=255)
                                                                                                                                         gATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGATCCATTCttt            <  1:1356996/150‑1 (MQ=255)
                                                                                                                                         gATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGAGTTGATGTCAACCTGCTCGCCACCCAATAGTTTTTCTTTCGGCTTTTCACAAACTTTGCCGTTAACGAGCACTCGCTGGTCGAGGATCCATTCttt            <  1:1645149/150‑1 (MQ=255)
                                                                                                                                           tatCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGATCCATTCTTTTa          <  2:198047/150‑1 (MQ=255)
                                                                                                                                           tatCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGATCCATTCTTTTa          <  2:1467846/150‑1 (MQ=255)
                                                                                                                                            atCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGATCCATTCTTTTAt         >  2:1991250/1‑150 (MQ=255)
                                                                                                                                            atCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGATCCATTCTTTTAt         <  1:1414022/150‑1 (MQ=255)
                                                                                                                                             tCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGATCCATTCTTTTAtt        <  2:323832/150‑1 (MQ=255)
                                                                                                                                              ccAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCACCCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGATCCATTCTTTTATTc       >  1:1980438/1‑150 (MQ=255)
                                                                                                                                              ccAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGATCCATTCTTTTATTc       >  2:1919681/1‑150 (MQ=255)
                                                                                                                                               cAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGATCCATTCTTTTATTcg      >  1:852393/1‑150 (MQ=255)
                                                                                                                                               cAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGATCCATGCTTTTATTcg      >  1:548800/1‑150 (MQ=255)
                                                                                                                                                aGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGATCCATTCTTTTATTcgc     <  2:780473/150‑1 (MQ=255)
                                                                                                                                                 gCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGATCCATTCTTTTATTcgcg    <  1:1396817/150‑1 (MQ=255)
                                                                                                                                                   gggATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGATCCATTCTTTTATTCGCGaa  <  1:1919681/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCACGCGGTACATCGGCAATGGGTGGATAGTAATGAAGCAACGCATTCAGTACCGTGCCATCCGGGTTACCCGCGCCAGGATGTACCACCAGGTCGCGCGGTTTATTAATGATAATAATGTCTTCATCTTCATAGACGATATCCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTCTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGATCCATTCTTTTATTCGCGAA  >  NC_000913/2735617‑2735913

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