Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,469,455 A→G 7.8% T611A (ACC→GCC) ‡ mgtA → Mg(2(+)) importing P‑type ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,469,4550AG 7.8% 222.1 / 7.7 78T611A (ACC→GCC) ‡mgtAMg(2(+)) importing P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base A (38/33);  new base G (2/4);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

GTGGAATTACCGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCT  >  NC_000913/4469306‑4469595
                                                                                                                                                     |                                                                                                                                            
gTGGAATTACCGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTTCAGCGAAACAGTGCCATGAAGTGGGGCTGGACTCGCGAGAGGAGGTCAATGGAAGTGATATTGAAACGCTATCAGACGACGAACTGGCAAATCATGCACAGCGTa                                                                                                                                              >  2:242850/1‑150 (MQ=255)
       tACCGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCt                                                                                                                                       >  2:1324882/1‑150 (MQ=255)
        aCCGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTg                                                                                                                                      >  2:1220149/1‑150 (MQ=255)
          cGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGtt                                                                                                                                    >  1:1325887/1‑150 (MQ=255)
           gTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGTTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTGGTACGCTGtgc                                                                                                                                   <  2:680499/150‑3 (MQ=255)
           gTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTGGTACGCTGtgc                                                                                                                                   <  2:841428/150‑3 (MQ=255)
           gTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTGGTACGCTGtgc                                                                                                                                   <  1:14981/150‑3 (MQ=255)
           gTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGttt                                                                                                                                   <  1:1172042/150‑1 (MQ=255)
           gTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGttt                                                                                                                                   >  1:1502316/1‑150 (MQ=255)
            tAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTGGTACGCTGtgca                                                                                                                                  <  2:464693/150‑4 (MQ=255)
             aaaaaTCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGc                                                                                                                                 <  1:1324882/150‑1 (MQ=255)
                    cTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTg                                                                                                                          <  1:1021648/150‑1 (MQ=255)
                          ggCGACAGTGAGTTATTTGCAGCGAAAGTGTGCCATGAAGTTGGGCTGGACGCGTGAGAGGTTGTCATTGTTAGTGATATTGAAACGCTATCTGTCGACGAACTGCCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCg                                                                                                                    <  2:1407856/150‑1 (MQ=255)
                          ggCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCg                                                                                                                    >  1:1191587/1‑150 (MQ=255)
                            cGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGAt                                                                                                                  >  2:1112037/1‑150 (MQ=255)
                              aCAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGc                                                                                                                <  2:440878/150‑1 (MQ=255)
                               cAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGc                                                                                                                >  1:204686/1‑149 (MQ=255)
                               cAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGc                                                                                                                <  2:204686/149‑1 (MQ=255)
                                aGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCAt                                                                                                              <  2:1285796/150‑1 (MQ=255)
                                aGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCAt                                                                                                              <  1:863568/150‑1 (MQ=255)
                                 gTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATa                                                                                                             <  2:1227559/150‑1 (MQ=255)
                                 gTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGATGTTTGCCCGCCTGACGCCGATGCATa                                                                                                             <  1:1009283/150‑1 (MQ=255)
                                   gAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATaaa                                                                                                           >  2:718386/1‑150 (MQ=255)
                                     gttagttGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGa                                                                                                         >  1:462215/1‑150 (MQ=255)
                                       tagttGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAAc                                                                                                       <  2:295745/150‑1 (MQ=255)
                                           tGCAGCGAAAGTGTGCCATGAAGTGGGCCTGGACGCCGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTCGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCAt                                                                                                   >  2:1295833/1‑150 (MQ=255)
                                            gCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCAtt                                                                                                  >  2:346458/1‑150 (MQ=255)
                                            gCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCAtt                                                                                                  <  1:971128/150‑1 (MQ=255)
                                            gCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTCGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTCCCACGCTGTTTGCCCGCCTGACGCCGATGTATAAAGAACGCAtt                                                                                                  >  1:1499645/1‑150 (MQ=255)
                                             cAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGCCGCCGATGCATAAACACCGCATtg                                                                                                 >  1:1019806/1‑150 (MQ=255)
                                             cAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATtg                                                                                                 >  1:955865/1‑150 (MQ=255)
                                             cAGCGAAACTGTCACATGAACTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACTACGAACTGGCAAATCTTGCACAGCGTACAACGCTGATTCCCCGCCTGACGCCGATGCATAAACAACGCATtg                                                                                                 >  2:369828/1‑150 (MQ=255)
                                              aGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATtgt                                                                                                <  1:1112037/150‑1 (MQ=255)
                                               gCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATtgtg                                                                                               >  1:572903/1‑150 (MQ=255)
                                               gCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATtgtg                                                                                               <  1:1113526/150‑1 (MQ=255)
                                                  aaaGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACGACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGAcc                                                                                            <  2:1545505/150‑1 (MQ=255)
                                                       gtgCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACt                                                                                       >  1:163993/1‑150 (MQ=255)
                                                               aaGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGcgcg                                                                               <  1:346458/150‑1 (MQ=255)
                                                                 gTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGaa                                                                             >  1:955958/1‑150 (MQ=255)
                                                                      gCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCa                                                                        <  2:1325887/150‑1 (MQ=255)
                                                                        tGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCAtg                                                                      <  2:187701/150‑1 (MQ=255)
                                                                        tGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCAtg                                                                      >  1:86509/1‑150 (MQ=255)
                                                                         ggACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCAtgt                                                                     <  2:1502316/150‑1 (MQ=255)
                                                                           aCGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTgg                                                                   >  1:390179/1‑150 (MQ=255)
                                                                            cgcgGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGt                                                                  <  2:983354/150‑1 (MQ=255)
                                                                              cgGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTg                                                                >  1:1358821/1‑150 (MQ=255)
                                                                              cgGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTg                                                                <  1:1451243/150‑1 (MQ=255)
                                                                               gggAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTgg                                                               >  2:289261/1‑150 (MQ=255)
                                                                                ggAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGc                                                              <  2:874525/150‑1 (MQ=255)
                                                                                   gagGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCttt                                                           >  1:886979/1‑150 (MQ=255)
                                                                                     ggtggtCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGCCGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCGTTAt                                                         >  2:1226541/1‑150 (MQ=255)
                                                                                     ggtggtCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTAt                                                         <  1:718386/150‑1 (MQ=255)
                                                                                       tggtCATTCGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATgg                                                       <  1:253969/150‑1 (MQ=255)
                                                                                        ggtCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATggg                                                      >  1:697174/1‑150 (MQ=255)
                                                                                         gtCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGc                                                     <  1:50926/150‑1 (MQ=255)
                                                                                               ggTAGTGATATTGAAACGCTATCTTACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGt                                               <  1:1610980/150‑1 (MQ=255)
                                                                                               ggTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGt                                               >  1:227828/1‑150 (MQ=255)
                                                                                               ggTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGt                                               >  2:829640/1‑150 (MQ=255)
                                                                                                gTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTa                                              >  1:1581247/1‑150 (MQ=255)
                                                                                                   gTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTa                                           >  1:1296347/1‑150 (MQ=255)
                                                                                                       tatTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATgg                                                       <  1:736409/134‑1 (MQ=255)
                                                                                                       tatTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATgg                                                       >  2:736409/1‑134 (MQ=255)
                                                                                                         ttGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAatgatg                                     >  2:27843/1‑150 (MQ=255)
                                                                                                         ttGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCAGTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAatgatg                                     >  2:1216546/1‑150 (MQ=255)
                                                                                                               cGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGcgcccg                               >  2:682012/1‑150 (MQ=255)
                                                                                                                gCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGcgcccgc                              >  1:1202737/1‑150 (MQ=255)
                                                                                                                    tCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTa                          >  1:669183/1‑150 (MQ=255)
                                                                                                                    tCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCATACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTa                          <  1:1147003/150‑1 (MQ=255)
                                                                                                                      tGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAatg                                        >  1:844810/1‑134 (MQ=255)
                                                                                                                      tGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAatg                                        <  2:844810/134‑1 (MQ=255)
                                                                                                                         cgacgaACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACgcgc                     >  1:41414/1‑150 (MQ=255)
                                                                                                                         cgacgaACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACgcgc                     >  2:1156197/1‑150 (MQ=255)
                                                                                                                              aaCTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCgg                <  2:1255387/150‑1 (MQ=255)
                                                                                                                               aCTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGa               <  1:437402/150‑1 (MQ=255)
                                                                                                                                     aaaTCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCgg         <  2:924039/150‑1 (MQ=255)
                                                                                                                                        gggcaaacaGCGTGGTACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCAt      >  1:1287811/7‑150 (MQ=255)
                                                                                                                                         ggcaaacaGCGTGGTACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCAtt     >  1:1177000/6‑150 (MQ=255)
                                                                                                                                           tGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTc   >  2:148680/1‑150 (MQ=255)
                                                                                                                                            gCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCt  <  1:1220149/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
GTGGAATTACCGTAAAAATCCTTACTGGCGACAGTGAGTTAGTTGCAGCGAAAGTGTGCCATGAAGTGGGGCTGGACGCGGGAGAGGTGGTCATTGGTAGTGATATTGAAACGCTATCTGACGACGAACTGGCAAATCTTGCACAGCGTACCACGCTGTTTGCCCGCCTGACGCCGATGCATAAAGAACGCATTGTGACCTTACTGAAGCGCGAAGGGCATGTGGTTGGCTTTATGGGCGATGGTATTAATGATGCGCCCGCCTTACGCGCTGCGGATATCGGCATTTCT  >  NC_000913/4469306‑4469595

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