Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,703,191 C→T 5.3% A320V (GCT→GTT)  add → adenosine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,703,1910CT 5.3% 256.1 / 3.6 76A320V (GCT→GTT) addadenosine deaminase
Reads supporting (aligned to +/- strand):  ref base C (32/40);  new base T (2/2);  total (34/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

GTGCCAGCATTAACACTGACGATCCCGGCGTACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTC  >  NC_000913/1703047‑1703327
                                                                                                                                                |                                                                                                                                        
gTGCCAGCATTAACACTGACGATCCCGGCGTACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTTAg                                                                                                                                      <  4:998178/149‑1 (MQ=255)
gTGCCAGCATTAACACTGACGATCCCGGCGTACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAg                                                                                                                                      <  3:836721/149‑1 (MQ=255)
  gCCAGCATTAACACTGACGATCCCGGCGTACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGa                                                                                                                                    >  2:1364080/1‑149 (MQ=255)
    cAGCATTAACACTGACGATCCCGGCGTACAGGGAGTGGATATCATTGACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGaaa                                                                                                                                  <  4:1542611/149‑1 (MQ=255)
        aTTAACACTGACGATCCCGGCGTACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATCAATGGTCTAGAAATGGCTTTCCTCAGCGCTGAGGAAAAAcg                                                                                                                              >  4:1549696/1‑149 (MQ=255)
         ttAACACTGACGATCCCGGCGTACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGa                                                                                                                                       >  4:421006/1‑139 (MQ=255)
         ttAACACTGACGATCCCGGCGTACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGa                                                                                                                                       <  3:421006/139‑1 (MQ=255)
         ttAACACTGACGATCCCGGCGTACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAAcgc                                                                                                                             >  2:538279/1‑149 (MQ=255)
             cacTGACGATCCCGGCGTACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCAc                                                                                                                         <  1:1364080/149‑1 (MQ=255)
             cacTGACGATCCCGGCGTACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCAc                                                                                                                         <  1:226900/149‑1 (MQ=255)
              acTGACGATCCCGGCGTACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACt                                                                                                                        <  4:1090383/149‑1 (MQ=255)
                    gATCCCGGCGTACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCgaga                                                                                                                  >  4:1032235/1‑149 (MQ=255)
                       cccGGCGTACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGaaaa                                                                                                               <  3:100800/149‑1 (MQ=255)
                        ccGGCGTACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGaaaaa                                                                                                              <  2:1278141/149‑1 (MQ=255)
                            cGTACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTcg                                                                                                          <  3:1673079/149‑1 (MQ=255)
                              tACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAATAAGTcgcc                                                                                                        >  2:1618828/1‑149 (MQ=255)
                               aCAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTcgccg                                                                                                       <  3:1549696/149‑1 (MQ=255)
                                 aGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGt                                                                                                            >  2:162258/1‑142 (MQ=255)
                                 aGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGt                                                                                                            <  1:162258/142‑1 (MQ=255)
                                  gggAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGaaa                                                                                                                                  <  4:839780/119‑1 (MQ=255)
                                  gggAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGaaa                                                                                                                                  >  3:839780/1‑119 (MQ=255)
                                   ggAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGaa                                                                                                   <  2:1767877/149‑1 (MQ=255)
                                      gTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGaaa                                                                                                                                  <  4:1317611/115‑1 (MQ=255)
                                      gTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGaaa                                                                                                                                  >  3:1317611/1‑115 (MQ=255)
                                      gTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTa                                                                                                >  3:1052325/1‑149 (MQ=255)
                                        ggATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAAc                                                                                              <  4:1052325/149‑1 (MQ=255)
                                         gATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGaaa                                                                                                                                  >  1:935954/1‑112 (MQ=255)
                                         gATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGaaa                                                                                                                                  <  2:935954/112‑1 (MQ=255)
                                         gATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACa                                                                                             >  1:19844/1‑149 (MQ=255)
                                           tatCATTCACGAATATACTGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGa                                                                                                                                       <  3:832787/105‑1 (MQ=255)
                                           tatCATTCACGAATATACTGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGa                                                                                                                                       >  4:832787/1‑105 (MQ=255)
                                            atCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACaaaa                                                                                          <  3:1032235/149‑1 (MQ=255)
                                               aTTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCt                                                                                                                                         >  3:681539/1‑99 (MQ=255)
                                               aTTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCt                                                                                                                                         <  4:681539/99‑1 (MQ=255)
                                                 tCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAatggat                                                                                     <  1:538279/149‑1 (MQ=255)
                                                 tCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAatggat                                                                                     <  2:19844/149‑1 (MQ=255)
                                                   aCGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGTCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTcgccg                                                                                                       <  4:1394260/129‑1 (MQ=255)
                                                   aCGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGaaaaa                                                                                                              >  2:55383/1‑122 (MQ=255)
                                                   aCGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGaaaaa                                                                                                              <  1:55383/122‑1 (MQ=255)
                                                   aCGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTcgccg                                                                                                       >  3:1394260/1‑129 (MQ=255)
                                                         ataCCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCaaa                                                                             <  1:1618828/149‑1 (MQ=255)
                                                          taCCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAAcgcg                                                                                                                            >  4:35914/1‑101 (MQ=255)
                                                          taCCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAAcgcg                                                                                                                            <  3:35914/101‑1 (MQ=255)
                                                           aCCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGaaaaa                                                                                                              >  1:1472175/1‑114 (MQ=255)
                                                           aCCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGaaaaa                                                                                                              <  2:1472175/114‑1 (MQ=255)
                                                                 gCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGaaa                                                                                                                                  >  1:1705418/1‑88 (MQ=255)
                                                                 gCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGaaa                                                                                                                                  <  2:1705418/88‑1 (MQ=255)
                                                                    gcgcCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTaa                                                                                               <  3:387413/120‑1 (MQ=255)
                                                                    gcgcCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTaa                                                                                               >  4:387413/1‑120 (MQ=255)
                                                                         aGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAatggat                                                                                     >  3:429659/1‑125 (MQ=255)
                                                                         aGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAatggat                                                                                     <  4:429659/125‑1 (MQ=255)
                                                                              ctgGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGaaaaa                                                                                                              <  1:1363344/95‑1 (MQ=255)
                                                                              ctgGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGaaaaa                                                                                                              >  2:1363344/1‑95 (MQ=255)
                                                                               tgGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAg                                                                                                             <  3:1759409/95‑1 (MQ=255)
                                                                               tgGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAg                                                                                                             >  4:1759409/1‑95 (MQ=255)
                                                                               tgGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCaaa                                                                             <  4:1025070/127‑1 (MQ=255)
                                                                               tgGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCaaa                                                                             >  3:1025070/1‑127 (MQ=255)
                                                                                        ccGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGTTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAc                                                                        <  2:1790981/123‑1 (MQ=255)
                                                                                        ccGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGTTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAc                                                                        >  1:1790981/1‑123 (MQ=255)
                                                                                        ccGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGTTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAc                                                                        >  1:974984/1‑123 (MQ=255)
                                                                                        ccGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGTTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAc                                                                        <  2:974984/123‑1 (MQ=255)
                                                                                            gAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAatggatg                                                                                    >  3:973375/1‑107 (MQ=255)
                                                                                            gAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAatggatg                                                                                    <  4:973375/107‑1 (MQ=255)
                                                                                            gAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAatggatg                                                                                    >  1:1143600/1‑107 (MQ=255)
                                                                                            gAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAatggatg                                                                                    <  2:1143600/107‑1 (MQ=255)
                                                                                                   tCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAg                                                                                                                                      >  3:482285/1‑50 (MQ=255)
                                                                                                   tCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAg                                                                                                                                      <  4:482285/50‑1 (MQ=255)
                                                                                                               aGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAccatc                                                                    >  4:51949/1‑104 (MQ=255)
                                                                                                               aGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAccatc                                                                    <  3:51949/104‑1 (MQ=255)
                                                                                                                        gTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGa                     >  4:1826137/1‑142 (MQ=255)
                                                                                                                        gTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGa                     <  3:1826137/142‑1 (MQ=255)
                                                                                                                                    cTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTc  >  1:499687/1‑149 (MQ=255)
                                                                                                                                     tttCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCa                                                                         <  3:506942/77‑1 (MQ=255)
                                                                                                                                     tttCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCa                                                                         >  4:506942/1‑77 (MQ=255)
                                                                                                                                        ccTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAccatcca                                                                  >  2:258426/1‑81 (MQ=255)
                                                                                                                                        ccTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAccatcca                                                                  <  1:258426/81‑1 (MQ=255)
                                                                                                                                               gcTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGc             <  4:314471/127‑1 (MQ=255)
                                                                                                                                               gcTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGc             >  3:314471/1‑127 (MQ=255)
                                                                                                                                                |                                                                                                                                        
GTGCCAGCATTAACACTGACGATCCCGGCGTACAGGGAGTGGATATCATTCACGAATATACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTC  >  NC_000913/1703047‑1703327

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: