Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,396,083 T→A 6.1% V3D (GTC→GAC)  insH1 → IS5 transposase and trans‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,396,0830TA 6.1% 195.4 / 3.5 66V3D (GTC→GAC) insH1IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base T (31/29);  new base A (2/2);  total (34/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 = 4.71e-01

CTCCAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTT  >  NC_000913/1395938‑1396185
                                                                                                                                                 |                                                                                                      
ctccAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCAt                                                                                                     <  6:263207/149‑1 (MQ=255)
   cAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGATCAAGTCGCGGATATGAGATCATGTTTGTCATcag                                                                                                  >  6:285967/1‑147 (MQ=255)
   cAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCt                                                                                                   <  3:80035/148‑1 (MQ=255)
    aGATAATTGTCGTCATTCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTg                                                                                                  >  2:3630/1‑148 (MQ=255)
    aGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTg                                                                                                  >  4:25489/1‑148 (MQ=255)
    aGATAATTGTCGTCAAGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAAAATTAGCTAAGGAACATGCGAACAAGTCCCTGATAAGAGATCATGTTTGTCATCTg                                                                                                  >  4:323463/1‑148 (MQ=39)
      aTAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAg                                                                                               >  3:308042/1‑149 (MQ=255)
         aTTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGTAAGGAAGTAGCAACATTAGCTAAGTAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCa                                                                                            <  3:25037/149‑1 (MQ=255)
          ttGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTAATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCa                                                                                            >  5:153420/1‑148 (MQ=255)
               gtcATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGaa                                                                                       <  6:227242/148‑1 (MQ=255)
               gtcATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGaa                                                                                       <  6:232413/148‑1 (MQ=255)
               gtcATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAAc                                                                                      <  5:191434/149‑1 (MQ=255)
                    gCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGt                                                                                 >  2:249420/1‑149 (MQ=255)
                     acGGAAACCTCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGt                                                                                 <  5:55521/147‑1 (MQ=255)
                      cGGAAACCCCTGGCTGGGCTATTTTACAGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTtc                                                                               <  1:3630/149‑1 (MQ=255)
                      cGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGtt                                                                                >  6:323720/1‑148 (MQ=255)
                      cGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGtt                                                                                >  5:102805/1‑148 (MQ=255)
                      cGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGtgc                                                                               <  1:38551/149‑3 (MQ=255)
                      cGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTtc                                                                               <  7:75789/149‑1 (MQ=255)
                       ggAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTtca                                                                              >  6:224191/1‑149 (MQ=255)
                        gAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCACCATTAGCTAAGGAAGGTGCGAACAAGTCCCTGAGATGAGATCATGTTTGTCATCTGGAGCCATAGATCAGGGTtcat                                                                             >  5:240261/1‑149 (MQ=255)
                        gAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTtcat                                                                             >  6:337838/1‑149 (MQ=255)
                        gAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTtcat                                                                             >  4:210142/1‑149 (MQ=255)
                        gAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTtcat                                                                             >  1:131295/1‑149 (MQ=255)
                         aaaCCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTtcatc                                                                            <  2:292399/149‑1 (MQ=255)
                             cccTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTcatcat                                                                          >  4:225751/1‑147 (MQ=255)
                             cccTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTcatcat                                                                          >  4:278539/1‑147 (MQ=255)
                             cccTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTcatcat                                                                          >  2:44672/1‑147 (MQ=255)
                             cccTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGa                                                                        <  1:249420/149‑1 (MQ=255)
                             cccTGGCGGGGCTATTTTACCGCGACAATACATTCAGATAATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAATGAAGGTGCGAACAAGACCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTcatcat                                                                          >  7:18337/1‑147 (MQ=255)
                               cTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGt                                                                      >  8:205703/1‑149 (MQ=255)
                                   cGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGtcatc                                                                  <  1:44672/149‑1 (MQ=255)
                                       gCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACt                                                              <  4:111986/149‑1 (MQ=255)
                                        cTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACt                                                              >  7:191806/1‑148 (MQ=255)
                                           ttttACCTCGACAATTCATTAAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTAcc                                                          >  4:127922/1‑149 (MQ=255)
                                                ccGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGc                                                     >  7:260663/1‑149 (MQ=255)
                                                   cgACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAg                                                                                         <  7:42110/110‑1 (MQ=255)
                                                   cgACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAg                                                                                         >  8:42110/1‑110 (MQ=255)
                                                   cgACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGa                                                  <  6:41051/149‑1 (MQ=255)
                                                       aaTTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGt                                              <  5:337838/149‑1 (MQ=255)
                                                         ttcattcaGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGACATCTGGAGCCATAg                                                                                         >  3:145991/1‑104 (MQ=255)
                                                         ttcattcaGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGACATCTGGAGCCATAg                                                                                         >  3:151250/1‑104 (MQ=255)
                                                         ttcattcaGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGACATCTGGAGCCATAg                                                                                         <  4:145991/104‑1 (MQ=255)
                                                         ttcattcaGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGACATCTGGAGCCATAg                                                                                         <  4:151250/104‑1 (MQ=255)
                                                            attcaGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGt                                                                                                        <  1:305571/86‑1 (MQ=255)
                                                            attcaGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGt                                                                                                        >  2:305571/1‑86 (MQ=255)
                                                            attcaGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAAtt                                         <  6:115077/149‑1 (MQ=255)
                                                             ttcaGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGtcatc                                                                  >  8:219322/1‑123 (MQ=255)
                                                             ttcaGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGtcatc                                                                  <  7:219322/123‑1 (MQ=255)
                                                                atatcatcAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCTAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTcag                                      <  6:55125/146‑1 (MQ=255)
                                                                 gATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTc                                                                                                       <  8:270126/82‑1 (MQ=255)
                                                                 gATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTc                                                                                                       >  7:270126/1‑82 (MQ=255)
                                                                      tcaATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTc                                                       <  8:263515/125‑1 (MQ=255)
                                                                      tcaATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTc                                                       >  7:263515/1‑125 (MQ=255)
                                                                        aaTAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGGCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAgc                              <  6:274640/148‑1 (MQ=255)
                                                                               aGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCa                       <  3:210142/148‑1 (MQ=255)
                                                                                 ggaaggaaGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCCTCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGAc                    >  4:307212/1‑149 (MQ=255)
                                                                                 ggaaggaaGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGGCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTACGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGAc                    >  8:116539/1‑149 (MQ=255)
                                                                                   aaggaagTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAAtt                                         <  6:355433/126‑1 (MQ=255)
                                                                                   aaggaagTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAAtt                                         >  5:355433/1‑126 (MQ=255)
                                                                                   aaggaagTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAgcc                    >  6:188306/1‑145 (MQ=255)
                                                                                   aaggaagTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCa                  <  8:191806/149‑1 (MQ=255)
                                                                                      gaagTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCATTATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCATCAGTAAGCGCCGTCAGACCAGaa               <  6:289378/149‑1 (MQ=255)
                                                                                         gTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGaa               >  8:31234/1‑146 (MQ=39)
                                                                                           aGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGCCAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGaaaa             >  3:286623/1‑146 (MQ=38)
                                                                                            gCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAagag          <  4:311429/148‑1 (MQ=37)
                                                                                              aaCATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCTACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGAtt       <  7:116539/149‑1 (MQ=25)
                                                                                                    aGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCtt  <  1:54869/148‑1 (MQ=11)
                                                                                                                                                 |                                                                                                      
CTCCAGATAATTGTCGTCATGCCGGAAACCCCTGGCGGGGCTATTTTACCGCGACAATTCATTCAGATCATCAATAGTCAGGGAAGGAAGTAGCAACATTAGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAAAGAGATTTTCTT  >  NC_000913/1395938‑1396185

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: ‑