Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,377,923 A→G 5.2% R95R (CGT→CGC)  rpsI ← 30S ribosomal subunit protein S9

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,377,9230AG 5.2% 235.4 / 2.5 77R95R (CGT→CGC) rpsI30S ribosomal subunit protein S9
Reads supporting (aligned to +/- strand):  ref base A (33/40);  new base G (2/2);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

GGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGTTTCTCAACCATGTCGACCAGTTCCAGCGGCTGACG  >  NC_000913/3377775‑3378063
                                                                                                                                                    |                                                                                                                                            
ggTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGaa                                                                                                                                              <  3:5295/149‑1 (MQ=255)
ggTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGaa                                                                                                                                              <  2:200566/149‑1 (MQ=255)
  ttttttCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACg                                                                                                                                            >  1:339317/1‑149 (MQ=255)
  ttttttCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACg                                                                                                                                            >  6:182775/1‑149 (MQ=255)
      ttCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAg                                                                                                                                         <  2:226975/148‑1 (MQ=255)
       tCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAgg                                                                                                                                        <  1:233796/148‑1 (MQ=255)
       tCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAgg                                                                                                                                        <  5:444952/148‑1 (MQ=255)
          aaaaTTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACt                                                                                                                                    >  2:29961/1‑149 (MQ=255)
              ttGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACtcgtc                                                                                                                                <  2:300957/149‑1 (MQ=255)
                 ttttCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTa                                                                                                                             >  3:106711/1‑149 (MQ=255)
                 ttttCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTa                                                                                                                             >  8:81563/1‑149 (MQ=255)
                    tCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTAtt                                                                                                                           <  5:226958/148‑1 (MQ=255)
                     cTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTc                                                                                                                          <  5:353961/148‑1 (MQ=255)
                        ccGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCAt                                                                                                                       <  8:345316/148‑1 (MQ=255)
                         cGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCa                                                                                                                     <  2:339317/149‑1 (MQ=255)
                         cGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCa                                                                                                                     <  6:61725/149‑1 (MQ=255)
                                gAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTa                                                                                                                             >  2:75732/1‑134 (MQ=255)
                                gAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTa                                                                                                                             <  1:75732/134‑1 (MQ=255)
                                gAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGcgcg                                                                                                              >  5:338952/1‑149 (MQ=255)
                                    ccAATTAACGTTTGGTGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGt                                                                                                           >  6:282317/1‑148 (MQ=255)
                                    ccAATTAACGTTTGGCGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGt                                                                                                           <  5:282317/148‑1 (MQ=255)
                                     cAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGa                                                                                                         <  6:42180/149‑1 (MQ=255)
                                      aaTTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGa                                                                                                         >  3:256493/1‑148 (MQ=255)
                                             gTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGcgcg                                                                                                              <  1:244282/136‑1 (MQ=255)
                                             gTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGcgcg                                                                                                              >  2:244282/1‑136 (MQ=255)
                                             gTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGa                                                                                                 <  2:251951/149‑1 (MQ=255)
                                             gTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGa                                                                                                 >  8:67838/1‑149 (MQ=255)
                                                  gagaACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGa                                                                                             <  4:256493/148‑1 (MQ=255)
                                                  gagaACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGAt                                                                                            <  6:1455/149‑1 (MQ=255)
                                                  gagaACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGAt                                                                                            <  6:443244/149‑1 (MQ=255)
                                                      aCTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCa                                                                                        <  7:155069/149‑1 (MQ=255)
                                                       cTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCAc                                                                                       <  2:450748/149‑1 (MQ=255)
                                                           ggACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTg                                                                                                  >  7:129647/1‑134 (MQ=255)
                                                           ggACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTg                                                                                                  <  8:129647/134‑1 (MQ=255)
                                                                    cgTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCaga                                                                          <  4:106711/149‑1 (MQ=255)
                                                                       gCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAGCGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAgaga                                                                        >  7:186413/1‑148 (MQ=255)
                                                                       gCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAGCGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAgaga                                                                        >  7:192633/1‑148 (MQ=255)
                                                                       gCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAgaga                                                                        >  6:366289/1‑148 (MQ=255)
                                                                         tttACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGCTACCGTGACGGATCGCACCAGCCTGACCAGAGATac                                                                     >  6:265415/1‑149 (MQ=255)
                                                                           tACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATac                                                                     <  7:81563/147‑1 (MQ=255)
                                                                                    cGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCATCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCtt                                                            <  2:80085/147‑1 (MQ=255)
                                                                                    cGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCttt                                                           >  6:365354/1‑148 (MQ=255)
                                                                                    cGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCtt                                                            >  1:80085/1‑147 (MQ=255)
                                                                                      aCTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAAc                                                        >  3:116178/1‑149 (MQ=255)
                                                                                      aCTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAAc                                                        >  5:111451/1‑149 (MQ=255)
                                                                                          tCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGt                                                     <  8:8188/148‑1 (MQ=255)
                                                                                            ttACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCaga                                                                                                                   <  6:112631/84‑1 (MQ=255)
                                                                                            ttACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCaga                                                                                                                   >  5:112631/1‑84 (MQ=255)
                                                                                            ttACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACTGTGAt                                                  >  5:83261/1‑149 (MQ=255)
                                                                                            ttACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGAt                                                  >  2:133318/1‑149 (MQ=255)
                                                                                            ttACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGAt                                                  >  7:537/1‑149 (MQ=255)
                                                                                              aCGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGa                                                   <  7:67838/146‑1 (MQ=255)
                                                                                              aCGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATg                                                 <  8:186821/148‑1 (MQ=255)
                                                                                                gTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTa                                               >  7:74678/1‑148 (MQ=255)
                                                                                                gTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTa                                               >  7:19863/1‑148 (MQ=255)
                                                                                                gTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTa                                               >  7:26664/1‑148 (MQ=255)
                                                                                                          aCGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGtt                                    >  6:447059/1‑149 (MQ=255)
                                                                                                            gAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGTTtc                                  <  5:365354/149‑1 (MQ=255)
                                                                                                                gTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGTTTCTCaa                              <  6:83261/149‑1 (MQ=255)
                                                                                                                gTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGTTTCTCaa                              <  5:265415/149‑1 (MQ=255)
                                                                                                                 tCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGCTCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGTTTCTCAAc                             >  7:449967/1‑149 (MQ=255)
                                                                                                                 tCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGAc                                                                              <  6:213023/100‑1 (MQ=255)
                                                                                                                 tCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGAc                                                                              >  5:213023/1‑100 (MQ=255)
                                                                                                                     gAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCCCCTTTAACGGTGATGTACAGGTCCAGTTTCTCAACCATg                         >  5:453422/1‑149 (MQ=255)
                                                                                                                       gTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGTTTCTCAACCATGTc                       <  8:19863/149‑1 (MQ=255)
                                                                                                                       gTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGTTACTCAACCATGTc                       <  8:26664/149‑1 (MQ=255)
                                                                                                                          aCGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGTTTCTCAACCATGTCGAc                    >  6:4964/1‑149 (MQ=255)
                                                                                                                           cGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGTTTCTCAACCATGTCGAcc                   <  4:116178/149‑1 (MQ=255)
                                                                                                                            gAAGCCAGCTTTACGCAGTTCAGAGCGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGTTTCTCAACCATGTCGACCa                  <  8:186413/149‑1 (MQ=255)
                                                                                                                            gAAGCCAGCTTTACGCAGTTCAGAGCGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGTTTCTCAACCATGTCGACCa                  <  8:192633/149‑1 (MQ=255)
                                                                                                                            gAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGTTTCTCAACCATGTCGACCa                  <  8:74678/149‑1 (MQ=255)
                                                                                                                                    cTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGTTTCTCAACCATGTCGACCAGTTCCAGc          <  5:366289/149‑1 (MQ=255)
                                                                                                                                    cTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGTTTCTCAACCATGTCGACCAGTTCCAGc          <  6:338952/149‑1 (MQ=255)
                                                                                                                                             gTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGTTTCTCAACCATGt                        <  1:416527/126‑1 (MQ=255)
                                                                                                                                             gTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGTTTCTCAACCATGt                        >  2:416527/1‑126 (MQ=255)
                                                                                                                                             gTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGTTTCTCAACCATGTCGACCAGTTCCAGCGGCTGACg  >  7:250924/1‑148 (MQ=255)
                                                                                                                                              ttCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGTTTCTCAACCATGTCGACCAGTTCCAGCGGCt      >  3:78318/1‑143 (MQ=255)
                                                                                                                                              ttCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGTTTCTCAACCATGTCGACCAGTTCCAGCGGCt      <  4:78318/143‑1 (MQ=255)
                                                                                                                                              ttCAGAACGCAGGGACTCGTCGTAATCCAGCAGAGCGCGGGGGATACCGTGACGGATCGCACCAGCCTGACCAGAGTTACCACCCCCTTTCTCGGGGATGTACAGGTCCAGTTTCTCAACCATGTCGACCAGTTCCAGCGGCTGACg  <  6:409450/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
GGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTCAGAACGCAGGGACTCGTCGTATTCCATCAGAGCGCGGGTGATACCGTGACGGATCGCACCAGCCTGACCAGAGATACCACCACCTTTAACGGTGATGTACAGGTCCAGTTTCTCAACCATGTCGACCAGTTCCAGCGGCTGACG  >  NC_000913/3377775‑3378063

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