Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,425,646 A→G 6.9% pseudogene (17/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,6460AG 6.9% 262.4 / 6.6 88pseudogene (17/453 nt)ydaYpseudogene, Rac prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base A (40/41);  new base G (2/4);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

GCTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTAT  >  NC_000913/1425498‑1425791
                                                                                                                                                    |                                                                                                                                                 
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                                                                                                                                 <  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: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:769753/1‑149 (MQ=255)
                                                                   ggcTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGa                                                                                >  1:775020/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:918330/1‑149 (MQ=255)
                                                                                        ttCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCa                                                           >  1:1574679/1‑149 (MQ=255)
                                                                                          cTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCaaa                                                         <  2:1393812/149‑1 (MQ=255)
                                                                                          cTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCaaa                                                         <  2:12117/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                                                                                                           >  2:614801/1‑90 (MQ=255)
                                                                                                   tCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTgag                                                                                                           <  1:614801/90‑1 (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:129514/149‑1 (MQ=255)
                                                                                                                       aCGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGAtt                            >  1:63802/1‑149 (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)
                                                                                                                                                    |                                                                                                                                                 
GCTGGACTAATCCGGTGCTGATTGACGAAAAGGGCGAAATTATTGCGGGTCACGGTCGTGTTATGGCGGCTGAAATGCTCAAAATGGATTCTGTTCCGGTCATTGTTCTGTCTGGCCTGACGGATGAGCAGAAGCAGCGATAACGATCAGTATCGCTCCCGTAATGCATTAATCCGTCGCCACATTGAGAAAATGGATGCCAGTTTGCACGTCGGAACGAAGGAGTTTGATATTTCAAAGGTTTCCGAGGTGGATTCTGTTGATGATTTACTCATTGATAATGCCGCTCGTTAT  >  NC_000913/1425498‑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: