Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,574,062 A→T 17.3% C273S (TGC→AGC)  asd ← aspartate‑semialdehyde dehydrogenase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,574,0620AT17.3% 449.6 / 82.3 192C273S (TGC→AGC) asdaspartate‑semialdehyde dehydrogenase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (80/77);  new base T (17/16);  total (98/94)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.06e-01

Only 100 of 194 total aligned reads displayed.

GAACGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCACTCTTCGCG  >  NC_000913/3573963‑3574161
                                                                                                   |                                                                                                   
gAACGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCa                                                                                                     >  1:2953015/1‑100 (MQ=255)
   cGACTTTCGCCCACGGATTGTGCGCGGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTTAATGCCTGGCTGTGGCAgcg                                                                                                  <  2:1408207/100‑1 (MQ=255)
    gACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAgcgc                                                                                                 <  1:2241812/100‑1 (MQ=255)
       tttCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAAt                                                                                              <  1:3712908/100‑1 (MQ=255)
         tCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTTATAGTGAATGCCTGGCTGTGGCTGCGCAATGc                                                                                            <  1:3958481/100‑1 (MQ=255)
          cGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGcc                                                                                           <  1:445240/100‑1 (MQ=255)
          cGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGcc                                                                                           <  2:2208368/100‑1 (MQ=255)
           gCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCTTGGCTGTGGCTGCGCAATGccc                                                                                          >  2:58972/1‑100 (MQ=255)
           gCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGccc                                                                                          >  2:777043/1‑100 (MQ=255)
             ccACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCg                                                                                        <  1:1773683/100‑1 (MQ=255)
             ccACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCg                                                                                        >  1:3450290/1‑100 (MQ=255)
              cACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCg                                                                                        <  2:1148032/99‑1 (MQ=255)
               aCGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGac                                                                                      >  1:154214/1‑100 (MQ=255)
                cGGATTGTGCGCAGCCAGCAGTTCTTCCGCGGGCGGAATGGACACATCTTTTGTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGaca                                                                                     <  2:2931798/100‑1 (MQ=255)
                  gATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACg                                                                                   >  1:2429194/1‑100 (MQ=255)
                    ttGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGca                                                                                 <  2:1304656/100‑1 (MQ=255)
                    ttGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGca                                                                                 <  2:2751767/100‑1 (MQ=255)
                     tgtgCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGcac                                                                                <  2:4018937/100‑1 (MQ=255)
                       tgCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGcacac                                                                              >  1:3779550/1‑100 (MQ=255)
                        gcgcAGCCAGCAGTTCTTCCACGGTCGGAATAGCCACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATTCCCCGACACGaacaca                                                                             >  2:3472121/1‑100 (MQ=255)
                        gcgcAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGcacacn                                                                             <  2:3323926/100‑2 (MQ=255)
                          gcAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATa                                                                           >  1:170095/1‑100 (MQ=255)
                            agccagcAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATaaa                                                                         <  1:2457827/100‑1 (MQ=255)
                            agccagcAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATaaa                                                                         <  1:1958133/100‑1 (MQ=255)
                             gccagcAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAAc                                                                        >  2:3087522/1‑100 (MQ=255)
                             gccagcAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAACGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAAc                                                                        >  1:1040689/1‑100 (MQ=255)
                             gacagcagTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGACTGCCTGGCTGTGGCAGCGCAATGGCCCGACACGCACACATAAAc                                                                        <  2:2853732/98‑1 (MQ=255)
                                 gcagTTCTTCCACGGTCGGAATAGACACATCTTTTTTTTTTTTATTAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATc                                                                    >  2:1111301/1‑100 (MQ=255)
                                 gcagTTCTTCCACGGTCGGAATAGACACATCTTTTTTCTCTTTAATCGTGAATGCCTGGCTGTGGCCGCGCAATGCCCCGCCCCGCACACATAAACCATc                                                                    >  2:1042296/1‑100 (MQ=255)
                                   agTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTa                                                                  >  2:3021591/1‑100 (MQ=255)
                                   agTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTa                                                                  >  1:983116/1‑100 (MQ=255)
                                    gTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTAc                                                                 <  1:748911/100‑1 (MQ=255)
                                    gTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTAc                                                                 >  1:2215883/1‑100 (MQ=255)
                                     ttcttcCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTAcc                                                                <  1:926549/100‑1 (MQ=255)
                                      tcttcCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCg                                                               >  1:2082168/1‑100 (MQ=255)
                                       cttcCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTAcc                                                                >  2:4054513/1‑98 (MQ=255)
                                       cttcCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCgg                                                              <  1:4153764/100‑1 (MQ=255)
                                          ccACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAAt                                                           <  1:2597230/100‑1 (MQ=255)
                                          ccACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGACTGTCGCACCGCCATGCCCAGCCAAGCACACATAAACCATCTACCGCAAt                                                           >  2:2480747/1‑100 (MQ=255)
                                            aCGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTa                                                         >  2:1516545/1‑100 (MQ=255)
                                             cGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTAc                                                        >  1:594191/1‑100 (MQ=255)
                                                 cGGAATCGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGCCACGCACACATAAACCATCTACCGGAATTACGGaa                                                    >  2:802248/1‑100 (MQ=255)
                                                  ggAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGa                                                     >  2:1929955/1‑98 (MQ=255)
                                                  ggAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAg                                                   >  1:3868615/1‑100 (MQ=255)
                                                  ggAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACAGAAg                                                   >  1:1438284/1‑100 (MQ=255)
                                                     aTAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGAtg                                                <  2:1664108/100‑1 (MQ=255)
                                                     aTAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGAtg                                                <  2:1224789/100‑1 (MQ=255)
                                                      tAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTTCCGTAATTCCGGTAGAtgg                                               <  1:1765401/100‑2 (MQ=255)
                                                       aGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGAtgtg                                              <  2:4113015/100‑1 (MQ=255)
                                                        gACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATgtgt                                             >  2:3646356/1‑100 (MQ=255)
                                                        gACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTTCCGTAATTCCGGTAGAtggtt                                             <  1:291041/100‑4 (MQ=38)
                                                        gACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATgtgt                                             <  2:2104930/100‑1 (MQ=255)
                                                        gACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATgtgt                                             >  1:728048/1‑100 (MQ=255)
                                                        gACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATgtgt                                             >  1:2034940/1‑100 (MQ=255)
                                                        gACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATgtgt                                             >  1:3250938/1‑100 (MQ=255)
                                                         acacATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATgtgt                                             <  1:2326226/99‑1 (MQ=255)
                                                          cacaTCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTg                                           >  1:567339/1‑100 (MQ=255)
                                                          cacaTCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTTCCGTAATTCCGGTAGAtggttta                                           <  1:3605140/100‑6 (MQ=37)
                                                            caTCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAg                                         >  2:266987/1‑100 (MQ=255)
                                                            caTCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAg                                         >  1:3175064/1‑100 (MQ=255)
                                                            caTCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTTCCGTAATTCCGGTAGAtggtttatg                                         <  2:205315/100‑8 (MQ=37)
                                                             aTCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAgg                                        >  1:1184246/1‑100 (MQ=255)
                                                             aTCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAgg                                        >  1:1227196/1‑100 (MQ=255)
                                                             aTCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAgg                                        <  1:2984753/100‑1 (MQ=255)
                                                              tCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGa                                       >  1:3035530/1‑100 (MQ=255)
                                                               cTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGAt                                      >  1:137442/1‑100 (MQ=255)
                                                                 tttttCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCt                                    <  1:2039967/100‑1 (MQ=255)
                                                                 tttttCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCt                                    <  2:3955014/100‑1 (MQ=255)
                                                                    ttCAATTTAATAGTGAATGCCTGGCTGTGGCTGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCttgt                                 <  1:266990/100‑1 (MQ=255)
                                                                     tCAATTTAATAGTGACTGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCttgtt                                >  2:1801947/1‑100 (MQ=255)
                                                                     tCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCttgtt                                >  1:1854772/1‑100 (MQ=255)
                                                                     tCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCttgtt                                <  2:3815238/100‑1 (MQ=255)
                                                                      cAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACCCATAAACAATCTACCGGAATTACGGAAAATGTGTTGAGGATCttgttg                               >  1:3754565/1‑100 (MQ=255)
                                                                      cAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCttgttg                               >  2:2085678/1‑100 (MQ=255)
                                                                      cAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCttgttg                               <  2:875858/100‑1 (MQ=255)
                                                                       aaTTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTgg                              >  2:738969/1‑100 (MQ=255)
                                                                       aaTTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTgg                              <  2:1171001/100‑1 (MQ=255)
                                                                         tttAATAGTGAATGACTGGCTGTGGCAGCGCAATGCCCCGACACGCAGACATAAATCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGtt                            <  2:3124816/100‑1 (MQ=255)
                                                                         tttAATAGGAAATGCCCGGCTGTCGGCCCGCAATGCCCCGATACGTTCACATAAACTATCTGCCGGAATTACGGAAGAGGTGTTGAGGATCttgttg                               <  2:3347472/97‑1 (MQ=255)
                                                                          ttAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGttt                           >  2:483835/1‑100 (MQ=255)
                                                                          ttAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGttt                           <  1:2927098/100‑1 (MQ=255)
                                                                           tAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTc                          >  2:1624941/1‑100 (MQ=255)
                                                                               aGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGcc                      <  2:2429179/100‑1 (MQ=255)
                                                                                 tGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTg                    <  2:23865/100‑1 (MQ=255)
                                                                                 tGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTg                    <  2:1743200/100‑1 (MQ=255)
                                                                                 tGAATGACTGGCTGTGGCTGCGCAATGCCCCGAGACGCAGACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTg                    <  2:3585585/100‑1 (MQ=255)
                                                                                   aaTGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGcc                  <  2:3412069/100‑1 (MQ=255)
                                                                                    aTGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGccc                 <  2:1065616/100‑1 (MQ=255)
                                                                                    aTGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGccc                 <  1:3455082/100‑1 (MQ=255)
                                                                                       ccTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCttt              <  1:3975867/100‑1 (MQ=255)
                                                                                       ccTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCttt              >  1:2717503/1‑100 (MQ=255)
                                                                                         tGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCTCCTGCCCTTTcc            <  1:2738163/100‑1 (MQ=255)
                                                                                         tGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTcc            <  1:3807490/100‑1 (MQ=255)
                                                                                         tGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTcc            <  1:679467/100‑1 (MQ=255)
                                                                                         tGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGAGCTTGTTGGTTTCCGCCTGCCCTTTcc            >  2:3967908/1‑100 (MQ=255)
                                                                                          ggCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCa           >  2:928298/1‑100 (MQ=255)
                                                                                           gCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCAc          <  1:2443310/100‑1 (MQ=255)
                                                                                               tgGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCActct       >  2:3249703/1‑99 (MQ=255)
                                                                                                ggCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCACTCtt      >  1:1053110/1‑99 (MQ=255)
                                                                                                ggCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCACTCTTc     >  1:1476735/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GAACGACTTTCGCCCACGGATTGTGCGCAGCCAGCAGTTCTTCCACGGTCGGAATAGACACATCTTTTTTCAATTTAATAGTGAATGCCTGGCTGTGGCAGCGCAATGCCCCGACACGCACACATAAACCATCTACCGGAATTACGGAAGATGTGTTGAGGATCTTGTTGGTTTCCGCCTGCCCTTTCCACTCTTCGCG  >  NC_000913/3573963‑3574161

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