Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,579,072 T→A 6.0% F142L (TTT→TTA)  PP_5589 → rhs repeat‑associated core domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,579,0720TA 6.0% 190.3 / 2.6 67F142L (TTT→TTA) PP_5589rhs repeat‑associated core domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (29/34);  new base A (2/2);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

TACCGCCGCAGGATATAATTCAGTCCAGCCCTGAGTTAAGGCCTTTTGTAGAACGTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAGCGCTAGAATTCCTGAGTTGATGAAAGTTACTGCTTGGATTAGGGCTGAACG  >  NC_002947/3578927‑3579212
                                                                                                                                                 |                                                                                                                                            
tACCGCCGCAGGATATAATTCAGTCCAGCCCTGAGTTAAGGCCTTTTGTAGAACGTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATc                                                                                                                                           <  7:248224/149‑1 (MQ=255)
    gccgcAGGATATAATTCAGTCCAGCCCTGAGTTAAGGCCTTTTGTAGAACGTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGc                                                                                                                                       <  3:313147/149‑1 (MQ=255)
      cgcAGGATATAATTCAGTCCAGCCCTGAGTTAAGGCCTTTTGTAGAACGTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATAAAATCAAGcac                                                                                                                                     <  4:139265/149‑3 (MQ=255)
      cgcAGGATATAATTCAGTCCAGCCCTGAGTTAAGGCCTTTTGTAGAACGTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATAAAATCAAGcac                                                                                                                                     <  5:165342/149‑3 (MQ=255)
         aGGATATAATTCAGTCCAGCCCTGAGTTAAGGCCTTTTGTAGAACGTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGtgtt                                                                                                                                  <  4:137870/149‑1 (MQ=255)
               taATTCAGTCCAGCCCTGAGTTAAGGCCTTTTGTGGAACGTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGc                                                                                                                             <  8:74737/148‑1 (MQ=255)
                aaTTCAGTCCAGCCCTGAGTTAAGGCCTTTTGTAGAACGTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGccc                                                                                                                           >  2:185225/1‑149 (MQ=255)
                 aTTCAGTCCAGCCCTGAGTTAAGGCCTTTTGTAGAACGTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGcccc                                                                                                                          >  2:37711/1‑149 (MQ=255)
                 aTTCAGTCCAGCCCTGAGTTAAGGCCTTTTGTAGAACGTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGccc                                                                                                                           <  7:135339/148‑1 (MQ=255)
                 aTTCAGTCCAGCCCTGAGTTAAGGCCTTTTGTAGAACGTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGccc                                                                                                                           <  7:322134/148‑1 (MQ=255)
                       tCCAGCCCTGAGTTAAGGCCTTTTGTAGAACGTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGt                                                                                                                    >  3:147269/1‑149 (MQ=255)
                         cAGCCCTGAGTTAAGGCCTTTTGTAGAACGTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTc                                                                                                                   <  1:37711/148‑1 (MQ=255)
                             ccTGAGTTAAGGCCTTTTGTAGAACGTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACgg                                                                                                              <  1:185225/149‑1 (MQ=255)
                                 aGTTAAGGCCTTTTGTAGAACGTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACgggg                                                                                                            <  8:204035/147‑1 (MQ=255)
                                  ggtAAGGCCTTTTGTAGAACGTAATTCGAAGTTGCCTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTCGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGggtg                                                                                                         <  6:255629/147‑1 (MQ=255)
                                  gTTAAGGCCTTTTGTAGAACGTAATTCGAAGTTGACTGTTCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGggtg                                                                                                         <  7:163747/149‑1 (MQ=255)
                                  gTTAAGGCCTTTTGTAGAACGTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGggtg                                                                                                         <  5:111502/149‑1 (MQ=255)
                                   ttAAGGCCTTTTGTAGAACGTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGggtgg                                                                                                        >  2:12843/1‑149 (MQ=255)
                                     aaGGCCTTTTGTAGAACGTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGgtggtg                                                                                                      >  4:273644/1‑149 (MQ=255)
                                      aGGCCTTTTGTAGAACGTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGgtgt                                                                                                     <  8:124875/149‑1 (MQ=255)
                                           ttttGTAGAACGTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGt                                                                                                                    <  6:223096/129‑1 (MQ=255)
                                           ttttGTAGAACGTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGt                                                                                                                    >  5:223096/1‑129 (MQ=255)
                                             ttGTAGAACGTAATTCGAAGTTGAATGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAgatga                                                                                               <  7:218707/148‑1 (MQ=255)
                                               gTAGAACGTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAgatga                                                                                               >  8:218707/1‑146 (MQ=255)
                                                tAGAACGTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGAttgt                                                                                           >  7:83011/1‑149 (MQ=255)
                                                      gTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGGTGATTgtttgttt                                                                                     >  4:288624/1‑149 (MQ=255)
                                                         aTTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTtgct                                                                                  >  7:303251/1‑149 (MQ=255)
                                                            cGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTtgctgt                                                                                <  6:349102/148‑2 (MQ=255)
                                                                gTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTAt                                                                            >  7:140641/1‑148 (MQ=255)
                                                                         gcttgctCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAat                                                                   <  1:12843/148‑1 (MQ=255)
                                                                           ttgctCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAat                                                                   >  5:135388/1‑146 (MQ=255)
                                                                           ttgctCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAat                                                                   <  6:135388/146‑1 (MQ=255)
                                                                             gctCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATAtt                                                                <  4:147269/147‑1 (MQ=255)
                                                                               tCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTTATATTGCct                                                            >  8:162293/1‑149 (MQ=255)
                                                                               tCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCct                                                            >  3:47813/1‑149 (MQ=255)
                                                                                 gCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCctct                                                          >  7:170432/1‑149 (MQ=255)
                                                                                        ggggTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAgcg                                                   >  6:347736/1‑149 (MQ=255)
                                                                                           gTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAgcgc                                                  >  6:56596/1‑147 (MQ=255)
                                                                                               ccTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGgtggtg                                                                                                      <  2:186159/91‑1 (MQ=255)
                                                                                               ccTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGgtggtg                                                                                                      >  1:186159/1‑91 (MQ=255)
                                                                                               ccTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGc                                                                          >  5:287617/1‑119 (MQ=255)
                                                                                               ccTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGc                                                                          <  6:287617/119‑1 (MQ=255)
                                                                                                    aTCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAGCGCTAGAATTCCt                                        <  4:47813/148‑1 (MQ=255)
                                                                                                     tCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAGCGCTAGAATTCCTGa                                      >  3:20782/1‑149 (MQ=255)
                                                                                                       tAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAGCGCTAGAATTCCTGAGt                                    <  8:83011/149‑1 (MQ=255)
                                                                                                         gCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGGTGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAGCGCTAGAATTCCTGa                                      <  3:288624/145‑1 (MQ=255)
                                                                                                          cAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAGCGCTAGAATTCCTGAGTtga                                 >  6:281651/1‑149 (MQ=255)
                                                                                                            gAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAGCGCTAGAATTCCTGAGTtgat                                >  6:91709/1‑148 (MQ=255)
                                                                                                                  tatTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAGCGCTAGAATTCCTGAGTTGATGAAAGt                          >  4:276585/1‑148 (MQ=255)
                                                                                                                       ccGTAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAGCGCTAGAATTCCTGAGTTGATGAAAGTTACTGc                    <  8:140641/149‑1 (MQ=255)
                                                                                                                        cGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTa                                                                             >  5:255660/1‑91 (MQ=255)
                                                                                                                        cGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTa                                                                             <  6:255660/91‑1 (MQ=255)
                                                                                                                          gaggaTCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAGCGCTAGAATTCCTTAGTTGATGAAATTTACTGCTTg                 >  3:96452/1‑149 (MQ=255)
                                                                                                                          gaggaTCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAGCGCTAGAATTCCTGa                                      <  4:83574/128‑1 (MQ=255)
                                                                                                                          gaggaTCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAGCGCTAGAATTCCTGa                                      >  3:83574/1‑128 (MQ=255)
                                                                                                                             gaTCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAGCGCTAGAATTCCTGAGTTGATGaaa                            >  5:24752/1‑135 (MQ=255)
                                                                                                                             gaTCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAGCGCTAGAATTCCTGAGTTGATGaaa                            <  6:24752/135‑1 (MQ=255)
                                                                                                                              aTCTGGAAGTGGTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAGCGCTAGAATTCCTGAGTTGATGAAAGTTACTGCTTGGAt              <  8:170432/148‑1 (MQ=255)
                                                                                                                              aTCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAGCGCTAGAATTCCTGAGTTGATGAAAGTTACTGCTTGGAt              <  3:276585/148‑1 (MQ=255)
                                                                                                                                   gAAGTGCTTGATTTTATCATGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAGCGCTAGAATTCCTGAGTTGATGAAAGTTACTGCTTGGATTAgg          <  4:20782/147‑1 (MQ=255)
                                                                                                                                   gAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAGCGCTAGAATTCCTGAGTTGATGAAAGTTACTGCTTGGATTAGGGc        <  7:162293/149‑1 (MQ=255)
                                                                                                                                    aaGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAGCGCTAGAATTCCTGAGTTGATGAAAGTTACTGCTTGGATTAGGGCt       <  5:56596/149‑1 (MQ=255)
                                                                                                                                      acgCTTGATAAAATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAGCGCTAGAATTCCTGAGTTGATGAAAGTTACTGCTTGGATTAGGGCTg      >  7:289164/3‑148 (MQ=255)
                                                                                                                                        gCTTGATAAAATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAGCGCTAGATTTCCTGAGTTGATGAAAGTTTCTGCGTGGATTAGGGCTGAAc   >  4:11364/1‑149 (MQ=255)
                                                                                                                                         cTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAgcgc                                                  >  7:115743/1‑101 (MQ=255)
                                                                                                                                         cTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAgcgc                                                  <  8:115743/101‑1 (MQ=255)
                                                                                                                                         cTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAGCGCTAGAATTCCTGAGTTGATGAAAGTTACTGCTTGGATTAGGGCTGAACg  >  2:183610/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
TACCGCCGCAGGATATAATTCAGTCCAGCCCTGAGTTAAGGCCTTTTGTAGAACGTAATTCGAAGTTGACTGTGCTTGCTCGCTTGAGGGGGTCGCCTTCATCTAGCAGAAATATATTACCGGAGGATCTGGAAGTGCTTGATTTTATCAAGCGTGTTGTCGCCCCTCCTGTCGACGGGGGTGGTGTAGATGATTGTTTGTTTGCTGCTTATGCTAGTAATATTGCCTCTGGCAGCGCTAGAATTCCTGAGTTGATGAAAGTTACTGCTTGGATTAGGGCTGAACG  >  NC_002947/3578927‑3579212

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