Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,074,945 T→G 15.6% N1313T (AAC→ACC)  putA ← fused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,074,9450TG15.6% 274.0 / 46.2 115N1313T (AAC→ACC) putAfused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (57/39);  new base G (5/13);  total (63/52)
Fisher's exact test for biased strand distribution p-value = 1.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.31e-01

Only 100 of 117 total aligned reads displayed.

AGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAA  >  NC_000913/1074846‑1075044
                                                                                                   |                                                                                                   
aGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGt                                                                                                     <  2:1118371/100‑1 (MQ=255)
  cGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTa                                                                                                   <  2:1104969/100‑1 (MQ=255)
   ggATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTAc                                                                                                  <  2:1794543/100‑1 (MQ=255)
      ttattattGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGc                                                                                               >  2:901986/1‑100 (MQ=255)
       tattatcGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGAGGTAACAGCCTCCAGCTACCTGTTTAACTTATAGTCATTAAGCTGGCGTTACCg                                                                                                >  1:840492/1‑98 (MQ=255)
          tattGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGcc                                                                                              <  1:2295271/97‑1 (MQ=255)
          tattGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCcagc                                                                                           >  1:225333/1‑100 (MQ=255)
          tattGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGGTAACGCcagc                                                                                           <  2:444044/100‑1 (MQ=255)
           attGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCcagcg                                                                                          >  2:302567/1‑100 (MQ=255)
           attGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGGTAACGCcagct                                                                                          <  1:2162591/100‑2 (MQ=255)
           attGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGGTAACGCcagct                                                                                          <  2:819810/100‑2 (MQ=255)
            ttGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGGTAACGCcagctt                                                                                         <  2:174971/100‑3 (MQ=255)
            ttGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGGTAACGCcagctt                                                                                         <  1:2027758/100‑3 (MQ=255)
             tGAGGGGATTCGACCACGTTAAAGATGCTGGAGGATGTTGTAACATCCTCGGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCTGTAACGCcagctt                                                                                         <  1:2336995/99‑3 (MQ=255)
             tGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCcagcggc                                                                                        >  2:521071/1‑100 (MQ=255)
              gAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGGTAACGCcagcttaa                                                                                       <  2:2331098/100‑5 (MQ=255)
              gAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGGTAACGCcagcttaa                                                                                       <  1:748771/100‑5 (MQ=255)
                ggtgATTAGACCACGTTAAACATGCGTGAGGTGGTTGTAACATCCTTTGGCTATCTGTTTAACCTATTGTCATTAAGCTGGCGTTACCGTcagcggcagc                                                                                     <  1:1281607/97‑1 (MQ=255)
                ggggATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGGTAACGCcagcttaatg                                                                                     <  1:1112566/100‑7 (MQ=255)
                ggggATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGGTAACGCcagcttaatg                                                                                     <  1:1806161/100‑7 (MQ=255)
                 gggATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGGTAACGCcagcttaatga                                                                                    <  2:54284/100‑8 (MQ=255)
                  ggATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCcagcggcagcgg                                                                                   <  1:1267458/100‑1 (MQ=255)
                  ggATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGGTAACGCcagcttaatgac                                                                                   <  1:1067993/100‑9 (MQ=38)
                   gATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGt                                                                                  >  1:720625/1‑100 (MQ=255)
                      tCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTAtt                                                                               <  1:1504296/100‑1 (MQ=255)
                        gACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTca                                                                             >  1:1528944/1‑100 (MQ=255)
                        gACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTca                                                                             <  2:56969/100‑1 (MQ=255)
                         aCCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTcac                                                                            <  1:2050784/100‑1 (MQ=255)
                           cACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCacac                                                                          <  2:1100678/100‑1 (MQ=255)
                              gTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCa                                                                       >  2:474723/1‑100 (MQ=255)
                              gTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACATATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCa                                                                       >  2:859416/1‑100 (MQ=255)
                                tAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAg                                                                      <  1:2365748/99‑1 (MQ=255)
                                tAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGc                                                                     <  1:1767509/100‑1 (MQ=255)
                                   aGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGaa                                                                  >  1:2116295/1‑100 (MQ=255)
                                   aGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGaa                                                                  <  2:1706477/100‑1 (MQ=255)
                                    gtagCCGGAGGATGTTACAACCTCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAAc                                                                 >  2:232112/4‑100 (MQ=38)
                                      agCCGGAGGATGTTACAACCTCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGc                                                               >  2:2564270/2‑100 (MQ=38)
                                       gCCGGAGGATGTTACAACCTCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGc                                                               >  1:334485/1‑99 (MQ=38)
                                       gCCGGAGGATGTTACAACCTCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCt                                                              >  2:1366405/1‑100 (MQ=38)
                                        ccGGAGGATGTTACAACCTCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTc                                                             >  2:474053/1‑100 (MQ=38)
                                        ccGGAGGATGTTACAACCTCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTc                                                             >  1:1793804/1‑100 (MQ=38)
                                        ccGGAGGATGTTACAACCTCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTc                                                             >  2:2567669/1‑100 (MQ=38)
                                        ccGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTc                                                             <  1:1034081/100‑1 (MQ=255)
                                          ggaggaGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGa                                                           >  2:1273816/1‑100 (MQ=255)
                                             ggaggTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGatat                                                        >  1:456865/1‑100 (MQ=255)
                                               aggTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATAc                                                      >  2:2322385/1‑100 (MQ=255)
                                                ggTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACa                                                     >  2:1236428/1‑100 (MQ=255)
                                                ggTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTAACACTCAGAGAACGCTCGATATAc                                                      >  1:1058754/1‑99 (MQ=255)
                                                  ttGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGc                                                   <  1:2414508/100‑1 (MQ=255)
                                                    gTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCg                                                 >  2:1695463/1‑100 (MQ=255)
                                                    gTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCg                                                 >  1:2035008/1‑100 (MQ=255)
                                                     tAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGt                                                >  1:2498477/1‑100 (MQ=255)
                                                       aCATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTc                                              <  2:294836/100‑1 (MQ=255)
                                                        cATCCTCCGGCTACCTGTTTGACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTcc                                             <  1:78210/100‑1 (MQ=255)
                                                        cATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTcc                                             >  2:395182/1‑100 (MQ=255)
                                                        cATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTcc                                             >  1:2265739/1‑100 (MQ=255)
                                                        cATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTcc                                             <  2:2265754/100‑1 (MQ=255)
                                                         aTCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCa                                            >  2:480147/1‑100 (MQ=255)
                                                          tcctccGGCTACCTGTTTAACCTATAGTCATTAATCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGCACGCTCGATATACAGCCGTTCCAg                                           >  1:389820/1‑100 (MQ=255)
                                                           cctccGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGa                                          >  1:669526/1‑100 (MQ=255)
                                                           cctccGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGa                                          >  2:1478965/1‑100 (MQ=255)
                                                           cctccGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGa                                          >  1:685944/1‑100 (MQ=255)
                                                           cctccGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGa                                          >  1:2480441/1‑100 (MQ=255)
                                                             tccGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAg                                        <  2:1028631/100‑1 (MQ=255)
                                                             tccGGCTACCTGTGTAACCTATAGTCATTAAGCTGGCGCTACCGCCAGCGGCAGCGGTATTCACACTCAGCAAACGCCCGATATAAAGCCGt                                                >  1:1595650/1‑92 (MQ=255)
                                                               cGGCTACCTGTTTAACCTATGGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGa                                      <  1:1195439/100‑1 (MQ=255)
                                                               cGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGa                                      >  2:2214627/1‑100 (MQ=255)
                                                                ggCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGat                                     <  1:1880744/100‑1 (MQ=255)
                                                                ggCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGat                                     >  2:1272115/1‑100 (MQ=255)
                                                                 ncTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGata                                    >  2:1266518/2‑100 (MQ=255)
                                                                 gCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGata                                    >  2:257567/1‑100 (MQ=255)
                                                                    aCCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAACGATATTg                                 <  2:1356124/100‑1 (MQ=255)
                                                                       tGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCtt                              <  2:2580145/100‑1 (MQ=255)
                                                                          ttAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCg                           >  1:1145367/1‑100 (MQ=255)
                                                                          ttAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCg                           <  2:675814/100‑1 (MQ=255)
                                                                          ttAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCg                           <  2:2234913/100‑1 (MQ=255)
                                                                           tAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGc                          <  2:383405/100‑1 (MQ=255)
                                                                           tAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGc                          <  1:987849/100‑1 (MQ=255)
                                                                             aCCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGcc                         >  1:2345199/1‑99 (MQ=255)
                                                                             aCCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCa                        >  1:1373136/1‑100 (MQ=255)
                                                                             aCCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCa                        >  1:781710/1‑100 (MQ=255)
                                                                             aCCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCa                        >  1:1735232/1‑100 (MQ=255)
                                                                                tataGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACgg                     <  2:1891405/100‑1 (MQ=255)
                                                                                 ataGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACggg                    <  2:1275302/100‑1 (MQ=255)
                                                                                 ataGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACggg                    <  2:2006450/100‑1 (MQ=255)
                                                                                  taGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGc                   >  2:1599023/1‑100 (MQ=255)
                                                                                  taGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGc                   <  1:208158/100‑1 (MQ=255)
                                                                                      gctgccgCTGGCGGTAACGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCaaa                >  2:1239244/7‑99 (MQ=255)
                                                                                       ctgccgCTGGCGGTAACGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAAc              >  2:692870/6‑100 (MQ=255)
                                                                                       aTTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAAc              >  2:1597751/1‑100 (MQ=255)
                                                                                         tAAGTTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAACGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAAccc            <  1:2099168/100‑1 (MQ=255)
                                                                                         tAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAAccc            >  2:152745/1‑100 (MQ=255)
                                                                                         tAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAAccc            <  2:2038838/100‑1 (MQ=255)
                                                                                           cgCTGGCGGTAACGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTg          >  2:575204/2‑100 (MQ=255)
                                                                                           cgCTGGCGGCAACGCCAGCGGCAGCGGTATTCACACGCAGCGCACGCTCGATATACAGCCGGTCCAGCAGGATATAGCTTTCGCCACGGCCCAAACCCTg          >  1:2156851/2‑100 (MQ=255)
                                                                                            gCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGc         <  2:558308/100‑1 (MQ=255)
                                                                                            gCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGc         <  2:1831695/100‑1 (MQ=255)
                                                                                               ggCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCAcc      >  1:1797869/1‑100 (MQ=255)
                                                                                               ggCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCAcc      >  2:708911/1‑100 (MQ=255)
                                                                                                   ttACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGaaa  <  1:2540108/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
AGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAA  >  NC_000913/1074846‑1075044

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: