Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,678,347 T→C V256A (GTG→GCG)  b1605 → putative arginine/ornithine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,678,3470TC94.4% 244.0 / 3.9 90V256A (GTG→GCG) b1605putative arginine/ornithine antiporter
Reads supporting (aligned to +/- strand):  ref base T (4/1);  new base C (48/37);  total (52/38)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

ATAAACGTGATGTTGGCAAAGCGACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGTGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAATGATGGGACCATGGGGCGAAATCATCAT  >  NC_000913/1678257‑1678434
                                                                                          |                                                                                       
ataaACGTGATGTTGGCAAAGCGACACTGCTGGCGGTTCTCTCCGCCCTGGGCGTTTACTTGCTGGGAACGCTGCTTTCGCTGGGTGtggcgg                                                                                       >  2:3796539/1‑93 (MQ=255)
 taaACGTGATGTTGGCAAAGCGACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTggcggc                                                                                      >  1:2819589/1‑93 (MQ=255)
 taaACGTGATGTTGGCAAAGCGACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTggcggc                                                                                      >  1:1526023/1‑93 (MQ=255)
 taaACGTGATGTTGGCAAAGCGACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGTGGc                                                                                      >  1:281388/1‑93 (MQ=255)
     cGTGATGTTGGCAAAGCGACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGt                                                                                  >  2:4286618/1‑93 (MQ=255)
     cGTGATGTTGGCAAAGCGACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGt                                                                                  >  2:3193416/1‑93 (MQ=255)
        gATGTTGGCAAAGCGACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCt                                                                               >  1:3821267/1‑93 (MQ=255)
          tGTTGGCAAAGCGACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGa                                                                             >  1:2138075/1‑93 (MQ=255)
          tGTTGGCAAAGCGACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGa                                                                             >  2:987538/1‑93 (MQ=255)
           ggtGGCAAAGGGACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGaa                                                                            <  1:4006066/91‑1 (MQ=255)
            ttGGCAAAGCGACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAAc                                                                           >  2:1998781/1‑93 (MQ=255)
             tGGCAAAGCGACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACt                                                                          <  1:2257155/93‑1 (MQ=255)
              ggCAAAGCGACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTg                                                                         >  2:899443/1‑93 (MQ=255)
                cAAAGCGACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGc                                                                       >  1:3549921/1‑93 (MQ=255)
                 aaaGCGACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCt                                                                      >  2:3366949/1‑93 (MQ=255)
                 aaaGCGACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCt                                                                      <  2:4226228/93‑1 (MQ=255)
                 aaaGCGACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCt                                                                      >  2:103862/1‑93 (MQ=255)
                 aaaGCGACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCt                                                                      >  2:3511427/1‑93 (MQ=255)
                  aaGCGACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTg                                                                     >  2:527951/1‑93 (MQ=255)
                  aaGCGACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTg                                                                     >  1:4677701/1‑93 (MQ=255)
                   aGCGACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGTGGCACGTCCTGAACTGGCTGa                                                                    >  2:3975851/1‑93 (MQ=255)
                     cGACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGaaa                                                                  <  1:2179127/93‑1 (MQ=255)
                     cGACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGaaa                                                                  >  2:1857935/1‑93 (MQ=255)
                      gACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAAt                                                                 >  2:4067483/1‑93 (MQ=255)
                       acacTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGTGGCACGTCCTGAACTGGCTGAAAtt                                                                >  1:4031782/1‑93 (MQ=255)
                       acacTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGTACGTCCTGAACTGGCTGAAAtt                                                                >  1:2764552/1‑93 (MQ=255)
                       acacTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAAtt                                                                >  2:2014949/1‑93 (MQ=255)
                         acTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGTGGCACGTCCTGAACTGGCTGAAATTCg                                                              >  2:2349271/1‑93 (MQ=255)
                           tgctgGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTa                                                            >  2:4145187/1‑93 (MQ=255)
                            gctgGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTaa                                                           >  1:264214/1‑93 (MQ=255)
                             ctgGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAAc                                                          >  2:3562047/1‑93 (MQ=255)
                              tgGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAAcc                                                         <  2:106384/93‑1 (MQ=255)
                               ggCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAAccc                                                        <  2:1246510/93‑1 (MQ=255)
                                  ggTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTc                                                     <  1:2546283/93‑1 (MQ=255)
                                   gTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGTGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCa                                                    >  2:3232860/1‑93 (MQ=255)
                                   gTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCa                                                    >  1:4485058/1‑93 (MQ=255)
                                    ttCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCaa                                                   >  2:3935808/1‑93 (MQ=255)
                                    ttCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCaa                                                   <  2:1175383/93‑1 (MQ=255)
                                       tctcCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATgg                                                >  2:4072909/1‑93 (MQ=255)
                                       tctcCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATgg                                                >  1:1003795/1‑93 (MQ=255)
                                        ctcCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGc                                               <  2:3530319/93‑1 (MQ=255)
                                        ctcCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGc                                               >  1:2289514/1‑93 (MQ=255)
                                        ctcCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGc                                               <  1:2062456/93‑1 (MQ=255)
                                         tcCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGcc                                              >  1:524207/1‑93 (MQ=255)
                                         tcCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGcc                                              <  2:450043/93‑1 (MQ=255)
                                         tcCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGcc                                              >  1:1183725/1‑93 (MQ=255)
                                         tcCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGcc                                              <  1:4276155/93‑1 (MQ=255)
                                         tcCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGcc                                              <  2:3774494/93‑1 (MQ=255)
                                           cGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCgg                                            >  1:2300182/1‑93 (MQ=255)
                                           cGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCgg                                            >  2:584890/1‑93 (MQ=255)
                                           cGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCgg                                            <  1:195797/93‑1 (MQ=255)
                                            gCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGt                                           <  2:3772925/93‑1 (MQ=255)
                                             ctctGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTc                                          >  2:1940841/1‑93 (MQ=255)
                                              tctGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCt                                         >  1:262766/1‑93 (MQ=255)
                                                tGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGa                                       >  2:720733/1‑93 (MQ=255)
                                                 gggCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGAt                                      >  1:3327520/1‑93 (MQ=255)
                                                  ggCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGAtg                                     >  1:435246/1‑93 (MQ=255)
                                                  ggCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGAtg                                     <  2:2518312/93‑1 (MQ=255)
                                                   gCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGTGGCACGTCCTGAACTGGCTGAAATTCGTAAACCGTCAATTGCCGGCCTGAtgt                                    >  2:817431/1‑92 (MQ=255)
                                                     gTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGTGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGAtggtg                                  <  1:3500989/93‑1 (MQ=255)
                                                     gTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGAtggtg                                  <  1:1545116/93‑1 (MQ=255)
                                                     gTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGAtggtg                                  <  1:2696311/93‑1 (MQ=255)
                                                     gTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTACTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGGCTGAtggtg                                  >  2:4220439/1‑93 (MQ=255)
                                                      tttACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCCGAAATTCGGAACCCCTCCATCGCCGGTCCGAtggggg                                 >  1:3980599/1‑93 (MQ=255)
                                                       ttACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGa                                >  1:3399529/1‑93 (MQ=255)
                                                        tACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGaa                               >  1:3846338/1‑93 (MQ=255)
                                                          cTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAat                             <  1:1942410/93‑1 (MQ=255)
                                                           ttGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAatg                            <  1:1421274/93‑1 (MQ=255)
                                                            tGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAatga                           <  1:4642990/93‑1 (MQ=255)
                                                            tGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAatga                           >  2:3716871/1‑93 (MQ=255)
                                                            tGCTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAatga                           >  1:39855/1‑93 (MQ=255)
                                                              cTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGNAatgatg                         <  1:2941685/93‑1 (MQ=255)
                                                              cTGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAatgatg                         <  2:4434456/93‑1 (MQ=255)
                                                               tGGTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAATGATgg                        <  1:3205278/93‑1 (MQ=255)
                                                                 gTAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTACCCCGTCAATGGCCGGTCTGATGGTGGAAATGATgggg                      <  2:3473770/93‑2 (MQ=255)
                                                                  tAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAATGATggggc                     <  1:431258/93‑3 (MQ=255)
                                                                  tAACGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAATGATggggc                     <  1:2867418/93‑3 (MQ=255)
                                                                   aaCGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAATGATGGGGcc                    >  1:2851819/1‑93 (MQ=255)
                                                                    aCGCTGCTTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAATGATGGGGCCa                   <  1:753850/93‑1 (MQ=255)
                                                                      gctgctTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAATGATGGGGCCATg                 <  1:2263401/93‑1 (MQ=255)
                                                                      cctgctTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAATGATGGGGCCATg                 <  1:560248/92‑1 (MQ=255)
                                                                       ctgctTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAATGATGGGGCCATgg                >  2:3228795/1‑93 (MQ=255)
                                                                       ctgctTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAATGATGGGGCCATgg                <  2:1414576/93‑1 (MQ=255)
                                                                       ctgctTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAATGATGGGGCCATgg                <  1:2305426/93‑1 (MQ=255)
                                                                        tgctTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAATGATGGGGCCATggg               >  2:4226585/1‑93 (MQ=255)
                                                                        tgctTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAATGATGGGGCCATggg               <  2:4390730/93‑1 (MQ=255)
                                                                        tgctTTCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAATGAAGGGGCCATggg               >  2:3978269/1‑93 (MQ=255)
                                                                           tttCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAATGATGGGGCCATGGGGCg            >  1:404610/1‑93 (MQ=255)
                                                                             tCGCTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAATGATGGGGCCATGGGGCGaa          <  1:1155289/93‑1 (MQ=255)
                                                                                cTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAATGATGGGGCCATGGGGCGAAatc       >  2:3554163/1‑93 (MQ=255)
                                                                                cTGGGTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAATGATGGGGCCATGGGGCGAAatc       >  2:2271774/1‑93 (MQ=255)
                                                                                  gggTGTGGCGGCGCGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAATGATGGGGCCATGGGGCGAAatcat     <  1:2944905/93‑1 (MQ=255)
                                                                                  gggTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAATGATGGGGCCATGGGGCGCAatcat     >  2:98001/1‑93 (MQ=255)
                                                                                  gggTGTGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAATGATGGGGCCATGGGGCGAAatcat     <  2:1338472/93‑1 (MQ=255)
                                                                                    gtgtGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAATGATGGGGCCATGGGGCGAAatcatcc   >  2:3276459/1‑92 (MQ=255)
                                                                                    gtgtGGCGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAATGATGGGGCCATGGGGCGAAatcatca   <  1:4103202/93‑1 (MQ=255)
                                                                                     tgtgGCGGCACGTCTTGATCGGGCTGAACTTCGTTACCGGTCGATGGCCGGTCTGATGGGGGAAATGATGGGTCCATGGGGCGAAatcatcat  <  2:433596/93‑1 (MQ=255)
                                                                                          |                                                                                       
ATAAACGTGATGTTGGCAAAGCGACACTGCTGGCGGTTCTCTCCGCTCTGGGCGTTTACTTGCTGGTAACGCTGCTTTCGCTGGGTGTGGTGGCACGTCCTGAACTGGCTGAAATTCGTAACCCGTCAATGGCCGGTCTGATGGTGGAAATGATGGGACCATGGGGCGAAATCATCAT  >  NC_000913/1678257‑1678434

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑