Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,194,811 G→T 5.3% P47P (CCC→CCA rluE ← 23S rRNA pseudouridine(2457) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,194,8110GT 5.3% 283.8 / 2.8 85P47P (CCC→CCArluE23S rRNA pseudouridine(2457) synthase
Reads supporting (aligned to +/- strand):  ref base G (49/31);  new base T (2/3);  total (51/34)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCC  >  NC_000913/1194663‑1194958
                                                                                                                                                    |                                                                                                                                                   
gcgcCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAg                                                                                                                                                     >  2:298155/1‑149 (MQ=255)
    cTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGtt                                                                                                                                                 >  8:416187/1‑149 (MQ=255)
      gCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTa                                                                                                                                               >  3:178675/1‑149 (MQ=255)
       cAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTATGGTTTAt                                                                                                                                              >  8:208848/1‑149 (MQ=255)
           gcgcCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGaa                                                                                                                                          <  4:178675/149‑1 (MQ=255)
            cgcCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAAc                                                                                                                                         <  5:190831/149‑1 (MQ=255)
             gcCGTTATTGGTCAGCACCAGCAAGCCTTCGCTATCGCGGTCAAGGCGACCTGCCGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGaa                                                                                                                                          <  4:259441/147‑1 (MQ=255)
                 ttATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGAt                                                                                                                                    <  2:245192/149‑1 (MQ=255)
                      ggTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGTCCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATcaca                                                                                                                                <  7:416187/148‑1 (MQ=255)
                       gTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTg                                                                                                                                            >  8:381546/1‑135 (MQ=255)
                       gTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTg                                                                                                                                            <  7:381546/135‑1 (MQ=255)
                            cACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACATGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTgg                                                                                                                         >  6:76126/1‑149 (MQ=255)
                                   aaCCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATGAACAGGATCACACGCGTGGGCTGGttt                                                                                                                 >  2:11645/1‑149 (MQ=255)
                                        ttCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGTCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAg                                                                                                             >  3:122793/1‑149 (MQ=255)
                                        ttCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTc                                                                                                               <  1:427537/147‑1 (MQ=255)
                                        ttCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTc                                                                                                               >  2:427537/1‑147 (MQ=255)
                                        ttCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAg                                                                                                             >  3:165342/1‑149 (MQ=255)
                                          cGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAg                                                                                                             <  1:298155/147‑1 (MQ=255)
                                             tATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTa                                                                                                        >  8:130350/1‑149 (MQ=255)
                                               tCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTATGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTa                                                                                                        <  7:208848/147‑1 (MQ=255)
                                                 gcgGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAgg                                                                                                                                                    <  2:424107/101‑1 (MQ=255)
                                                 gcgGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAgg                                                                                                                                                    >  1:424107/1‑101 (MQ=255)
                                                   ggTCAAGGCGACCTGCTGCATAAACACCCTGAACCTGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTg                                                                                                  <  8:55471/149‑1 (MQ=255)
                                                    gTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTgg                                                                                                 <  5:76126/149‑1 (MQ=255)
                                                    gTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTgg                                                                                                 <  3:35865/149‑1 (MQ=255)
                                                    gTCAAGGCGACCTGCTGCATAAACACCCAGAACCGGGATGAATTCTTTTAATGTGTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGCGTTTATTGAACAGGATCACACGCGTGGGCTGCTTTTCAGGTTTACGCCTgg                                                                                                 >  2:223089/1‑149 (MQ=255)
                                                       aaGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGt                                                                                                                   <  7:145007/128‑1 (MQ=255)
                                                       aaGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGt                                                                                                                   >  8:145007/1‑128 (MQ=255)
                                                         ggCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACAcgc                                                                                                                             <  2:434700/116‑1 (MQ=255)
                                                         ggCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACAcgc                                                                                                                             >  1:434700/1‑116 (MQ=255)
                                                                tgctgcATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGc                                                                                      >  7:422438/1‑148 (MQ=255)
                                                                 gctgcATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGc                                                                                      >  2:36397/1‑147 (MQ=255)
                                                                  ctgcATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGa                                                                                    <  8:174940/148‑1 (MQ=255)
                                                                     cATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCt                                                                                 <  7:194123/148‑1 (MQ=255)
                                                                     cATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCt                                                                                 >  2:434682/1‑148 (MQ=255)
                                                                         aaCACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTa                                                                                                        <  3:288063/121‑1 (MQ=255)
                                                                         aaCACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTa                                                                                                        >  4:288063/1‑121 (MQ=255)
                                                                          acacCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCg                                                                           >  1:25592/1‑149 (MQ=255)
                                                                              ccTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGt                                                                                                           <  7:308000/113‑1 (MQ=255)
                                                                              ccTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGt                                                                                                           >  8:308000/1‑113 (MQ=255)
                                                                              ccTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTgg                                                                                                 >  1:459040/1‑123 (MQ=255)
                                                                              ccTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTgg                                                                                                 <  2:459040/123‑1 (MQ=255)
                                                                                tGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAAc                                                                     >  3:212562/1‑149 (MQ=255)
                                                                                       aaaTGAATTCTTTTAATGTTTTGAGTCCGGCTTCATCGGAGAACTGCGGAATAACATCGTAGGGTTAATTCAACAGGATCACACGCGTGGGCTGATTTTCAGGTTTCCGCCTGGTAGAACGATGCGAGATGAATCTCTTAACCTGGTGa                                                              >  4:219738/3‑149 (MQ=255)
                                                                                          tGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGAttt                                                           <  1:87716/149‑1 (MQ=255)
                                                                                          tGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGAttt                                                           >  6:333648/1‑149 (MQ=255)
                                                                                            aaTTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTc                                                          <  1:434682/148‑1 (MQ=255)
                                                                                            aaTTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTc                                                          >  7:288413/1‑148 (MQ=255)
                                                                                             aTTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCt                                                         >  3:37321/1‑148 (MQ=255)
                                                                                              ttCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTa                                                        >  5:180084/1‑148 (MQ=255)
                                                                                                cTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTaaaa                                                     <  1:11645/148‑1 (MQ=255)
                                                                                                 ttttAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGCACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTaa                                                       >  2:25494/1‑146 (MQ=255)
                                                                                                   ttAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTATGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTg                                                                                                  <  5:454657/101‑1 (MQ=255)
                                                                                                   ttAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTATGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTg                                                                                                  >  6:454657/1‑101 (MQ=255)
                                                                                                   ttAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAaagaa                                                  >  8:102439/1‑149 (MQ=255)
                                                                                                    tAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGTACTTCGTATGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAaagaa                                                  <  4:37321/148‑1 (MQ=255)
                                                                                                     aaTGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGt                                                >  7:49302/1‑149 (MQ=255)
                                                                                                      aTGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACAcgc                                                                                                                             <  1:46588/71‑1 (MQ=255)
                                                                                                      aTGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACAcgc                                                                                                                             >  2:46588/1‑71 (MQ=255)
                                                                                                       tGTTTTGCGTCCGGCTGCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGc                                                                          <  7:328761/121‑1 (MQ=255)
                                                                                                       tGTTTTGCGTCCGGCTGCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGc                                                                          >  8:328761/1‑121 (MQ=255)
                                                                                                             gCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCa                                        >  4:322932/1‑149 (MQ=255)
                                                                                                               gTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATg                                      >  8:252860/1‑149 (MQ=255)
                                                                                                                   ggCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGATCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTAt                                  >  5:287430/1‑149 (MQ=255)
                                                                                                                    gCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTAtt                                 >  6:28714/1‑149 (MQ=255)
                                                                                                                       tCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTAtttt                               >  8:266279/1‑148 (MQ=255)
                                                                                                                        cATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTAtttt                               <  3:322932/147‑1 (MQ=255)
                                                                                                                         aTCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTAtttt                               >  4:243865/1‑146 (MQ=255)
                                                                                                                            ggTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTTCATGGTATTTTCAgcga                         >  6:399559/1‑149 (MQ=255)
                                                                                                                               gAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTATTTTCAGAGATTa                      >  4:65060/1‑149 (MQ=255)
                                                                                                                                aaCTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTATTTTCAGAGATTAt                     >  7:127036/1‑149 (MQ=255)
                                                                                                                                aaCTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTATTTTCAGAGATTAt                     <  7:266279/149‑1 (MQ=255)
                                                                                                                                aaCTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTATTTTCAGAGATTAt                     <  8:288413/149‑1 (MQ=255)
                                                                                                                                aaCTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTATTGCATGGTATTTTCAGAGAtt                       >  8:372828/1‑147 (MQ=255)
                                                                                                                                aaCTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGAGGGTTTTCAGGTGTGCGGCTGGTAGATCGTTGCGCGCTGAATCGCTTAACCTGGTGATTTGGAAAAGATGGTTTTTGGGTGGTATTTTCAGAGAGTAt                     >  4:359671/1‑149 (MQ=255)
                                                                                                                                 aCTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTATTTTCAGAGATTAt                     <  5:333648/148‑1 (MQ=255)
                                                                                                                                     cGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTATTTTCAGAGATTATGaa                  <  6:180084/147‑1 (MQ=255)
                                                                                                                                      ggAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGt                                                                >  7:474702/1‑100 (MQ=255)
                                                                                                                                      ggAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGt                                                                <  8:474702/100‑1 (MQ=255)
                                                                                                                                          gaaCATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTATTTTCAGAGATTATGAATTGCCg            >  5:1301/1‑148 (MQ=255)
                                                                                                                                          gaaCATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTATTTTCAGAGATTATGAATTGCCg            >  7:146942/1‑148 (MQ=255)
                                                                                                                                          gaaCATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTATTTTCAGAGATTATGAATTGCCg            >  3:214811/1‑148 (MQ=255)
                                                                                                                                          gaaCATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTATTTTCAGAGATTATGAATTGCCg            >  2:431391/1‑148 (MQ=255)
                                                                                                                                          gaaCATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGACTTCTAAAAGAAGTTTTTTGCATGGTATTTTCAGAGATTATGAATTGCCg            >  6:405491/1‑148 (MQ=255)
                                                                                                                                            aCATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTATTTTCAGAGATTATGAATTGCCGCAt         >  2:268968/1‑149 (MQ=255)
                                                                                                                                               tCGTAGGGTTTTTTGAATAGTATGACTCGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGATCGTTGCGATGTTAATCGATTAACCTGGTGATTTCTAAAAGAATTTTTTTGGTTTTTATTTTCAGAGTTTATGAATTGCGGAATtat      <  8:422438/149‑1 (MQ=255)
                                                                                                                                                 gTAGGGTTTATTGAACAGGAGCACTCGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTTCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTATTTTCAGAGATTATGAATTGCCGCATtata     <  1:223089/148‑1 (MQ=255)
                                                                                                                                                  tAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGc   >  5:17479/1‑149 (MQ=255)
                                                                                                                                                    gggTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGcc  <  1:268968/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCGCCTGCAACGCGCCGTTATTGGTCAGCACCAGCAACCCTTCGCTATCGCGGTCAAGGCGACCTGCTGCATAAACACCCTGAACCGGGATGAATTCTTTTAATGTTTTGCGTCCGGCTTCATCGGTGAACTGCGGAAGAACATCGTAGGGTTTATTGAACAGGATCACACGCGTGGGCTGGTTTTCAGGTTTACGCCTGGTAGAACGTTGCGAGCTGAATCGCTTAACCTGGTGATTTCTAAAAGAAGTTTTTTGCATGGTATTTTCAGAGATTATGAATTGCCGCATTATAGCC  >  NC_000913/1194663‑1194958

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