Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,323,313 G→C 7.7% intergenic (‑83/+23) phnE ← / ← phnD defective phosphonate ABC transporter permease;transport; Transport of small molecules: Anions; membrane channel protein component of Pn transporter; membrane component of ABC superfamily/phosphonate ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,323,3130GC 7.7% 492.1 / 27.2 142intergenic (‑83/+23)phnE/phnDdefective phosphonate ABC transporter permease;transport; Transport of small molecules: Anions; membrane channel protein component of Pn transporter; membrane component of ABC superfamily/phosphonate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (64/67);  new base C (6/5);  total (70/72)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.50e-01

Only 100 of 144 total aligned reads displayed.

AACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTT  >  NC_000913/4323168‑4323458
                                                                                                                                                 |                                                                                                                                                 
aaCGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGc                                                                                                                                                >  4:572164/1‑149 (MQ=255)
   ggCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTc                                                                                                                                             <  3:951826/149‑1 (MQ=255)
       cAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTAt                                                                                                                                         <  3:1492395/149‑1 (MQ=255)
            cAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGAc                                                                                                                                    >  1:1326981/1‑149 (MQ=255)
            cAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGAc                                                                                                                                    <  3:572164/149‑1 (MQ=255)
             aGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGAcc                                                                                                                                   <  3:1491819/149‑1 (MQ=255)
                 gCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCggg                                                                                                                                                  >  3:518174/1‑130 (MQ=255)
                 gCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATc                                                                                                                                        <  4:1648249/140‑1 (MQ=255)
                 gCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATc                                                                                                                                        >  3:1648249/1‑140 (MQ=255)
                  cGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGa                                                                                                                              >  1:1441626/1‑149 (MQ=255)
                   gagaACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATc                                                                                                                                        >  1:1031612/1‑138 (MQ=255)
                   gagaACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATc                                                                                                                                        <  2:1031612/138‑1 (MQ=255)
                   gagaACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGAt                                                                                                                             >  4:602186/1‑149 (MQ=255)
                    agaACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGAtt                                                                                                                            >  4:776079/1‑149 (MQ=255)
                       aCCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATcc                                                                                                                                       >  2:936088/1‑135 (MQ=255)
                       aCCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATcc                                                                                                                                       <  1:936088/135‑1 (MQ=255)
                        ccAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATAAGACGCCCCGGCGTCGCATCCGa                                                                                                                                     <  4:1105040/136‑1 (MQ=255)
                        ccAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATAAGACGCCCCGGCGTCGCATCCGa                                                                                                                                     >  3:1105040/1‑136 (MQ=255)
                         cAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTAc                                                                                                                         <  3:1265635/147‑1 (MQ=255)
                                 cTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAGGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCtt                                                                                                               >  2:1309362/1‑149 (MQ=255)
                                  ttGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCCCGGCGTcg                                                                                                                                            >  2:414740/1‑118 (MQ=255)
                                           gCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGa                                                                                                     >  1:1805006/1‑149 (MQ=255)
                                             gatggtgatggtTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGAACCGCTTTACTCACCGAAc                                                                                                   <  2:1326981/149‑1 (MQ=255)
                                                  tgatggtTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATc                                                                                              <  3:602186/149‑1 (MQ=255)
                                                    atggtTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCAcc                                                                                                       >  4:936032/1‑138 (MQ=255)
                                                         ttGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTa                                                                                       <  2:1454083/149‑1 (MQ=255)
                                                               gTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCgt                                                                                 <  1:1275402/149‑1 (MQ=255)
                                                                 ttGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCa                                                                                                                     <  1:360538/111‑1 (MQ=255)
                                                                 ttGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCa                                                                                                                     >  2:360538/1‑111 (MQ=255)
                                                                  tGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCAACGCGCTAAGCGCgatga                                                                              <  1:726931/149‑2 (MQ=255)
                                                                  cttatCCGGCCTACAGACCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCt                                                                                         >  1:1269685/5‑138 (MQ=255)
                                                                    gCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTc                                                                            >  1:614461/1‑149 (MQ=255)
                                                                      tCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCgtt                                                                                <  1:1527828/143‑1 (MQ=255)
                                                                      tCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCgtt                                                                                >  2:1527828/1‑143 (MQ=255)
                                                                      tCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAg                                                                          >  4:1692075/1‑149 (MQ=255)
                                                                      tCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAg                                                                          <  3:1692075/149‑1 (MQ=255)
                                                                      tCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGAGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAg                                                                          >  2:449016/1‑149 (MQ=255)
                                                                       ccGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAgg                                                                         <  2:614461/149‑1 (MQ=255)
                                                                         ggTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCAACGCGCTAAGCGCgatga                                                                              >  1:1043983/1‑141 (MQ=255)
                                                                         ggTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCAACGCGCTAAGCGCgatga                                                                              <  2:1043983/142‑2 (MQ=255)
                                                                          gTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGAtt                                                                                                                            <  4:1805545/95‑1 (MQ=255)
                                                                          gTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGAtt                                                                                                                            >  3:1805545/1‑95 (MQ=255)
                                                                          gTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAAc                                                                                                   <  2:238369/120‑1 (MQ=255)
                                                                          gTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAAc                                                                                                   >  1:238369/1‑120 (MQ=255)
                                                                           tCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCa                                                                           >  1:1521871/1‑143 (MQ=255)
                                                                           tCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCa                                                                           <  2:1521871/143‑1 (MQ=255)
                                                                               tAGTCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAg                                                                 <  1:449016/149‑1 (MQ=255)
                                                                               tAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAg                                                                 >  2:1811053/1‑149 (MQ=255)
                                                                                  gCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTa                                                                                                             >  3:1244691/1‑102 (MQ=255)
                                                                                  gCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTc                                                              >  2:3320/1‑149 (MQ=255)
                                                                                  gCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTc                                                              <  3:776079/149‑1 (MQ=255)
                                                                                    cGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCAt                                                            >  1:1601497/1‑149 (MQ=255)
                                                                                     ggATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATc                                                           >  1:908225/1‑149 (MQ=255)
                                                                                       aTAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACGGAACTCATCGCGCTTAGCGCGTTGTTCAGGCg                                                                       >  1:553926/1‑135 (MQ=255)
                                                                                         aaGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAgg                                                                         >  1:981007/1‑131 (MQ=255)
                                                                                         aaGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAgg                                                                         <  2:981007/131‑1 (MQ=255)
                                                                                         aaGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGc                                                       <  2:282215/149‑1 (MQ=255)
                                                                                           ggTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTa                                                                                                                          >  3:141164/1‑80 (MQ=255)
                                                                                           ggTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTa                                                                                       <  2:988328/115‑1 (MQ=255)
                                                                                            gTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAg                                                                 <  3:1253991/136‑1 (MQ=255)
                                                                                            gTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAg                                                                 >  4:1253991/1‑136 (MQ=255)
                                                                                            gTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTgc                                                    >  3:491182/1‑149 (MQ=255)
                                                                                                ttGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGt                                                                                                                                              >  1:247316/1‑55 (MQ=255)
                                                                                                    aCCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATcc                                                          <  3:1805696/135‑1 (MQ=255)
                                                                                                    aCCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATcc                                                          >  4:1805696/1‑135 (MQ=255)
                                                                                                     ccGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCa                                                             >  4:828788/1‑131 (MQ=255)
                                                                                                     ccGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCa                                                             <  3:828788/131‑1 (MQ=255)
                                                                                                         aTCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCAt                                                            >  2:1667792/1‑128 (MQ=255)
                                                                                                         aTCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCgg                                       <  4:491182/149‑1 (MQ=255)
                                                                                                           ccGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCCGAATCTCcgtg                                     <  1:1619936/149‑1 (MQ=255)
                                                                                                             ggCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAgg                                                                         <  4:1356801/111‑1 (MQ=255)
                                                                                                             ggCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCggtggt                                   >  4:1026575/1‑149 (MQ=255)
                                                                                                              gCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATc                                                           <  4:565536/124‑1 (MQ=255)
                                                                                                                aTCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGtt                                  >  2:1577001/1‑147 (MQ=255)
                                                                                                                 tCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGttttt                               >  3:590195/1‑149 (MQ=255)
                                                                                                                 tCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGttttt                               >  1:63258/1‑149 (MQ=255)
                                                                                                                    aCGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCg                                        <  4:1646177/137‑1 (MQ=255)
                                                                                                                        cTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAgg                                                                         <  1:404326/100‑1 (MQ=255)
                                                                                                                          gCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAgg                                                                >  3:575711/1‑107 (MQ=255)
                                                                                                                          gCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAgg                                                                <  4:575711/107‑1 (MQ=255)
                                                                                                                           cacaATGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTTCTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCgg                                       <  4:1252030/131‑1 (MQ=255)
                                                                                                                            acaTTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTc                    <  1:3320/149‑1 (MQ=255)
                                                                                                                                tGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGc                <  2:1601497/149‑1 (MQ=255)
                                                                                                                                 gTCGGATAAGACGCCCCGGCGTCGCATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTgcgc                                                  >  2:614426/1‑114 (MQ=38)
                                                                                                                                 gTCGGATAAGACGCCCCGGCGTCGCATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCggtg                                     <  2:709241/127‑1 (MQ=38)
                                                                                                                                 gTCGGATAAGACGCCCCGGCGTCGCATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCggtg                                     >  1:709241/1‑127 (MQ=38)
                                                                                                                                    ggATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTc                    <  1:1798348/141‑1 (MQ=255)
                                                                                                                                    ggATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTc                    >  2:1798348/1‑141 (MQ=255)
                                                                                                                                     gATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCAtt           <  4:590195/149‑1 (MQ=255)
                                                                                                                                      ataagACGCCCCGGCGTCGCATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCt                                                      <  4:1887060/101‑1 (MQ=38)
                                                                                                                                      ataagACGCCCCGGCGTCGCATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCt                                                      >  3:1887060/5‑105 (MQ=38)
                                                                                                                                      ataagACGCCCCGGCGTCGCATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTc          >  1:1102708/5‑149 (MQ=255)
                                                                                                                                          gACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGtt                                  >  4:324789/1‑121 (MQ=255)
                                                                                                                                          gACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGtt                                  <  3:324789/121‑1 (MQ=255)
                                                                                                                                          gACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTc                    <  2:820794/135‑1 (MQ=255)
                                                                                                                                          gACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTc      >  4:66969/1‑149 (MQ=255)
                                                                                                                                           aCGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTcc     >  4:1400608/1‑149 (MQ=255)
                                                                                                                                              ccGGGGCGTCTTATCCGACCTACGATTTACTGCAACGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGATGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCttt  <  1:1535725/149‑1 (MQ=255)
                                                                                                                                               cGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTc                    >  3:1858136/1‑130 (MQ=255)
                                                                                                                                               cGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTc                    <  4:1858136/130‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
AACGGCCCAGCTCAGAAGCGAGAACCAGCTGCGCTTGGGTGGGGCGATGGTGATGGTTTGCATGTTTGGCTCCGGTCTGTAGGCCGGATAAGGTGCTTGCACCGCATCCGGCATCAACGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTT  >  NC_000913/4323168‑4323458

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