Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,425,641 C→T 6.6% pseudogene (12/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,6410CT 6.6% 277.4 / 6.9 90pseudogene (12/453 nt)ydaYpseudogene, Rac prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (39/45);  new base T (2/4);  total (41/49)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.78e-01

ATTCGGCTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCG  >  NC_000913/1425493‑1425783
                                                                                                                                                    |                                                                                                                                              
aTTCGGCTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGTGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAAc                                                                                                                                                >  2:1075010/1‑149 (MQ=255)
   cGGCTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATTGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGAt                                                                                                                                             <  1:533512/149‑1 (MQ=255)
   cGGCTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGAt                                                                                                                                             <  2:1253942/149‑1 (MQ=255)
    ggCTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGATCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATc                                                                                                                                            <  1:435213/149‑1 (MQ=255)
     gCTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTATGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCa                                                                                                                                           >  1:281513/1‑149 (MQ=255)
      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                                                                                                                         >  2:34690/1‑149 (MQ=255)
                       tGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGc                                                                                                                         <  1:1513551/149‑1 (MQ=255)
                       tGATTGACGAAAAGCGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGc                                                                                                                         <  2:543495/149‑1 (MQ=255)
                        gATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCa                                                                                                                        <  2:552343/149‑1 (MQ=255)
                         aTTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCAt                                                                                                                       <  2:392929/149‑1 (MQ=255)
                         aTTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCAt                                                                                                                       >  2:428985/1‑149 (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:775020/1‑149 (MQ=255)
                                                                        ggcTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGa                                                                        >  1:769753/1‑149 (MQ=255)
                                                                         gcTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGaa                                                                       <  1:1075010/149‑1 (MQ=255)
                                                                             aaaTGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGaa                                                                   <  2:1154788/149‑1 (MQ=255)
                                                                             aaaTGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGaa                                                                   >  1:804004/1‑149 (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:1574679/1‑149 (MQ=255)
                                                                                             ttCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCa                                                   >  1:918330/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)
                                                                                                                                                    |                                                                                                                                              
ATTCGGCTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCG  >  NC_000913/1425493‑1425783

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: