Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,873,567 G→A 5.4% R142R (CGC→CGT cysC ← adenosine 5'‑phosphosulfate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,873,5670GA 5.4% 241.2 / 2.3 74R142R (CGC→CGTcysCadenosine 5'‑phosphosulfate kinase
Reads supporting (aligned to +/- strand):  ref base G (35/35);  new base A (2/2);  total (37/37)
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.62e-01

TAATAATTGCTGTACCAAATTTGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTC  >  NC_000913/2873423‑2873714
                                                                                                                                                |                                                                                                                                                   
taataaTTGCTGTACCAAATTTGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCACggg                                                                                                                                                 >  1:266100/1‑149 (MQ=255)
taataaTTGCTGTACCAAATTTGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCACggg                                                                                                                                                 >  1:271517/1‑149 (MQ=255)
  ataatTGCTGTACCAAATTTGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATcgc                                                                                                                                                    <  4:56601/144‑1 (MQ=255)
  ataatTGCTGTACCAAATTTGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCgg                                                                                                                                                  <  2:237332/146‑1 (MQ=255)
  ataatTGCTGTACCAAATTTGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCgg                                                                                                                                                  <  1:85646/146‑1 (MQ=255)
    aatTGCTGTACCAAATTTGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTc                                                                                                                                             >  1:30198/1‑149 (MQ=255)
     atTGCTGTACCAAATTTGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCg                                                                                                                                            <  2:177317/149‑1 (MQ=255)
     atTGCTGTACCAAATTTGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCATAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCGTc                                                                                                                                             <  5:54958/148‑1 (MQ=255)
            tACCAAATTTGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCg                                                                                                                                     >  1:88879/1‑149 (MQ=255)
            tACCAAATTTGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTACCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCg                                                                                                                                     >  1:297818/1‑149 (MQ=255)
              ccAAATTTGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGc                                                                                                                                    <  3:190282/148‑1 (MQ=255)
              ccAAATATGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATGAGCCTTTGGGATCGCGGGCTTCGCAAATCGc                                                                                                                                    >  3:352208/1‑148 (MQ=255)
                 aaTTTGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAgc                                                                                                                                >  7:257803/1‑149 (MQ=255)
                 aaTTTGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAgc                                                                                                                                >  6:311031/1‑149 (MQ=255)
                     tGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAgcggc                                                                                                                             <  5:311031/148‑1 (MQ=255)
                          cTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTc                                                                                                                                             <  7:278587/127‑1 (MQ=255)
                          cTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTc                                                                                                                                             >  8:278587/1‑127 (MQ=255)
                          cTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATc                                                                                                                        <  1:304674/148‑1 (MQ=255)
                           tAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAgcggcg                                                                                                                            <  7:157517/143‑1 (MQ=255)
                           tAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAgcggcg                                                                                                                            >  8:157517/1‑143 (MQ=255)
                            aaTTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGc                                                                                                                                           >  7:328106/1‑127 (MQ=255)
                            aaTTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGATTCGc                                                                                                                                           <  8:328106/127‑1 (MQ=255)
                             aTTGTTCAGCATTGAGAAGAATTTCTGCCGATTCAGGCGCTTCGTATACTGAATCTATTCCCGTATAGTTGCGCAGTTCACCGGCAGGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGAc                                                                                                                     <  4:105478/148‑1 (MQ=255)
                             aTTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCACGGGCTTCGCAAATCGCCAGCGGCGTATCGACa                                                                                                                    <  2:266100/149‑1 (MQ=255)
                             aTTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCACGGGCTTCGCAAATCGCCAGCGGCGTATCGACa                                                                                                                    <  2:271517/149‑1 (MQ=255)
                              ttGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACa                                                                                                                    >  8:139425/1‑148 (MQ=255)
                                 ttCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAaca                                                                                                                >  3:148026/1‑149 (MQ=255)
                                 ttCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAac                                                                                                                 <  5:17853/148‑1 (MQ=255)
                                  tCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAacac                                                                                                               >  6:127241/1‑149 (MQ=255)
                                    aCCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACtt                                                                                                             <  3:166346/149‑1 (MQ=255)
                                      cATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTc                                                                                                            >  2:221386/1‑148 (MQ=255)
                                           agaTGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATa                                                                                                        >  7:1822/1‑147 (MQ=255)
                                            gaTGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATc                                                                                                                                      >  1:152931/1‑116 (MQ=255)
                                            gaTGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATc                                                                                                                                      <  2:152931/116‑1 (MQ=255)
                                               gAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACtt                                                                                                             >  3:174823/1‑138 (MQ=255)
                                               gAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACtt                                                                                                             <  4:174823/138‑1 (MQ=255)
                                                         cGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTa                                                                                        <  6:303372/149‑1 (MQ=255)
                                                               aGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAgcg                                                                                                                               <  5:140983/104‑1 (MQ=255)
                                                               aGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAgcg                                                                                                                               >  6:140983/1‑104 (MQ=255)
                                                               aGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGtt                                                                                  >  6:247860/1‑149 (MQ=255)
                                                               aGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGt                                                                                   <  1:221386/148‑1 (MQ=255)
                                                                ggCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTc                                                                                 <  2:30198/149‑1 (MQ=255)
                                                                  cgcTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTcgc                                                                               <  1:7525/149‑1 (MQ=255)
                                                                     ttCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGaa                                                                            >  4:329332/1‑149 (MQ=255)
                                                                       cGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGaa                                                                            <  8:28438/147‑1 (MQ=255)
                                                                            aCGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGAACCATCTg                                                                     <  2:88879/149‑1 (MQ=255)
                                                                                 aTCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGAACCATCTGGCGtt                                                                >  4:80709/1‑149 (MQ=255)
                                                                                      ttCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAacac                                                                                                               >  4:188004/1‑97 (MQ=255)
                                                                                      ttCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAacac                                                                                                               <  3:188004/97‑1 (MQ=255)
                                                                                        cccGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGAACCATCTGGCGTTCGgcgc                                                          <  2:297818/148‑1 (MQ=255)
                                                                                         ccGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATc                                                                                                                        <  4:301616/85‑1 (MQ=255)
                                                                                         ccGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATc                                                                                                                        >  3:301616/1‑85 (MQ=255)
                                                                                         ccGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCgg                                                        >  7:48580/1‑149 (MQ=255)
                                                                                            tAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGgtgt                                                     >  3:328482/1‑149 (MQ=255)
                                                                                               aGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATGAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGc                                                  <  4:352208/149‑1 (MQ=255)
                                                                                                gTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGAt                                                                                                         >  2:114415/1‑93 (MQ=255)
                                                                                                gTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGAt                                                                                                         <  1:114415/93‑1 (MQ=255)
                                                                                                 ttGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGc                                                  >  5:68772/1‑147 (MQ=255)
                                                                                                    cgcAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCgaga                                              >  1:176418/1‑148 (MQ=255)
                                                                                                        gTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAAt                                         <  3:329332/149‑1 (MQ=255)
                                                                                                          tCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTAGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGc                                       >  1:361191/1‑149 (MQ=255)
                                                                                                            aCCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCgg                                     >  7:36597/1‑149 (MQ=255)
                                                                                                            aCCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCgg                                     >  3:78973/1‑149 (MQ=255)
                                                                                                            aCCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCgg                                     >  5:3721/1‑149 (MQ=255)
                                                                                                                 cACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGc                                <  2:176418/149‑1 (MQ=255)
                                                                                                                  aCGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGca                               <  4:328482/149‑1 (MQ=255)
                                                                                                                  aCGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGca                               <  4:148026/149‑1 (MQ=255)
                                                                                                                   cgcgCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGcac                              <  3:80709/149‑1 (MQ=255)
                                                                                                                     cgcTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGcacca                            <  8:1822/149‑1 (MQ=255)
                                                                                                                             tatataAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTtctc                                                                                           <  7:210760/78‑1 (MQ=255)
                                                                                                                             tatataAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTtctc                                                                                           >  8:210760/1‑78 (MQ=255)
                                                                                                                                ataAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGAACCATCTg                                                                     <  6:109617/97‑1 (MQ=255)
                                                                                                                                ataAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGAACCATCTg                                                                     >  5:109617/1‑97 (MQ=255)
                                                                                                                                    gCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCaa             >  6:213905/1‑149 (MQ=255)
                                                                                                                                    gCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCaa             >  4:10913/1‑149 (MQ=255)
                                                                                                                                                gcgGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTc  <  3:10913/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TAATAATTGCTGTACCAAATTTGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTTTGGGATCGCGGGCTTCGCAAATCGCCAGCGGCGTATCGACAAACACTTCGATAAAGCGCCCTTCTCCTACGCGTTCGCGAACCATCTGGCGTTCGGCGCGGTGTGGCGAGATAAATGCGGTCAGCACCACCAGTCCGGCTTCAACCATCAAATTC  >  NC_000913/2873423‑2873714

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