Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,629,470 (T)5→6 coding (69/2736 nt) rbbA ← ribosome‑associated ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,629,4651.T95.4% 214.1 / 2.1 65coding (74/2736 nt)rbbAribosome‑associated ATPase
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base T (33/29);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

ATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCAC  >  NC_000913/3629330‑3629598
                                                                                                                                         |                                                                                                                                      
aTCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGtttt                                                                                                                                    <  1:1465454/140‑1 (MQ=255)
 tCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGttttt                                                                                                                                   >  1:1395723/1‑140 (MQ=255)
 tCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGttttt                                                                                                                                   >  2:227359/1‑140 (MQ=255)
   acaTTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGttttttc                                                                                                                                 <  2:1405227/140‑3 (MQ=255)
     aTTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGttttttcca                                                                                                                               <  2:861567/140‑5 (MQ=255)
      ttGCCCTGTTCAATGACGCGGGCACCGGAAATCAACAACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGttttttccat                                                                                                                              >  1:819474/1‑135 (MQ=255)
         cccTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAAt                                                                                                                           >  2:1434753/1‑140 (MQ=255)
                  aTGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCa                                                                                                                  <  1:1434753/140‑1 (MQ=255)
                   tGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCAc                                                                                                                 >  1:1224219/1‑140 (MQ=255)
                      cgcgGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCAcgcc                                                                                                              >  1:1099555/1‑140 (MQ=255)
                        cgGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCAcgcccg                                                                                                            >  2:766972/1‑140 (MQ=255)
                        cgGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCAcgcccg                                                                                                            >  1:709216/1‑140 (MQ=255)
                        cgGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCAcgcccg                                                                                                            >  1:928092/1‑140 (MQ=255)
                        cgGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCAcgcccg                                                                                                            >  1:342119/1‑140 (MQ=255)
                        cgGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCAcgcccg                                                                                                            >  1:1466484/1‑140 (MQ=255)
                        cgGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCAcgcccg                                                                                                            >  1:1431126/1‑140 (MQ=255)
                               cGGAAATCAACGACAACAAGCTGGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCACCGTTTTTTTCCATAATGCTGGCTCCCGGCCGCCAGTTg                                                                                                     <  2:1425168/140‑1 (MQ=255)
                                ggAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTgc                                                                                                    >  1:129142/1‑140 (MQ=255)
                                 gAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTgcg                                                                                                   >  2:1330622/1‑140 (MQ=255)
                                 gAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTgcg                                                                                                   >  2:1001832/1‑140 (MQ=255)
                                        aCGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAgg                                                                                            <  2:1466484/140‑1 (MQ=255)
                                         cGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGc                                                                                           >  2:963428/1‑140 (MQ=255)
                                                aaGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGcgggacgg                                                                                    <  2:112972/140‑1 (MQ=255)
                                                   cTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGaa                                                                                 >  2:786391/1‑140 (MQ=255)
                                                   cTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGaa                                                                                 <  2:910934/140‑1 (MQ=255)
                                                    tCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAAc                                                                                <  2:928092/140‑1 (MQ=255)
                                                     cGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAAcc                                                                               <  1:1001832/140‑1 (MQ=255)
                                                     cGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAAcc                                                                               <  2:129142/140‑1 (MQ=255)
                                                     cGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAAcc                                                                               <  1:1330622/140‑1 (MQ=255)
                                                              gACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGa                                                                      >  1:404008/1‑140 (MQ=255)
                                                                     tCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTcc                                                                         >  1:1475797/1‑130 (MQ=255)
                                                                     tCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTcc                                                                         <  2:1475797/130‑1 (MQ=255)
                                                                     tCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCa                                                               >  1:1369277/1‑140 (MQ=255)
                                                                      ccGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCAt                                                              <  1:227359/140‑1 (MQ=255)
                                                                           ccAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCgg                                                         <  1:1293278/140‑1 (MQ=255)
                                                                            cAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGc                                                        <  1:1001493/140‑1 (MQ=255)
                                                                               tCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAAc                                                     <  1:666997/140‑1 (MQ=255)
                                                                               tCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAAc                                                     <  2:1395723/140‑1 (MQ=255)
                                                                                cAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAAcc                                                    <  1:379621/140‑1 (MQ=255)
                                                                                 aGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCt                                                   <  2:1369277/140‑1 (MQ=255)
                                                                                 aGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGACGCGTCATTGCGGCAACCt                                                   <  1:786391/140‑1 (MQ=255)
                                                                                   cccGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTca                                                 >  1:183620/1‑140 (MQ=255)
                                                                                        ccATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACg                                            >  2:957296/1‑140 (MQ=255)
                                                                                          aTACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAg                                          <  1:957296/140‑1 (MQ=255)
                                                                                           tACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGGACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAgg                                         >  2:173298/1‑140 (MQ=255)
                                                                                               gCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtcgt                                     <  2:183620/140‑1 (MQ=255)
                                                                                                  ggCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAg                                  >  1:417197/1‑140 (MQ=255)
                                                                                                  ggCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAg                                  >  2:855333/1‑140 (MQ=255)
                                                                                                  ggCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAg                                  >  2:1114076/1‑140 (MQ=255)
                                                                                                  ggCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAg                                  >  1:6587/1‑140 (MQ=255)
                                                                                                    ccGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGc                                <  2:1224219/140‑1 (MQ=255)
                                                                                                        aaTATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACg                            >  2:64936/1‑140 (MQ=255)
                                                                                                           atCGAGAGTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGaa                         <  1:1114076/140‑1 (MQ=255)
                                                                                                               agagTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGttc                     <  2:1431126/140‑1 (MQ=255)
                                                                                                                 agTGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGttctt                   >  1:901551/1‑140 (MQ=255)
                                                                                                                  gTGATATTGTTCAGCGCAACGGTTTTTACCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGaa                         >  2:160209/1‑133 (MQ=255)
                                                                                                                  gTGATATTGTTCAGCGCAACGGTTTTTACCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGaa                         <  1:160209/133‑1 (MQ=255)
                                                                                                                   tGATATTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttca                 <  2:1099555/140‑1 (MQ=255)
                                                                                                                     atatTGTTCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttcatt               >  1:553693/1‑140 (MQ=255)
                                                                                                                           ttCAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCg         <  2:553693/140‑1 (MQ=255)
                                                                                                                             cAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCAcc           >  2:767853/1‑136 (MQ=255)
                                                                                                                             cAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCAcc           <  1:213874/136‑1 (MQ=255)
                                                                                                                             cAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCAcc           >  2:213874/1‑136 (MQ=255)
                                                                                                                             cAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCAcc           <  1:767853/136‑1 (MQ=255)
                                                                                                                             cAGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCa       >  1:55039/1‑140 (MQ=255)
                                                                                                                              aGCGCAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCAc      >  1:469447/1‑140 (MQ=255)
                                                                                                                               gcgcAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCAcc     <  2:417197/140‑1 (MQ=255)
                                                                                                                               gcgcAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCAcc     <  2:901551/140‑1 (MQ=255)
                                                                                                                               gcgcAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCAcc     <  1:766972/140‑1 (MQ=255)
                                                                                                                                cgcAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGGCGGGAAACAGTTTCAGATGCGTCATTGCTGCGACCTCTCCACGCAGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCAccc    >  2:578721/1‑140 (MQ=255)
                                                                                                                                  cAACGGTTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCAc  >  2:1076049/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                      
ATCACATTGCCCTGTTCAATGACGCGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCAC  >  NC_000913/3629330‑3629598

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