Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,827,031 T→A D410V (GAT→GTT)  astB ← succinylarginine dihydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,827,0310TA98.9% 305.6 / ‑3.9 87D410V (GAT→GTT) astBsuccinylarginine dihydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (42/44);  total (42/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

CTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCACCGC  >  NC_000913/1826934‑1827125
                                                                                                 |                                                                                              
ctcGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACg                                                                                              <  1:2316143/100‑1 (MQ=255)
ctcGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACg                                                                                              <  2:282433/100‑1 (MQ=255)
  cgcgCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGgc                                                                                            <  1:2144892/100‑1 (MQ=255)
  cgcgCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGgc                                                                                            <  2:377408/100‑1 (MQ=255)
  cgcgCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGgc                                                                                            <  1:1625409/100‑1 (MQ=255)
   gcgcTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGgca                                                                                           >  1:910098/1‑100 (MQ=255)
    cgcTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCACTTCGCGCAGCAATTGCGGGTCGGCCAGAACGgcag                                                                                          >  1:1641986/1‑100 (MQ=255)
      ntGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGgcagca                                                                                        >  1:1300909/2‑100 (MQ=255)
      cTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGgcagca                                                                                        <  2:1190662/100‑1 (MQ=255)
       tGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGcagcag                                                                                       <  1:1399721/100‑1 (MQ=255)
       tGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGcagcag                                                                                       >  2:82715/1‑100 (MQ=255)
        ggAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGt                                                                                      <  2:290525/100‑1 (MQ=255)
         gAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTa                                                                                     <  2:236958/100‑1 (MQ=255)
         gAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTa                                                                                     <  2:984952/100‑1 (MQ=255)
            cGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAgg                                                                                  <  2:332545/100‑1 (MQ=255)
              gATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCCGCAATTGAGGGTCGGCCAGAACGGCAGCAGTACGGCg                                                                                >  2:1612903/1‑100 (MQ=255)
               aTAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGa                                                                               >  2:254857/1‑100 (MQ=255)
                tAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGAt                                                                              <  1:1016743/100‑1 (MQ=255)
                  aaccgaaccgAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATcg                                                                            >  1:1034212/1‑100 (MQ=255)
                    ccgaaccgaGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATcgcg                                                                          >  1:782020/1‑100 (MQ=255)
                    ccgaaccgaGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATcgcg                                                                          >  2:284527/1‑100 (MQ=255)
                         ccgaGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGgtagt                                                                     >  1:1813502/1‑100 (MQ=255)
                         ccgaGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGgtagt                                                                     <  1:97644/100‑1 (MQ=255)
                          cgaGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGgtagta                                                                    >  1:389581/1‑100 (MQ=255)
                          ccagaTTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGgtagt                                                                     <  1:1669003/97‑1 (MQ=255)
                           gagaTTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTaa                                                                   >  2:500698/1‑100 (MQ=255)
                           gagaTTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTaa                                                                   >  1:801148/1‑100 (MQ=255)
                            agaTTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAAc                                                                  <  1:163878/100‑1 (MQ=255)
                             gaTTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACg                                                                 <  2:140778/100‑1 (MQ=255)
                             gaTTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACg                                                                 >  2:633189/1‑100 (MQ=255)
                             gaTTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACg                                                                 >  1:2265065/1‑100 (MQ=255)
                              aTTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGa                                                                <  2:451025/100‑1 (MQ=255)
                               ttCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGAt                                                               >  2:150412/1‑100 (MQ=255)
                                 cAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATcc                                                             >  1:1290944/1‑100 (MQ=255)
                                 cAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATcc                                                             >  1:1227241/1‑100 (MQ=255)
                                  aGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCa                                                            >  1:411134/1‑100 (MQ=255)
                                    tAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCAcc                                                          <  2:2279506/100‑1 (MQ=255)
                                     aaTTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCAccc                                                         <  2:1290941/100‑1 (MQ=255)
                                        tGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAg                                                       <  1:683633/99‑1 (MQ=255)
                                        tGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGt                                                      <  1:1786759/100‑1 (MQ=255)
                                          tcTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCa                                                    <  1:447135/99‑1 (MQ=255)
                                          gCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCa                                                    <  1:1909137/100‑1 (MQ=255)
                                           cTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCAt                                                   >  2:1718231/1‑100 (MQ=255)
                                             cAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCaaga                                                 >  2:447133/1‑97 (MQ=255)
                                             cAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTg                                                 >  1:964843/1‑100 (MQ=255)
                                             cAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTg                                                 <  2:1009532/100‑1 (MQ=255)
                                             cAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTg                                                 >  1:1190748/1‑100 (MQ=255)
                                              aaTACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGa                                                >  2:2273526/1‑100 (MQ=255)
                                               aTACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAg                                               <  1:1956116/100‑1 (MQ=255)
                                                 aCATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAgcg                                             <  2:1034210/100‑1 (MQ=255)
                                                   aTCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCa                                           >  2:1232079/1‑100 (MQ=255)
                                                     ccAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCAtt                                         >  2:876910/1‑100 (MQ=255)
                                                     ccAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCAtt                                         >  2:1734864/1‑100 (MQ=255)
                                                     ccAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCAtt                                         <  1:500700/100‑1 (MQ=255)
                                                      cAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTa                                        <  2:910096/100‑1 (MQ=255)
                                                      cAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTa                                        >  1:1460716/1‑100 (MQ=255)
                                                      cAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTa                                        >  1:1299737/1‑100 (MQ=255)
                                                       aGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTaa                                       >  1:1515997/1‑100 (MQ=255)
                                                       aGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTaa                                       >  1:404650/1‑100 (MQ=255)
                                                         tGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAATGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTaaan                                     <  2:1300906/100‑2 (MQ=255)
                                                         tGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAAc                                     <  2:2316158/100‑1 (MQ=255)
                                                          gCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAAc                                     <  2:2217554/99‑1 (MQ=255)
                                                             tCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCg                                 <  1:254859/100‑1 (MQ=255)
                                                              cccGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGt                                >  1:2106681/1‑100 (MQ=255)
                                                              cccGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGt                                >  2:666150/1‑100 (MQ=255)
                                                              cccGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGt                                <  2:964841/100‑1 (MQ=255)
                                                               ccGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTa                               >  2:24611/1‑100 (MQ=255)
                                                                   ccTTCGCGCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGt                           >  1:2008071/1‑100 (MQ=255)
                                                                       cgcgCAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtcat                       >  2:1859190/1‑100 (MQ=255)
                                                                         cgcAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtcatca                     <  1:1037034/100‑1 (MQ=255)
                                                                         cgcAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtcatca                     <  1:905277/100‑1 (MQ=255)
                                                                         cgcAGCAATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtcatca                     <  1:876912/100‑1 (MQ=255)
                                                                           cagcaATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtcatcatc                   >  1:371431/1‑100 (MQ=255)
                                                                           cagcaATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTtcatcat                    <  1:464384/99‑1 (MQ=255)
                                                                             gcaATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCAc                 <  2:2186517/100‑1 (MQ=255)
                                                                             gcaATTGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCAc                 >  2:2174283/1‑100 (MQ=255)
                                                                                  tGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCAccgccg            >  1:1313073/1‑100 (MQ=255)
                                                                                  tGCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCAccgccg            <  2:1813497/100‑1 (MQ=255)
                                                                                   gCGGGTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCgg           >  1:1663898/1‑100 (MQ=255)
                                                                                     gggTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGAt         <  1:300507/100‑1 (MQ=255)
                                                                                      ggTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGAtt        >  2:1207080/1‑100 (MQ=255)
                                                                                      ggTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGAtt        <  1:2273534/100‑1 (MQ=255)
                                                                                       gTCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTc       <  2:2008064/100‑1 (MQ=255)
                                                                                        tCGGCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCa      <  2:801147/100‑1 (MQ=255)
                                                                                          ggCCAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCAcc    >  2:1312991/1‑100 (MQ=255)
                                                                                            ccAGAACGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCACCGc  <  1:284529/100‑1 (MQ=255)
                                                                                            ccAGAACGGCAGCAGTAAGGCGATCGCGGTAGCAACGATGCATCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCACCGc  <  2:1641981/100‑1 (MQ=255)
                                                                                                 |                                                                                              
CTCGCGCTGGAACGGATAAACCGAACCGAGATTCAGTAATTGGCTCAATACATCCAGTGCTTCCCGCCCTTCGCGCAGCAATTGCGGGTCGGCCAGATCGGCAGCAGTAAGGCGATCGCGGTAGTAACGATCCACCCAGTCATTGAGCGCATTAAACAGCGTATCGTTCATCATCACCGCCGGATTCACCGC  >  NC_000913/1826934‑1827125

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