Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,428,527 T→C 6.2% Q67R (CAG→CGG)  insH1 ← IS5 transposase and trans‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,428,5270TC 6.2% 193.6 / 3.1 64Q67R (CAG→CGG) insH1IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base T (30/30);  new base C (2/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.48e-01

ATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGT  >  NC_000913/1428382‑1428672
                                                                                                                                                 |                                                                                                                                                 
aTGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCTCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAAtgctgc                                                                                                                                                 <  3:28584/148‑1 (MQ=255)
aTGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCtgca                                                                                                                                                <  4:68170/149‑1 (MQ=255)
     gctCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAg                                                                                                                                           >  7:372877/1‑149 (MQ=255)
       tCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGt                                                                                                                                          >  8:83688/1‑148 (MQ=255)
         gAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGAGaa                                                                                                                                       <  6:76452/149‑1 (MQ=255)
               gTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACg                                                                                                                                   >  6:30303/1‑147 (MQ=255)
                  gCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACgcagca                                                                                                                              >  8:264469/1‑149 (MQ=255)
                  gCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACgcagca                                                                                                                              >  3:54597/1‑149 (MQ=255)
                  gCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACgcagc                                                                                                                               <  8:372877/148‑1 (MQ=255)
                    ggAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATg                                                                                                                            <  4:54597/149‑1 (MQ=255)
                             tgatgGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTTCAGGGCATCTTCCATGGCACCGTCGCTCAGGATGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGc                                                                                                                   <  8:68376/149‑1 (MQ=255)
                             tgatgGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCtgcatg                                                                                                                                              <  7:115720/122‑1 (MQ=255)
                             tgatgGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCtgcatg                                                                                                                                              >  8:115720/1‑122 (MQ=255)
                               atgGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCg                                                                                                                  >  5:162246/1‑148 (MQ=255)
                                  gtggtgCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATa                                                                                                              <  7:264469/149‑1 (MQ=255)
                                     gtgCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAtgg                                                                                                           <  3:240752/149‑1 (MQ=255)
                                      tgCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGt                                                                                                                                          <  3:238914/117‑1 (MQ=255)
                                      tgCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGt                                                                                                                                          >  4:238914/1‑117 (MQ=255)
                                        cGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGtc                                                                                                                         >  5:121278/1‑132 (MQ=255)
                                        cGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGtc                                                                                                                         <  6:121278/132‑1 (MQ=255)
                                              ggCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGc                                                                                                   >  2:96330/1‑148 (MQ=255)
                                                            cAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTgg                                                                                    <  5:30303/149‑1 (MQ=255)
                                                            cAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTgg                                                                                    <  6:162246/149‑1 (MQ=255)
                                                               ggATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGAt                                                                                 <  3:71992/149‑1 (MQ=255)
                                                                 aTAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTg                                                                                     <  5:143147/143‑1 (MQ=255)
                                                                 aTAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTg                                                                                     >  6:143147/1‑143 (MQ=255)
                                                                     tCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGTATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACg                                                                           >  3:55954/1‑149 (MQ=255)
                                                                     tCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAACACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACg                                                                           >  2:70311/1‑149 (MQ=255)
                                                                        ggCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATTGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCt                                                                        <  5:322789/149‑1 (MQ=255)
                                                                           aaaCAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTcc                                                                                                                      <  4:239270/100‑1 (MQ=255)
                                                                           aaaCAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTcc                                                                                                                      >  3:239270/1‑100 (MQ=255)
                                                                           aaaCAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGa                                                                     >  5:129972/1‑149 (MQ=255)
                                                                           aaaCAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGa                                                                     >  5:244814/1‑149 (MQ=255)
                                                                            aaCAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGAt                                                                    <  4:6715/149‑1 (MQ=255)
                                                                             aCAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCCGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCtt                                                                                      >  4:45185/1‑130 (MQ=255)
                                                                             aCAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCCGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCtt                                                                                      >  4:49902/1‑130 (MQ=255)
                                                                             aCAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCCGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCtt                                                                                      <  3:49902/130‑1 (MQ=255)
                                                                             aCAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCCGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCtt                                                                                      <  3:45185/130‑1 (MQ=255)
                                                                                             gATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCt                                                   <  6:129972/149‑1 (MQ=255)
                                                                                              aTTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTg                                                  >  4:80696/1‑149 (MQ=255)
                                                                                               tttCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGTAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTggg                                                                                   <  4:242659/115‑1 (MQ=255)
                                                                                               tttCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGTAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTggg                                                                                   >  3:242659/1‑115 (MQ=255)
                                                                                                   gTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATg                                             >  4:15326/1‑149 (MQ=255)
                                                                                                   gTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATg                                             <  4:91236/149‑1 (MQ=255)
                                                                                                    tACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATgg                                            >  5:99093/1‑149 (MQ=255)
                                                                                                      cAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATgg                                            >  1:53115/1‑147 (MQ=255)
                                                                                                       aGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGc                                                              <  3:323079/128‑1 (MQ=255)
                                                                                                       aGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGc                                                              >  4:323079/1‑128 (MQ=255)
                                                                                                         ggCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTg                                                                                     >  4:12771/1‑103 (MQ=255)
                                                                                                         ggCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTg                                                                                     <  3:12771/103‑1 (MQ=255)
                                                                                                         ggCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGaa                                       >  7:156723/1‑149 (MQ=255)
                                                                                                          gCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAACACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAAt                                      <  1:70311/149‑1 (MQ=255)
                                                                                                                 ccATGGCACCGTCGCTCAGGTTGTACCAATGCTGTATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATATga                                  <  4:55954/146‑2 (MQ=255)
                                                                                                                  cATGGCTCGGTCGCTGAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGAATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCGGCCTTTGGATAAAACGGCTCGATGACTGCGGTCATATTCTGCCATGGCAGAATCTGCTcc                               <  4:85302/148‑1 (MQ=18)
                                                                                                                        aCCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACa                                                                >  7:107389/1‑109 (MQ=255)
                                                                                                                        aCCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACa                                                                <  8:107389/109‑1 (MQ=255)
                                                                                                                        aCCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATg                                             <  3:337092/128‑1 (MQ=255)
                                                                                                                        aCCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATg                                             >  4:337092/1‑128 (MQ=255)
                                                                                                                          cGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCgg                         >  4:219344/1‑146 (MQ=255)
                                                                                                                           gTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGgaga                     <  6:99093/149‑1 (MQ=255)
                                                                                                                           gTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCACCATGCGGgaga                     <  3:15326/149‑1 (MQ=255)
                                                                                                                                tCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATg                                             >  2:260091/1‑120 (MQ=255)
                                                                                                                                tCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATg                                             <  1:260091/120‑1 (MQ=255)
                                                                                                                                tCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGaa                 >  5:231497/1‑148 (MQ=255)
                                                                                                                                        tACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCtt         <  2:53115/148‑1 (MQ=255)
                                                                                                                                          ccAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCttt        >  7:375512/1‑147 (MQ=255)
                                                                                                                                          ccAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTc      <  1:96330/149‑1 (MQ=255)
                                                                                                                                           cAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCg     >  4:154324/1‑149 (MQ=255)
                                                                                                                                              tgctgcATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGt  <  3:219344/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
ATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCATCTTCCATGGCACCGTCGCTCAGGTTGTACCAATGCTGCATGCAGTGAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGT  >  NC_000913/1428382‑1428672

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