Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,633,785 C→A intergenic (‑40/‑12) yjjX ← / → ytjC non‑canonical purine NTP phosphatase, ITPase/XTPase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,633,7850CA98.5% 196.4 / NA 67intergenic (‑40/‑12)yjjX/ytjCnon‑canonical purine NTP phosphatase, ITPase/XTPase/phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (33/33);  minor base T (1/0);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAACCCGTGC  >  NC_000913/4633648‑4633924
                                                                                                                                         |                                                                                                                                           
gCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACATACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGc                                                                                                                                  <  2:529620/149‑1 (MQ=255)
 cAATATGGCAGGATCGTTCGCAGAAGATCTCGTGATATGCATGCAGAATGGACTGAATTTTAGCGGGATTGGTGGTCGCACTGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAATAAAGCa                                                                                                                                 <  1:700569/149‑1 (MQ=255)
  aaTATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTTGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCAt                                                                                                                                >  2:156677/1‑149 (MQ=255)
  aaTATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGCCAACTTGGTGCATAATAAGCATTACTCAGAAAATTTACGTTACAGCAGTATAAGGACAAAAAGCAt                                                                                                                                >  1:835955/1‑149 (MQ=255)
  aaTATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCAt                                                                                                                                <  2:1077040/149‑1 (MQ=255)
   atatGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATg                                                                                                                               >  1:273039/1‑149 (MQ=255)
         cAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAg                                                                                                                         <  2:901596/149‑1 (MQ=255)
         cAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGAGGTCGCACAGCCAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGTAAAAAAAGCAAGTAACAg                                                                                                                         >  1:1037753/1‑149 (MQ=255)
          aGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAgg                                                                                                                        >  2:442101/1‑149 (MQ=255)
          aGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAgg                                                                                                                        <  2:560202/149‑1 (MQ=255)
                cTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATAcc                                                                                                                  <  1:1130845/149‑1 (MQ=255)
                cTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATAcc                                                                                                                  <  1:1141051/149‑1 (MQ=255)
                cTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCAAGTTACAGTTATAca                                                                                                                  >  2:23964/1‑148 (MQ=255)
                      cGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTcc                                                                                                            >  1:671840/1‑149 (MQ=255)
                      cGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTcc                                                                                                            <  1:515012/149‑1 (MQ=255)
                       gAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATATGGAAAAAAAGCATGTTACAGGTATACCTAGTCCg                                                                                                           >  1:79497/1‑149 (MQ=255)
                       gAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCg                                                                                                           >  2:479668/1‑149 (MQ=255)
                       gAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGAGTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCg                                                                                                           >  2:69674/1‑149 (MQ=255)
                              tcGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGt                                                                                                    <  2:615522/149‑1 (MQ=255)
                              tcGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGt                                                                                                    >  2:385133/1‑149 (MQ=255)
                                gTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGa                                                                                                  <  2:318297/149‑1 (MQ=255)
                                  gAATTGCCGGCAGCATTGCCTGAATTTTCGCGTGACTGGTGGTAGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATAAACGTTACAGCAGTATAAGGTAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGaaa                                                                                                <  2:403237/149‑1 (MQ=255)
                                         cTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTg                                                                                         <  2:1078088/149‑1 (MQ=255)
                                             aGAATGGCCTGAATTTTAGCGGGTTTTGTGGTCGCACAGACAACTTGGTGCATAATCATCATTACTCAGAAAATTACCGTTACAGCAGTATAAGGATAAACAGGATGTTACAGGTATCCCTAGTCCGCCACGGTGAAACGCAGTGGAAc                                                                                     <  1:218097/149‑1 (MQ=255)
                                                  ggCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGa                                                                                <  1:627034/149‑1 (MQ=255)
                                                   gCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAg                                                                               >  2:897745/1‑149 (MQ=255)
                                                     cTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCg                                                                             <  2:934082/149‑1 (MQ=255)
                                                       gAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGAc                                                                           <  2:229395/149‑1 (MQ=255)
                                                       gAATTTTAGCGGGATGGGTGGTCGCAGAGACAACTTGGTGCATAATCAGCATTTCTCAGAAAATTAACGGTACAGCAGTATAAGGAAAACAAGCATCTTACAGGTATACCTAGTCCGCCACGGGGAAACGCAGTGGAACGCCGAGCGAc                                                                           <  2:330435/149‑1 (MQ=255)
                                                       gAATGTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAACCGCAGTGGAACGCCGAGCGAc                                                                           <  2:590326/149‑1 (MQ=255)
                                                          ttttAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTa                                                                        >  1:1014002/1‑149 (MQ=255)
                                                                cGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAgg                                                                  >  2:204321/1‑149 (MQ=255)
                                                                cGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTAATAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAgg                                                                  <  2:854201/149‑1 (MQ=255)
                                                                      tggtggTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGt                                                            <  1:463566/149‑1 (MQ=255)
                                                                       ggtggtCGCACAGACAACTTGGTGCAGAATCAGCATTACTCAGAAAATTACCGTTACAGCTGTTTAAGGAAAAACAGCTTGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATACAGGGCCAGTc                                                           >  2:1011223/1‑149 (MQ=255)
                                                                        gtggtCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGGTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCAGAGCGACGTAATAAGGGCCAGTCt                                                          >  2:1033658/1‑149 (MQ=255)
                                                                         tggtCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTg                                                         >  2:1069680/1‑149 (MQ=255)
                                                                              gCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACATAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGc                                                    >  1:264364/1‑149 (MQ=255)
                                                                                acaGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCAAGTTGAAACTCAGTTGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGccc                                                  >  2:694288/1‑149 (MQ=255)
                                                                                 aaGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGGGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGGCTGACAGCCCg                                                 <  1:135066/148‑1 (MQ=255)
                                                                                  aGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGc                                                >  2:444531/1‑149 (MQ=255)
                                                                                  aGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCACTAAAAGGAAAAAAAGCTTGTTAAAGGTATAACTAGTCCGCCACGGTCAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCAAGTCTGCCAGccccc                                                >  2:305115/1‑147 (MQ=255)
                                                                                      aaCTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGAc                                            >  1:58151/1‑149 (MQ=255)
                                                                                       aCTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGAcc                                           <  1:775910/149‑1 (MQ=255)
                                                                                              gCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGCAAGCCCGCTGACCGCCAAAg                                    >  2:290454/1‑149 (MQ=255)
                                                                                              gCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAg                                    <  1:204321/149‑1 (MQ=255)
                                                                                              gCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAg                                    >  2:818134/1‑149 (MQ=255)
                                                                                               cATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGATCGACGTATTCAGGGCCACTCTGACAGCCCGCTGACCGCCAAAgg                                   >  1:320024/1‑149 (MQ=255)
                                                                                                aTAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGt                                  <  2:206380/149‑1 (MQ=255)
                                                                                                  aaTCCGCATTACTCAGAAAATTAACGTTACAGCAGTCTAAGGAAAAAACTCATGTTACAGGTATACCTAGTCCGCCACGGTGAACCGCCGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGa                                <  1:1151699/149‑1 (MQ=255)
                                                                                                  aaTCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAACCGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGa                                >  2:240003/1‑149 (MQ=255)
                                                                                                  aaTCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGa                                <  1:479668/149‑1 (MQ=255)
                                                                                                   aTCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAg                               <  2:273039/149‑1 (MQ=255)
                                                                                                     cAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGAAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCa                             >  1:33957/1‑149 (MQ=255)
                                                                                                      aGCATTACTCAGAAAATTAACGATACAGCAGTAAAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCaa                            <  1:376981/149‑1 (MQ=255)
                                                                                                       gCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGAAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAAc                           >  1:513632/1‑149 (MQ=255)
                                                                                                          ttACTCAGAAAATTAACGTTACAGCAGTATAAGGAAATAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTTAGCAACAgg                        <  1:277307/149‑1 (MQ=255)
                                                                                                            aCTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCg                      <  1:268794/149‑1 (MQ=255)
                                                                                                             cTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGa                     <  1:385133/149‑1 (MQ=255)
                                                                                                             cTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGa                     >  1:1080483/1‑149 (MQ=255)
                                                                                                                aGAAAATTAACGTTACAGCAGTCTAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGTCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCGGCTGACCGCCAAAGGTGAGCAACAGGCGATGc                  <  1:156677/149‑1 (MQ=255)
                                                                                                                aGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGc                  >  1:275389/1‑149 (MQ=255)
                                                                                                                    aaTTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGt              >  1:108747/1‑149 (MQ=255)
                                                                                                                     aTTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTg             >  1:218036/1‑149 (MQ=255)
                                                                                                                       tAACGTGACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTACGCCACGGTGAAACGCAGAGGAACGCCGAGCGAAGTATTCAGGGCCAGTCTTAAAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGTTGGc           >  1:699018/1‑149 (MQ=255)
                                                                                                                                agcagTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCCACCCGTGc  >  1:797733/1‑149 (MQ=255)
                                                                                                                                agcagTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAACCCGTGc  <  1:998058/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                           
GCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAACCCGTGC  >  NC_000913/4633648‑4633924

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: