Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,830,953 T→G N344H (AAT→CAT)  astC ← succinylornithine transaminase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,830,9530TG97.5% 404.1 / ‑2.6 120N344H (AAT→CAT) astCsuccinylornithine transaminase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base G (50/67);  total (53/67)
Fisher's exact test for biased strand distribution p-value = 8.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

CAAAACGCACCACGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATTCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTCAACAAACCAGTCGT  >  NC_000913/1830857‑1831048
                                                                                                |                                                                                               
 aaaaCGCACCACGTTGCCACCCGCAATCAGTACCATCACGCCCGCTTTCGCCGCTTCCTACGCCATCTGGTTCGCTTGCCCGCCGTACTCGTCACTCAGt                                                                                             >  2:2760453/1‑100 (MQ=255)
   aaCGCACCACGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATTCAGTac                                                                                           >  1:1853875/1‑100 (MQ=255)
   aaCGCACCACGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTac                                                                                           <  1:1976828/100‑1 (MQ=255)
     cGCACCACGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCGTGCAGTacac                                                                                         <  1:1035896/100‑1 (MQ=255)
      gccccacGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTAcaca                                                                                        <  1:2034798/97‑1 (MQ=255)
      gCACCACGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTAcaca                                                                                        <  1:2757161/100‑1 (MQ=255)
      gCACCACGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTAcaca                                                                                        >  2:1044510/1‑100 (MQ=255)
       caccacGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCACTACACAg                                                                                       >  2:2791289/1‑100 (MQ=255)
       caccacGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTCCGCCGCTACCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAg                                                                                       >  1:3870101/1‑100 (MQ=255)
        tccacGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGc                                                                                      <  2:3550454/99‑1 (MQ=255)
        accacGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGc                                                                                      <  1:2673560/100‑1 (MQ=255)
         ccacGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGcc                                                                                     >  2:3736113/1‑100 (MQ=255)
         ccacGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGca                                                                                     >  1:3628309/1‑99 (MQ=255)
          cacGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCa                                                                                    >  2:727295/1‑100 (MQ=255)
          cacGTTGCCACCCGCAATCAGGACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCa                                                                                    <  2:2562293/100‑1 (MQ=255)
           acGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCaa                                                                                   >  1:730222/1‑100 (MQ=255)
              ttGCCACCCGCAATCAGTACCATCGCGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCa                                                                                <  2:257266/100‑1 (MQ=255)
              ttGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCa                                                                                <  1:718698/100‑1 (MQ=255)
               tGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAg                                                                               <  2:231267/100‑1 (MQ=255)
                 ccACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCa                                                                             >  1:471175/1‑100 (MQ=255)
                  cACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCaa                                                                            >  1:3340871/1‑100 (MQ=255)
                   aCCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCaaa                                                                           >  2:3980385/1‑100 (MQ=255)
                    cccGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAAc                                                                          <  2:1611162/100‑1 (MQ=255)
                    cccGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAAc                                                                          <  2:3311670/100‑1 (MQ=255)
                    cccGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAAc                                                                          >  2:544530/1‑100 (MQ=255)
                    cccGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAAc                                                                          >  1:1014067/1‑100 (MQ=255)
                    cccGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAAc                                                                          >  1:1767923/1‑100 (MQ=255)
                     ccGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAAcc                                                                         <  2:3307359/100‑1 (MQ=255)
                      cGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCt                                                                        >  2:572038/1‑100 (MQ=255)
                       gCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCtc                                                                       >  1:4084820/1‑99 (MQ=255)
                         aaTCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAg                                                                     <  1:583644/100‑1 (MQ=255)
                         aaTCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAg                                                                     <  1:3435860/100‑1 (MQ=255)
                         aaTCAGTACCATCACGCCTGCTTTCGCCGCTCCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAg                                                                     <  2:3150814/100‑1 (MQ=255)
                         aaTCAGTACCATCACGCCTGCTTTCGCCGCCTCCTGAGAGATCTGTTGCGCTTGCCCGGCTCCATCGGCATGCAGTACACAGCCAAGCAGCACACCTAAg                                                                     <  2:2268251/100‑1 (MQ=255)
                          aTCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGc                                                                    <  1:463274/100‑1 (MQ=255)
                          aTCAATACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGc                                                                    >  2:2672452/1‑100 (MQ=255)
                            cAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCcg                                                                  >  1:2268045/1‑100 (MQ=255)
                             aGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCcgc                                                                 >  2:2922529/1‑100 (MQ=255)
                              gTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCcgcg                                                                >  1:3578735/1‑100 (MQ=255)
                                aCCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGaa                                                              >  2:1645233/1‑100 (MQ=255)
                                 ccATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTGCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAAc                                                             <  1:2824608/100‑1 (MQ=255)
                                    tCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTc                                                          <  1:1668706/100‑1 (MQ=255)
                                    tCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTc                                                          <  1:2330810/100‑1 (MQ=255)
                                      aCGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCAc                                                        <  1:907876/100‑1 (MQ=255)
                                      aCGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCAc                                                        <  2:2633405/100‑1 (MQ=255)
                                       cGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACt                                                       <  2:1302851/100‑1 (MQ=255)
                                       cGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACt                                                       >  2:2801681/1‑100 (MQ=255)
                                        gCCTGGTCTCGCCGCTTCCGGAGAGATTTGTTTCGTCTGCCCGGCGTACTCGGCATGTAGTACACAGCCAATTAACAAATCTACGCCGCGAACTTCACTg                                                      <  2:3639792/100‑1 (MQ=255)
                                        gCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTg                                                      <  1:1396416/100‑1 (MQ=255)
                                        gCCTGCTTTCGCCGCTTCCCGAGAGATCTCTTTCGCTTGCCCCGCGTAATCGGCATGCAGTACACAGCCAATCCGCAATCCTAAGCCGCGTACTTCACTg                                                      <  2:937838/100‑1 (MQ=255)
                                          cTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTgca                                                    >  1:3421220/1‑98 (MQ=255)
                                           tGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAAc                                                   <  2:338963/100‑1 (MQ=255)
                                           tGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAAc                                                   <  1:2945414/100‑1 (MQ=255)
                                           tGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAAc                                                   >  1:2449864/1‑100 (MQ=255)
                                             cTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACaa                                                 >  1:1093233/1‑100 (MQ=255)
                                              tttCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACaaa                                                <  1:583875/100‑1 (MQ=255)
                                               ttCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTCAGCCGCGAACTTCACTGAACAAAc                                               >  1:1006298/1‑100 (MQ=255)
                                               ttCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAAc                                               <  2:2733122/100‑1 (MQ=255)
                                                tCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAAcc                                              <  1:911113/100‑1 (MQ=255)
                                                 cgccgcTCCCTCAAAGATCTGTTTCGCTCGCCCGACGGACTCGGCATGCAGTACACAGCCAAACCGCAAACCTCAGCCGCGCACTTCACTGAACAAAccc                                             >  2:1189887/1‑99 (MQ=255)
                                                  gccgcTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCAt                                            >  1:2382686/1‑100 (MQ=255)
                                                  gccgcTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCAt                                            >  1:3097512/1‑100 (MQ=255)
                                                    cgcTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAg                                          >  1:3305993/1‑100 (MQ=255)
                                                     gcgtCCTGAGAGATCTGTGTCGCTTGCCCTGCGTAATCGGCATGCAGTACACATCCAATCAGCAAACCTAAGCCGCTAACTTCACTGAACAAACCATAGc                                         <  2:3360674/97‑1 (MQ=255)
                                                     gcTTCCTGAGCGATCTGTTTCGCTTGCCCGGCGTAATCGGCTTGCAGTACACAGCCAATCAGCAAACCTATGCCGCGAACTTCACTGAACCAACCATAGc                                         <  2:580596/100‑1 (MQ=255)
                                                     gcTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGc                                         <  2:2352376/100‑1 (MQ=255)
                                                      cTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCg                                        >  2:973932/1‑100 (MQ=255)
                                                        tCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGt                                      <  1:2769984/100‑1 (MQ=255)
                                                         ccTGAGAGATCTGGTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTg                                     <  2:2850013/100‑1 (MQ=255)
                                                          cTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGa                                    >  2:2161734/1‑100 (MQ=255)
                                                           tGAGAGATCTGTTTCGCTTGCCCGGCGTAAGCGGCATGCAGTACACAGACAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGAt                                   <  2:1827365/100‑1 (MQ=255)
                                                            gagagaTCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGAtt                                  <  2:3643781/100‑1 (MQ=255)
                                                             agagaTCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTa                                 <  2:3051302/100‑1 (MQ=255)
                                                              gagaTCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTaa                                <  2:663992/100‑1 (MQ=255)
                                                               agaTCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAAt                               <  1:421582/100‑1 (MQ=255)
                                                                 aTCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACACCCCATAGCGGTGACTCATAg                             >  2:3304862/1‑100 (MQ=255)
                                                                 aTCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAg                             <  1:568386/100‑1 (MQ=255)
                                                                 aTCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAg                             <  1:721468/100‑1 (MQ=255)
                                                                 aTCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAg                             >  2:736450/1‑100 (MQ=255)
                                                                 aTCTGGTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAg                             <  2:3328148/100‑1 (MQ=255)
                                                                 aTCTGCTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAg                             <  2:1518933/100‑1 (MQ=255)
                                                                   cTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTa                           <  1:2276262/100‑1 (MQ=255)
                                                                    tGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTAt                          <  1:2927805/100‑1 (MQ=255)
                                                                          gCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGAc                    >  2:1473457/1‑100 (MQ=255)
                                                                             tGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCt                 >  2:2652941/1‑100 (MQ=255)
                                                                              gCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCAGTGAACACACCATAGCGGTGATTAACAGTATAAAGACGcgc                >  2:165727/1‑98 (MQ=255)
                                                                              gCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTc                <  1:3410979/100‑1 (MQ=255)
                                                                               cccGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTca               <  1:1799707/100‑1 (MQ=255)
                                                                               cccGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTca               >  2:2247359/1‑100 (MQ=255)
                                                                                ccGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTcaa              >  2:1961000/1‑100 (MQ=255)
                                                                                  ggCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTcaaca            <  2:2882262/100‑1 (MQ=255)
                                                                                     gTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTCAACAAAc         >  2:3414501/1‑100 (MQ=255)
                                                                                     gTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTCAACAAAc         <  1:2495645/100‑1 (MQ=255)
                                                                                       aaTCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTCAACAAACCa       >  1:3100580/1‑100 (MQ=255)
                                                                                       aaTCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTCAACAAACCa       >  1:1096465/1‑100 (MQ=255)
                                                                                        aTCGGCATTCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTCAACAAACCAg      >  2:2806363/1‑100 (MQ=255)
                                                                                          cGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTCAACAAACCAGTc    >  1:2414393/1‑100 (MQ=255)
                                                                                           ggCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTCAACAAACCAGTCg   <  1:1371866/100‑1 (MQ=255)
                                                                                           ggCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTCAACAAACCAGTCg   >  1:2073236/1‑100 (MQ=255)
                                                                                            gCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTCAACAAACCAGTCGt  >  1:1516430/1‑100 (MQ=255)
                                                                                                |                                                                                               
CAAAACGCACCACGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATTCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTCAACAAACCAGTCGT  >  NC_000913/1830857‑1831048

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: