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.2% 741.0 / 3.6 214N344H (AAT→CAT) astCsuccinylornithine transaminase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base G (93/115);  total (96/118)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 214 total aligned reads displayed.

CAAAACGCACCACGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATTCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTCAACAAACCAGTCGTG  >  NC_000913/1830857‑1831049
                                                                                                |                                                                                                
cAAAACGCACCACGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCan                                                                                               <  1:4426647/100‑2 (MQ=255)
cAAAACGCACCACGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAg                                                                                               >  2:5163451/1‑100 (MQ=255)
cAAAACGCACCACGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAg                                                                                               <  1:3307655/100‑1 (MQ=255)
cAAAACGCACCACGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAg                                                                                               >  1:1699626/1‑100 (MQ=255)
 aaaaCGCACCACGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGt                                                                                              >  1:2825319/1‑100 (MQ=255)
 aaaaCGCACCACGTTCCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGt                                                                                              <  1:4942204/100‑1 (MQ=255)
   aaCGCACCACGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTac                                                                                            <  1:1572448/100‑1 (MQ=255)
   aaCGCACCACGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCAGCTTCCAGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTac                                                                                            <  1:5393400/100‑1 (MQ=255)
      gCACCACGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTAcaca                                                                                         >  2:4871767/1‑100 (MQ=255)
       caccacGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAg                                                                                        >  2:3846133/1‑100 (MQ=255)
        accacGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGc                                                                                       <  1:3510051/100‑1 (MQ=255)
         ccacGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGcc                                                                                      >  1:3932869/1‑100 (MQ=255)
         ccacGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGGTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGcc                                                                                      <  2:2931169/100‑1 (MQ=255)
             gTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATc                                                                                  <  1:2892130/100‑1 (MQ=255)
              ttGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCa                                                                                 >  2:1487457/1‑100 (MQ=255)
               tGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAg                                                                                >  2:3032076/1‑100 (MQ=255)
                  cACCCGAAATCAGTACCATCACGCCTGGTTTCGCCGCTTCCTGCGAGATCTGTTTCGCTTGCCCGTCGTAGTCTCCATGCAGTACACAGCCAATCAGCaa                                                                             >  2:5081601/1‑100 (MQ=255)
                   aCCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCaaa                                                                            <  2:1837998/100‑1 (MQ=255)
                    cccGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATTCAGTACACAGCCAATCAGCAAAc                                                                           >  1:4217235/1‑100 (MQ=255)
                    cccGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAAc                                                                           <  1:192336/100‑1 (MQ=255)
                     ccGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAAcc                                                                          <  2:3536631/100‑1 (MQ=255)
                     ccGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAAcc                                                                          <  1:658645/100‑1 (MQ=255)
                      cGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCt                                                                         <  1:1338679/100‑1 (MQ=255)
                      cGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCt                                                                         <  2:4074359/100‑1 (MQ=255)
                      cGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCt                                                                         <  1:1506555/100‑1 (MQ=255)
                       gCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTa                                                                        <  2:3588398/100‑1 (MQ=255)
                              gTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCcgcg                                                                 >  1:1848314/1‑100 (MQ=255)
                               tACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGa                                                                >  1:199371/1‑100 (MQ=255)
                                 ccATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAAc                                                              <  1:1011248/100‑1 (MQ=255)
                                 ccATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAAc                                                              >  2:1263839/1‑100 (MQ=255)
                                   aTCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACtt                                                            >  1:2366040/1‑100 (MQ=255)
                                    tCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTc                                                           <  2:5218874/100‑1 (MQ=255)
                                    tCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTc                                                           <  1:3089507/100‑1 (MQ=255)
                                    tCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTc                                                           >  2:854592/1‑100 (MQ=255)
                                     cACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCa                                                          <  2:4370684/100‑1 (MQ=255)
                                      aCGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCAc                                                         <  1:5258354/100‑1 (MQ=255)
                                        gCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTg                                                       <  1:3010390/100‑1 (MQ=255)
                                        gCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTg                                                       <  1:3633325/100‑1 (MQ=255)
                                        gCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTg                                                       <  1:1532723/100‑1 (MQ=255)
                                        gCCTGCTTTCGCCGCTTCCGGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTg                                                       <  1:2113439/100‑1 (MQ=255)
                                         ccTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGa                                                      >  1:686488/1‑100 (MQ=255)
                                           tGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAAc                                                    <  2:5285238/100‑1 (MQ=255)
                                           tGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAAc                                                    <  2:5348023/100‑1 (MQ=255)
                                             cTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCCTGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACaa                                                  >  1:3808955/1‑100 (MQ=255)
                                               ttCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAAc                                                <  2:2369804/100‑1 (MQ=255)
                                                tCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAAcc                                               <  2:61442/100‑1 (MQ=255)
                                                 cgccgcTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCa                                              >  2:4761598/1‑100 (MQ=255)
                                                    cgcTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAg                                           >  2:4697782/1‑100 (MQ=255)
                                                    cgcTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAg                                           <  1:3784062/100‑1 (MQ=255)
                                                       ttCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCgg                                        >  2:3547886/1‑100 (MQ=255)
                                                         ccTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTg                                      >  2:4453121/1‑100 (MQ=255)
                                                         ccTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTg                                      <  2:2512032/100‑1 (MQ=255)
                                                          cTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGa                                     >  1:138660/1‑100 (MQ=255)
                                                          cTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGa                                     >  1:3902409/1‑100 (MQ=255)
                                                           tGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGAt                                    >  1:4577519/1‑100 (MQ=255)
                                                           tGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGAt                                    <  2:2409779/100‑1 (MQ=255)
                                                            gagagaTCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGAtt                                   >  2:194350/1‑100 (MQ=255)
                                                            gagagaTCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGAtt                                   <  2:1999659/100‑1 (MQ=255)
                                                              tagaACTGATTAGCTTGTTCTGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTaa                                 >  2:652047/2‑100 (MQ=255)
                                                              gagaTCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTaa                                 <  1:3060609/100‑1 (MQ=255)
                                                               agaTCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAAt                                <  1:1414112/100‑1 (MQ=255)
                                                               agaTCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAAt                                <  2:3037735/100‑1 (MQ=255)
                                                               agaTCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAAt                                <  2:4779055/100‑1 (MQ=255)
                                                                gaTCTGTTTTGCTTGCCCGGCGTCATCGGCATGCAGTACACAACCAATCAGCAATCCTAAGCCGCGAACGTCACTGAACAAACCATAGCGGTGATTAATa                               <  1:4851726/100‑1 (MQ=255)
                                                                gaTCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATa                               >  2:3879078/1‑100 (MQ=255)
                                                                gaTCTGTTTCGCTCGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCCGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATATCGGCGATTGAtc                               >  2:3148966/1‑99 (MQ=255)
                                                                 aTCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAg                              <  2:3955628/100‑1 (MQ=255)
                                                                 aTCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAg                              <  1:1052446/100‑1 (MQ=255)
                                                                 aTCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAg                              <  1:3255800/100‑1 (MQ=255)
                                                                 aTCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAg                              <  1:3503299/100‑1 (MQ=255)
                                                                 aTCTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGAGAACTTCACTGAACAAACCATAGCGGTGATTAATAg                              <  2:4932535/100‑1 (MQ=255)
                                                                   cTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTa                            <  2:1672626/100‑1 (MQ=255)
                                                                   cTGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTa                            <  1:4319199/100‑1 (MQ=255)
                                                                    tgtggcGCTTGCCCGGCGTAATCGGCATTCAGTACACAGCCAATCAGCACACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTAt                           <  2:5406720/95‑1 (MQ=255)
                                                                    tGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTAt                           <  1:4518786/100‑1 (MQ=255)
                                                                    tGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTAt                           <  1:97770/100‑1 (MQ=255)
                                                                    tGTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTAt                           >  1:5369701/1‑100 (MQ=255)
                                                                    tGTTTCGCTTGCCCGGCGTAATCGGCAGGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTAt                           >  2:1524970/1‑100 (MQ=255)
                                                                     gTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTAtt                          <  2:2354759/100‑1 (MQ=255)
                                                                     gTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTAtt                          >  1:4658659/1‑100 (MQ=255)
                                                                     gTTTCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTAtt                          >  1:857959/1‑100 (MQ=255)
                                                                       ttCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTaa                        <  1:1936323/100‑1 (MQ=255)
                                                                        tCGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAACAGCGAACTTCACTGAACAAACCATAGCGGTGATTCATAGTATTAAg                       >  2:428440/1‑100 (MQ=255)
                                                                         cGCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGa                      >  1:22111/1‑100 (MQ=255)
                                                                          gCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGAc                     >  2:2912619/1‑100 (MQ=255)
                                                                          gCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGAc                     >  1:908298/1‑100 (MQ=255)
                                                                          gCTTGCCCGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGAc                     <  1:3172922/100‑1 (MQ=255)
                                                                            ttGCCCGGCGTAATCGGCTTGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGatcg                   >  2:1170090/1‑97 (MQ=255)
                                                                               cccGGCGTAATCGGCATTCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTca                >  2:4195142/1‑100 (MQ=255)
                                                                               cccGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTca                <  1:2497235/100‑1 (MQ=255)
                                                                                 cGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTcaac              >  2:5406101/1‑100 (MQ=255)
                                                                                 cGGCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTcaac              <  2:4528351/100‑1 (MQ=255)
                                                                                  ggCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTcaaca             >  2:170679/1‑100 (MQ=255)
                                                                                   gCGTAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTcaacaa            <  2:3889159/100‑1 (MQ=255)
                                                                                      tAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTCAACAAAcc         >  1:4502830/1‑100 (MQ=255)
                                                                                      tAATCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTCAACAAAcc         <  2:2654133/100‑1 (MQ=255)
                                                                                         tCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTCAACAAACCAGt      >  1:2732516/1‑100 (MQ=255)
                                                                                         tCGGCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTCAACAAACCAGt      <  2:2843877/100‑1 (MQ=255)
                                                                                           ggCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTCAACAAACCAGTCg    <  1:2980998/100‑1 (MQ=255)
                                                                                           ggCATGCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTCAACAAACCAGTCg    <  2:1040853/100‑1 (MQ=255)
                                                                                                |                                                                                                
CAAAACGCACCACGTTGCCACCCGCAATCAGTACCATCACGCCTGCTTTCGCCGCTTCCTGAGAGATCTGTTTCGCTTGCCCGGCGTAATCGGCATTCAGTACACAGCCAATCAGCAAACCTAAGCCGCGAACTTCACTGAACAAACCATAGCGGTGATTAATAGTATTAAGACGCTCAACAAACCAGTCGTG  >  NC_000913/1830857‑1831049

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