Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,425,647 G→T 6.9% pseudogene (18/453 nt) ydaY → pseudogene, Rac prophage;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,425,6470GT 6.9% 262.2 / 8.2 87pseudogene (18/453 nt)ydaYpseudogene, Rac prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (39/42);  new base T (2/4);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.16e-01

CTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTAT  >  NC_000913/1425499‑1425791
                                                                                                                                                    |                                                                                                                                                
cTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAg                                                                                                                                                  >  1:73136/1‑149 (MQ=255)
   gACTAATCCGGTGCTGATTGACGATAAGGGCGAAATTATTGCGGGGCACGGTCGTGTTATGGCTTCTGTAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGTCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTAt                                                                                                                                               <  1:1202221/149‑1 (MQ=255)
    aCTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATc                                                                                                                                              <  1:381521/149‑1 (MQ=255)
     cTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCg                                                                                                                                             >  2:763326/1‑149 (MQ=255)
     cTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCg                                                                                                                                             >  2:129514/1‑149 (MQ=255)
      tAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGc                                                                                                                                            >  1:1066581/1‑149 (MQ=255)
      tAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGc                                                                                                                                            >  1:1384890/1‑149 (MQ=255)
      tAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGATGAAATGCTCAAAAAGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGCCGGATGAGCAGAAGCAGCGATAACGATCAGTATCGc                                                                                                                                            >  1:1061221/1‑149 (MQ=255)
       aaTCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATACTGATCGTTATCGCt                                                                                                                                           <  1:1408801/149‑1 (MQ=255)
        aTCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCta                                                                                                                                          >  1:1137707/1‑148 (MQ=255)
        aTCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTc                                                                                                                                          >  2:1142429/1‑149 (MQ=255)
        aTCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGAAGGATGAGCAGAAGCAGCGATAACGAACAGTATCGCTc                                                                                                                                          >  1:1324357/1‑149 (MQ=255)
         tCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATACTGATCGTTATCGCtgc                                                                                                                                         <  2:881529/149‑3 (MQ=255)
         tCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATACTGATCGTTATCGCtgc                                                                                                                                         <  1:26139/149‑3 (MQ=255)
            ggTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATACTGATCGTTATCGCtgcttc                                                                                                                                      <  2:399945/149‑6 (MQ=255)
               gCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATTGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAAt                                                                                                                                   >  2:971157/1‑149 (MQ=255)
               gCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCATAATGGATTCTGTTCCGGTCTTTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAAt                                                                                                                                   <  1:364364/149‑1 (MQ=255)
                 tGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGc                                                                                                                                 <  1:1513551/149‑1 (MQ=255)
                 tGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGc                                                                                                                                 >  2:34690/1‑149 (MQ=255)
                 tGATTGACGAAAAGCGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGc                                                                                                                                 <  2:543495/149‑1 (MQ=255)
                  gATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCa                                                                                                                                <  2:552343/149‑1 (MQ=255)
                   aTTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCAt                                                                                                                               >  2:428985/1‑149 (MQ=255)
                   aTTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCAt                                                                                                                               <  2:392929/149‑1 (MQ=255)
                      gACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGTCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTATTGCATTaa                                                                                                                            <  2:1029398/149‑1 (MQ=255)
                                   aaaTTATTGCGCGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGATCCGGGCATTGTTCTGTCTGGCCTGACGGCTGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACAt                                                                                                               <  1:234624/149‑1 (MQ=255)
                                     attattGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTg                                                                                                             <  2:911953/149‑1 (MQ=255)
                                      ttattGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTga                                                                                                            <  2:370038/149‑1 (MQ=255)
                                       tattGCGGGGCACGGGCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGCTCCGGGCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTgag                                                                                                           <  2:214916/149‑1 (MQ=255)
                                            cGGGTCACGGTCGTGGTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTACGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAAt                                                                                                      >  1:720009/1‑149 (MQ=255)
                                             gggTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAAAGGATTCTGTTCCGGTCATTGTTATGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGCTCAGTATAGCACCCGTAATGCATTAATCAGTCGCCACATTGAGAAAATg                                                                                                     >  2:1217514/1‑149 (MQ=255)
                                               gTCACGGTCGTGTTATGGCGGCAGAAAAGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAACAGAAGCAGCGATAACGATCAGAATCGCTCCCGTAATGCATTAATCCGTAGCCACATTGATAAAATgct                                                                                                   >  2:215325/1‑147 (MQ=255)
                                               gTCACGGACGTGGTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCACTGTTCTGTCTGGCCTTACGGATGAGCAGTAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGa                                                                                                   <  1:296186/149‑1 (MQ=255)
                                                gcACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGAt                                                                                                  <  2:241821/148‑1 (MQ=255)
                                                  aCGGTCGTGTTATGGCGGCTGAAATGCTAAAAATGGATTCTGTTACGGTAATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGc                                                                                                >  1:1360871/1‑149 (MQ=255)
                                                     gTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAg                                                                                             >  2:891932/1‑149 (MQ=255)
                                                       cGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGtt                                                                                           <  1:800041/149‑1 (MQ=255)
                                                           ttATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCa                                                                                       >  2:340509/1‑149 (MQ=255)
                                                           ttATGGCGGCTGAAATGCTCAAAATGGATTCAGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCa                                                                                       >  1:855116/1‑149 (MQ=255)
                                                            tATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGGAc                                                                                      <  1:34690/149‑1 (MQ=255)
                                                            tATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCAc                                                                                      <  2:73136/149‑1 (MQ=255)
                                                                gcggcTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTTTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCg                                                                                  <  1:842761/149‑1 (MQ=255)
                                                                gcggcTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCg                                                                                  <  2:644693/149‑1 (MQ=255)
                                                                  ggcTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGa                                                                                >  1:769753/1‑149 (MQ=255)
                                                                  ggcTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGa                                                                                >  1:775020/1‑149 (MQ=255)
                                                                   gcTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGaa                                                                               <  1:1075010/149‑1 (MQ=255)
                                                                       aaaTGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGaa                                                                           >  1:804004/1‑149 (MQ=255)
                                                                       aaaTGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGaa                                                                           <  2:1154788/149‑1 (MQ=255)
                                                                        aaTGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAg                                                                          <  2:1328628/149‑1 (MQ=255)
                                                                           gCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATTCAAGTTTGCACGTCGGAACGAAGGat                                                                       >  2:469989/1‑148 (MQ=255)
                                                                             tCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCACGTAATGAATTAATCCGTCGCCACATTGAGAAAATCGATGCCAGTTTGCACGTCGGAACGAAGGAGtt                                                                     >  1:566096/1‑149 (MQ=255)
                                                                                aaaTGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCTTTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGa                                                                  >  1:912735/1‑149 (MQ=255)
                                                                                aaaTGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGACACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGa                                                                  >  1:310127/1‑149 (MQ=255)
                                                                                  aTGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGata                                                                <  2:1589812/149‑1 (MQ=255)
                                                                                     gATTCTGTTCCGGTAATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTGGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATAttt                                                             <  1:215325/149‑1 (MQ=255)
                                                                                       ttCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCa                                                           >  1:918330/1‑149 (MQ=255)
                                                                                       ttCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCa                                                           >  1:1574679/1‑149 (MQ=255)
                                                                                         cTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCaaa                                                         <  2:12117/149‑1 (MQ=255)
                                                                                         cTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCaaa                                                         <  2:1393812/149‑1 (MQ=255)
                                                                                         cTGTTCCGGTCATTGTTATTTCTGGCCAGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGACCGAAGGAGTTTGATATTTCaaa                                                         <  1:779337/149‑1 (MQ=255)
                                                                                          tGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAg                                                        <  1:971157/149‑1 (MQ=255)
                                                                                               cGGTCATTGTGCTGTCTGGCCTGACGGAGGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTc                                                   <  1:249073/149‑1 (MQ=255)
                                                                                                  tCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTgag                                                                                                           <  1:614801/90‑1 (MQ=255)
                                                                                                  tCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTgag                                                                                                           >  2:614801/1‑90 (MQ=255)
                                                                                                   cATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAg                                               >  1:1218740/1‑149 (MQ=255)
                                                                                                          ctgtctgGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGAtt                                        <  1:1142429/149‑1 (MQ=255)
                                                                                                           tgtctgACCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGCACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTc                                       >  1:1401782/1‑149 (MQ=255)
                                                                                                              ctgGCCTGACGGATGATCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAATGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGt                                    <  1:1368665/149‑1 (MQ=255)
                                                                                                              ctgGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGt                                    <  2:44144/149‑1 (MQ=255)
                                                                                                                 gCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTtga                                 >  2:897843/1‑149 (MQ=255)
                                                                                                                 gCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATCGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTtga                                 <  1:1271704/149‑1 (MQ=255)
                                                                                                                     gACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTgatgat                             >  1:727071/1‑149 (MQ=255)
                                                                                                                      aCGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGAtt                            >  1:63802/1‑149 (MQ=255)
                                                                                                                      aCGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGAtt                            <  1:129514/149‑1 (MQ=255)
                                                                                                                          aTGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACt                        >  2:39512/1‑149 (MQ=255)
                                                                                                                           ttatcAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTc                       <  2:1066581/145‑1 (MQ=255)
                                                                                                                           tGAGCAGAAGCAGCGATAACGATCTGTATCGCTCGCGTTATGCATTAATCCGTCGCCACATTTAGAAAATGGATGACAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCATAGGTTTCCGAGTTGGTTTCTGTTGATGATTTACTc                       <  2:1061221/149‑1 (MQ=255)
                                                                                                                           tGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGCTGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTc                       <  2:809865/149‑1 (MQ=255)
                                                                                                                             aGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGATGTGGATTCTGTTGATGATTTACTCAt                     <  2:281513/149‑1 (MQ=255)
                                                                                                                              gCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCAtt                    >  2:524928/1‑149 (MQ=255)
                                                                                                                                aGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGa                  <  1:1191354/149‑1 (MQ=255)
                                                                                                                                 gAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGAt                 >  2:12304/1‑149 (MQ=255)
                                                                                                                                  cgggagcGATACTGATCGTTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATa                >  1:490721/5‑149 (MQ=255)
                                                                                                                                    gcagcGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAAt              <  2:432669/149‑1 (MQ=255)
                                                                                                                                     gagcGATACTGATCGTTATCGCTCCCGTACTGCATTAATCCGTCGCCAGATTGAGTAAATGGATGCCAGTTTGCAGGTCGGAACGATGGATTTTGATATTTCAAAGGTTTCCGAGTTGGATTCTGTTGATGATTTACTGATTGATAATg             >  1:248107/2‑149 (MQ=255)
                                                                                                                                     caacGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAGCGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATg             <  1:129342/146‑1 (MQ=255)
                                                                                                                                        cGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTGTCCGAGGTGGGTTCTGGGGATGTTTTACTCTTTGGTAAGGCCg          >  2:272231/1‑149 (MQ=255)
                                                                                                                                               gATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTa   >  2:1368032/1‑149 (MQ=255)
                                                                                                                                                aTCAGTATCGCTCCCGTAATGCATTAATCCGTCGACACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAACGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTCATGATTTACTCATTCATAAAGCCGCTAGTTAt  >  2:319857/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTAT  >  NC_000913/1425499‑1425791

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