Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,803,104 A→G 24.8% pseudogene (11/426 nt) arpB → pseudogene, ankyrin repeats

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,803,1040AG24.8% 220.5 / 84.7 129pseudogene (11/426 nt)arpBpseudogene, ankyrin repeats
Reads supporting (aligned to +/- strand):  ref base A (40/57);  new base G (14/18);  total (54/75)
Fisher's exact test for biased strand distribution p-value = 8.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

Only 100 of 138 total aligned reads displayed.

TATAACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCG  >  NC_000913/1803005‑1803203
                                                                                                   |                                                                                                   
tataACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCaaaa                                                                                                     <  1:2557563/100‑1 (MQ=255)
 ataACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAAc                                                                                                    <  2:2123296/100‑1 (MQ=255)
 ataACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAAc                                                                                                    <  1:2000591/100‑1 (MQ=255)
 ataACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAAc                                                                                                    >  2:1879332/1‑100 (MQ=255)
  taACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACg                                                                                                   >  1:2109900/1‑100 (MQ=255)
  taACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATAGAGAAGGCGAACAGATTATGAGTCAAAACg                                                                                                   <  2:1341031/100‑1 (MQ=255)
  taAAAATTTTTCGTGCACATGTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAAc                                                                                                    <  2:956263/99‑1 (MQ=255)
   aaCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGa                                                                                                  <  1:3171455/100‑1 (MQ=255)
    aCAATTTTTCGTGCACATTTGGCTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGa                                                                                                  >  2:1349810/1‑99 (MQ=255)
      aaTTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGatat                                                                                               <  1:2452498/100‑1 (MQ=255)
          tttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCAtt                                                                                           <  1:1768443/100‑1 (MQ=255)
           ttCGTGCACATTTCGTTCTGGCAATAATAATTAATCATCATCCTCAAAACCTCAAGATATATATAAGGCGAACAGATTATGAGTCATAACGATATCATTa                                                                                          <  2:2026774/100‑1 (MQ=255)
                 cacaTTTCGTTCTGGCGATACTAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGaa                                                                                    >  1:3076123/1‑100 (MQ=255)
                    aTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTc                                                                                 >  1:1211287/1‑100 (MQ=255)
                      ttCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCAt                                                                               <  1:3115352/100‑1 (MQ=255)
                       tCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGAGTATGAGTCAAAACGATATCATTATCAGAACTCAtt                                                                              <  1:83384/100‑1 (MQ=255)
                        cGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATta                                                                             <  1:3235876/100‑1 (MQ=255)
                         gTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATtat                                                                            <  2:1517769/100‑1 (MQ=255)
                         gTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATtat                                                                            <  1:1000945/100‑1 (MQ=255)
                          ttCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATtata                                                                           <  1:669887/100‑1 (MQ=255)
                             tGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGt                                                                        <  1:2049349/100‑1 (MQ=255)
                              ggCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGtc                                                                       >  2:2322484/1‑100 (MQ=255)
                                cGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGtctc                                                                     >  2:1598112/1‑100 (MQ=255)
                                 gATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTctcc                                                                    <  1:2635908/100‑1 (MQ=255)
                                 gATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTctcc                                                                    <  1:70976/100‑1 (MQ=255)
                                  ataataATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTctcct                                                                   >  1:1545892/1‑100 (MQ=255)
                                  ataataATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTctcct                                                                   <  1:1630957/100‑1 (MQ=255)
                                  ataataATTAATCATCATCCGCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTctcc                                                                    <  2:686309/99‑1 (MQ=255)
                                    aataatTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCa                                                                 >  1:1154904/1‑100 (MQ=255)
                                     ataatTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCAt                                                                <  2:2541630/100‑1 (MQ=255)
                                     ataatTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCAt                                                                >  1:3077377/1‑100 (MQ=255)
                                      taattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGATACAAAACGATATCATTATCAGAAATCATTATAAGTCTCCTCTTa                                                               >  1:242586/1‑100 (MQ=255)
                                      taattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATa                                                               <  2:2696840/100‑1 (MQ=255)
                                      taattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATa                                                               >  1:1066453/1‑100 (MQ=255)
                                      taattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATa                                                               <  2:1548215/100‑1 (MQ=255)
                                       aattattCATCATCCTCAAAACCTCCAGATATATAAAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATa                                                               <  1:1015385/99‑1 (MQ=255)
                                       aattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCAt                                                                >  2:1975338/1‑98 (MQ=255)
                                       aattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAg                                                              <  1:252057/100‑1 (MQ=255)
                                         ttaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGAt                                                            <  2:1820274/100‑1 (MQ=255)
                                          taatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGAtt                                                           <  2:1435199/100‑1 (MQ=255)
                                           aatCATCATCCTCAAAACCTCCAGATATATATAAGGTGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTg                                                          <  2:1809138/100‑1 (MQ=255)
                                            atcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGc                                                         <  2:1448924/100‑1 (MQ=255)
                                            atcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGc                                                         >  1:150728/1‑100 (MQ=255)
                                            atcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGc                                                         <  1:2129256/100‑1 (MQ=255)
                                            atcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGc                                                         <  1:2794616/100‑1 (MQ=255)
                                            atcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGc                                                         <  1:727890/100‑1 (MQ=255)
                                             tcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGca                                                        >  2:2928531/1‑100 (MQ=255)
                                             tcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGca                                                        <  2:1916667/100‑1 (MQ=255)
                                             tcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGca                                                        <  1:2617549/100‑1 (MQ=255)
                                             tcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGca                                                        <  2:15152/100‑1 (MQ=255)
                                              catcatcCTCAAAACCTCCAGATATATATAATGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGca                                                        <  2:793190/99‑1 (MQ=255)
                                               atcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGcaca                                                      >  1:258083/1‑100 (MQ=255)
                                                tcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACAt                                                     >  1:546733/1‑100 (MQ=255)
                                                 catcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATc                                                    >  2:2005214/1‑100 (MQ=255)
                                                 catcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATc                                                    >  2:2048827/1‑100 (MQ=255)
                                                    ccTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGAt                                                 >  2:493000/1‑100 (MQ=255)
                                                    ccTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGAt                                                 <  2:2784031/100‑1 (MQ=255)
                                                       cAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGc                                              <  2:2673769/100‑1 (MQ=255)
                                                       cAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGc                                              <  1:2182406/100‑1 (MQ=255)
                                                        aaaaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCg                                             >  1:1716344/1‑100 (MQ=255)
                                                        aaaaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCg                                             <  1:2038519/100‑1 (MQ=255)
                                                        aaaaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCg                                             <  2:2828090/100‑1 (MQ=255)
                                                         aaaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGa                                            <  2:1145137/100‑1 (MQ=255)
                                                          aaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGAc                                           <  1:379222/100‑1 (MQ=255)
                                                              tCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATAc                                       <  1:2817510/100‑1 (MQ=255)
                                                               ccAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATAcc                                      >  1:1941420/1‑100 (MQ=255)
                                                               ccAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATAcc                                      <  1:958773/100‑1 (MQ=255)
                                                                cAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCa                                     >  1:2736274/1‑100 (MQ=255)
                                                                  gATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAac                                   <  1:2719679/100‑1 (MQ=255)
                                                                  gATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAac                                   <  2:2616364/100‑1 (MQ=255)
                                                                  gATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAac                                   <  2:1463431/100‑1 (MQ=255)
                                                                     atatataAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACAcc                                >  2:2923070/1‑100 (MQ=255)
                                                                      tatataAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCt                               >  1:440859/1‑100 (MQ=255)
                                                                       atataAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCtt                              <  2:647495/100‑1 (MQ=255)
                                                                       atataAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCtt                              >  2:1287060/1‑100 (MQ=255)
                                                                        tataAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtc                             >  1:2617811/1‑100 (MQ=255)
                                                                           aaGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACAGCTtcatc                          <  1:1715923/100‑1 (MQ=255)
                                                                           aaGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcatc                          <  1:858124/100‑1 (MQ=255)
                                                                           aaGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcatc                          <  1:2857881/100‑1 (MQ=255)
                                                                            aGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcatca                         <  1:549736/100‑1 (MQ=255)
                                                                            aGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcatca                         <  2:2815948/100‑1 (MQ=255)
                                                                             ggCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCag                        >  1:1044642/1‑100 (MQ=255)
                                                                             ggCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCag                        >  1:2033246/1‑100 (MQ=255)
                                                                             ggCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCag                        >  1:261461/1‑100 (MQ=255)
                                                                               cGAACAGATTATGAGTGAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCagag                      <  2:3077359/100‑1 (MQ=255)
                                                                               cGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCagag                      >  2:2563068/1‑100 (MQ=255)
                                                                                gAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGc                     <  1:794154/100‑1 (MQ=255)
                                                                                   cAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTa                  >  1:1916863/1‑100 (MQ=255)
                                                                                        tATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAAt             >  1:1983661/1‑100 (MQ=255)
                                                                                          tGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCa           <  2:1545888/100‑1 (MQ=255)
                                                                                            aGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAt         >  2:1765697/1‑100 (MQ=255)
                                                                                               cAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg      >  1:2527395/1‑100 (MQ=255)
                                                                                               cAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg      >  2:2425218/1‑100 (MQ=255)
                                                                                                aaaaCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTgc     <  1:3004761/100‑1 (MQ=255)
                                                                                                aaaaCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGAAATACCAACACCTTCATCAGAGCCTATTAATCAATTTgc     >  2:2065206/1‑100 (MQ=255)
                                                                                                   aCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTgcgcg  <  2:399573/100‑1 (MQ=255)
                                                                                                   aCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTgcgcg  <  1:878688/100‑1 (MQ=255)
                                                                                                   aCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTgcgcg  <  1:667483/100‑1 (MQ=255)
                                                                                                   aCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTgcgcg  >  1:556152/1‑100 (MQ=255)
                                                                                                   aCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTgcgcg  <  2:946780/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TATAACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCGCG  >  NC_000913/1803005‑1803203

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