Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,652,242 T→G 41.8% pseudogene (692/1158 nt) intQ → pseudogene, Qin prophage; phage integrase family;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,652,2420TG41.8% 52.4 / 175.0 136pseudogene (692/1158 nt)intQpseudogene, Qin prophage; phage integrase family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base T (39/40);  new base G (27/30);  total (66/70)
Fisher's exact test for biased strand distribution p-value = 8.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.70e-01

Only 100 of 144 total aligned reads displayed.

GGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAG  >  NC_000913/1652144‑1652337
                                                                                                  |                                                                                               
   ggaCATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGgcatt                                                                                             >  1:543439/1‑100 (MQ=255)
     aCATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtcattca                                                                                           <  2:2509354/100‑1 (MQ=255)
      cATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtcattcat                                                                                          >  2:2016947/1‑100 (MQ=255)
      cATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGtcattcat                                                                                          <  2:3941328/100‑1 (MQ=255)
      cATAGATCTGAAAAATGGAACAATGAGGATCAGGAGACATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGAAAGCACTCGAGtcattcat                                                                                          >  1:3286091/1‑100 (MQ=255)
       aTAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATc                                                                                         >  2:2136414/1‑100 (MQ=255)
       aTAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATc                                                                                         >  2:1518994/1‑100 (MQ=255)
        tAGATCTGGAAAATGGAACAATGCTGGACTGGAGAAATTTAGCAAAAGACCGTTTCACGGTACTAAAAAGATAGGCCGTAACCAATCAGTtcatcc                                                                                            >  1:71118/1‑94 (MQ=255)
         aGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCtt                                                                                       >  2:2759870/1‑100 (MQ=255)
         aGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCtt                                                                                       <  2:61665/100‑1 (MQ=255)
          tccgatctaaaTGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCTGGGCATTCATCtta                                                                                      <  1:4019138/92‑1 (MQ=255)
          gATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCtta                                                                                      >  2:502504/1‑100 (MQ=255)
          gATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCtta                                                                                      >  1:3092511/1‑100 (MQ=255)
           aTCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCttat                                                                                     >  1:985952/1‑100 (MQ=255)
            tCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCttatt                                                                                    <  2:853287/100‑1 (MQ=255)
            tCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCttatt                                                                                    <  2:3214892/100‑1 (MQ=255)
             cTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCttatta                                                                                   <  2:2783141/100‑1 (MQ=255)
             cTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCttatta                                                                                   >  2:3093436/1‑100 (MQ=255)
             cTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCttatta                                                                                   >  2:3093402/1‑100 (MQ=255)
               gAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAg                                                                                 >  2:4180043/1‑100 (MQ=255)
               gAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAg                                                                                 >  2:2475912/1‑100 (MQ=255)
               gAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAg                                                                                 <  1:973050/100‑1 (MQ=255)
               gAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAg                                                                                 >  2:4146229/1‑100 (MQ=255)
               gAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAg                                                                                 >  1:3315658/1‑100 (MQ=255)
               gAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAg                                                                                 >  1:102502/1‑100 (MQ=255)
                 aaaaTGGAACAATGATGCTTAGGAGAAATTTAGCAAAAGACCGAGTCCCGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTaa                                                                                  >  1:3606167/1‑97 (MQ=255)
                 aaaaTGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGcc                                                                               <  1:1986334/100‑1 (MQ=255)
                 aaaaTGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGcc                                                                               <  2:687845/100‑1 (MQ=255)
                  aaaTGGAACATTGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCa                                                                              >  1:1487619/1‑100 (MQ=255)
                  aaaTGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCTGGGCATTCATCttaagatcgga                                                                              >  2:4019158/1‑92 (MQ=255)
                    aTGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGc                                                                            <  2:1267645/100‑1 (MQ=255)
                     tGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCa                                                                           <  1:3099056/100‑1 (MQ=255)
                       gAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAAt                                                                         >  1:2742478/1‑100 (MQ=255)
                         aCAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCg                                                                       <  2:2923286/100‑1 (MQ=255)
                             tgatgaTCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGc                                                                   <  2:1865142/100‑1 (MQ=255)
                             tgatgaTCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGc                                                                   <  1:1141656/100‑1 (MQ=255)
                             tgatgaTCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGc                                                                   <  1:2206706/100‑1 (MQ=255)
                                     aGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCgg                                                            >  1:1313195/1‑99 (MQ=255)
                                      ggAGACATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGa                                                           >  1:1271969/1‑99 (MQ=255)
                                      ggAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTTCGGAg                                                          >  2:1558306/1‑100 (MQ=255)
                                      ggAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAg                                                          >  1:4185971/1‑100 (MQ=255)
                                      ggAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAACAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTATCCgaa                                                           >  1:1061660/1‑97 (MQ=255)
                                       gagaAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAg                                                          >  1:3171711/1‑99 (MQ=255)
                                       gagaAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGt                                                         <  1:1866116/100‑1 (MQ=255)
                                       gagaAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGt                                                         <  1:1587302/100‑1 (MQ=255)
                                        agaAATTTAGCAAAAGACCGTTTCACGGTACCAAAGACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTc                                                        <  1:1523358/100‑1 (MQ=255)
                                          aaaTTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAg                                                      <  2:3278111/100‑1 (MQ=255)
                                          aaaTTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCan                                                      <  1:872036/100‑2 (MQ=255)
                                               tAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGAc                                                 >  2:2489896/1‑100 (MQ=255)
                                               tAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGAc                                                 <  1:288105/100‑1 (MQ=255)
                                                 gCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACAt                                               <  1:4241200/100‑1 (MQ=255)
                                                 gCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCaccaaattt                                               >  1:1514310/1‑92 (MQ=255)
                                                  cAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACAtt                                              >  1:3068465/1‑100 (MQ=255)
                                                   aaaaGACCGTTTCACGGTACCAAAAACACCGGCGGGAACCAATCGGGGCATTCATCTTCTTAAGCCAGCAATCCACGCTCTCCGGAGTCAGATGACATTa                                             >  2:657243/1‑100 (MQ=255)
                                                   aaaaGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTa                                             <  1:535032/100‑1 (MQ=255)
                                                     aaGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAAc                                           <  2:1124574/100‑1 (MQ=255)
                                                      aGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACg                                          <  1:1087265/100‑1 (MQ=255)
                                                         ccGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACgaga                                       >  2:3645726/1‑100 (MQ=255)
                                                         ccGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACgaga                                       >  1:2512190/1‑100 (MQ=255)
                                                         ccGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACgaga                                       >  1:2570740/1‑100 (MQ=255)
                                                           gTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACt                                     >  2:2716769/1‑100 (MQ=255)
                                                              tCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAg                                  >  2:1266890/1‑100 (MQ=255)
                                                                aCGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGc                                 <  1:2444816/99‑1 (MQ=255)
                                                                aCGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCa                                >  1:3254502/1‑100 (MQ=255)
                                                                  ggTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCaaa                              <  2:2569644/100‑1 (MQ=255)
                                                                    tACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAag                             <  1:197163/99‑1 (MQ=255)
                                                                     aCCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTCAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAaga                            <  1:1906325/99‑1 (MQ=255)
                                                                      ccAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGc                          <  2:4223070/100‑1 (MQ=255)
                                                                      ccAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGc                          <  2:2966032/100‑1 (MQ=255)
                                                                      accaaaaCACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGc                          >  2:2337597/4‑100 (MQ=255)
                                                                       cAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCa                         >  1:682440/1‑100 (MQ=255)
                                                                       cAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATAACATTAACGAGACTGAGCAAAGAGCa                         >  1:680382/1‑100 (MQ=255)
                                                                        aaaaaCACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAagagatc                        >  2:494263/1‑97 (MQ=255)
                                                                        aaaaaCACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCat                        >  1:3663367/1‑100 (MQ=255)
                                                                          aaaCACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCatat                      <  2:3348091/100‑1 (MQ=255)
                                                                           aaCACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATc                     <  2:1129865/100‑1 (MQ=255)
                                                                           aaCACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATc                     >  1:2803866/1‑100 (MQ=255)
                                                                             cacaGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCAt                   <  1:720512/100‑1 (MQ=255)
                                                                             cacaGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCAt                   <  1:2856413/100‑1 (MQ=255)
                                                                             cacaGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCAt                   <  2:1378957/100‑1 (MQ=255)
                                                                             cacaGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCAt                   <  2:2291145/100‑1 (MQ=255)
                                                                              acaGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCAtt                  <  1:4247235/100‑1 (MQ=255)
                                                                               caGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTg                 <  1:2390500/100‑1 (MQ=255)
                                                                               caGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTg                 >  1:3068860/1‑100 (MQ=255)
                                                                                aGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGa                <  2:2005980/100‑1 (MQ=255)
                                                                                aGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGa                <  2:2108878/100‑1 (MQ=255)
                                                                                  gCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATg              >  2:1637799/1‑100 (MQ=255)
                                                                                  gCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATg              <  1:1515979/100‑1 (MQ=255)
                                                                                  gCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAACATATCATTGATg              <  1:4146207/100‑1 (MQ=255)
                                                                                    gggAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGtt            >  2:1380037/1‑100 (MQ=255)
                                                                                     ggAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTc           <  2:887473/100‑1 (MQ=255)
                                                                                      gAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCa          <  2:543441/100‑1 (MQ=255)
                                                                                         ccAGTCGGGGCATTCGGATTATTAGGCCAGTAATCGACGCTATCCGGAGTAAGATGACATTAACGAGAATGAGCAAAGAGCATATCATTGATGTTCan         <  1:1067920/98‑2 (MQ=255)
                                                                                         ccAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACtt       <  1:2663971/100‑1 (MQ=255)
                                                                                           aaTCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCa     <  2:3396588/100‑1 (MQ=255)
                                                                                           aaTCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCa     <  1:2522081/100‑1 (MQ=255)
                                                                                            aTCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCag    <  1:3168555/100‑1 (MQ=255)
                                                                                            aTCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCa     <  1:2296238/99‑1 (MQ=255)
                                                                                             tCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCaga   <  2:2081602/100‑1 (MQ=255)
                                                                                              cGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCagag  <  2:1315401/100‑1 (MQ=255)
                                                                                                  |                                                                                               
GGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAG  >  NC_000913/1652144‑1652337

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