Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 430,702 G→C 8.1% Q539E (CAG→GAG)  PP_0354 ← CBS domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947430,7020GC 8.1% 193.4 / 8.2 62Q539E (CAG→GAG) PP_0354CBS domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (31/26);  new base C (2/3);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

GTAGGCCTCTTCATAGGCGGCGCCGTCCAACGGTTCGATCACGCCCTTGGCTACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAGGGCGTTGTCAGCCATCATGCGCTGGAAGATGCGGTTGTCGGCTACCTGTGCCAGCAG  >  NC_002947/430556‑430849
                                                                                                                                                  |                                                                                                                                                    
gTAAGCCTCTTCATAGGCGGCGCCGTCCAACGGTTCGATCACGCCCTTGGCTACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGAc                                                                                                                                                   >  2:101633/1‑150 (MQ=255)
      ctctTCATAGGCGGCGCCGTCCAACGGTTCGATCACGCCCTTGGCTACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAg                                                                                                                                             >  1:112743/1‑150 (MQ=255)
      ctctTCATAGGCGGCGCCGTCCAACGGTTCGATCACGCCCTTGGCTACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAg                                                                                                                                             >  2:847093/1‑150 (MQ=255)
      ctctTCATAGGCGGCGCCGTCCAACGGTTCGATCACGCCCTTGGCTACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAg                                                                                                                                             >  1:323301/1‑150 (MQ=255)
            aTAGGCGGCGCCGTCCAACGGTTCGATCACGCCCTTGGCTACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAg                                                                                                                                       <  1:847093/150‑1 (MQ=255)
               ggcggcGCCGTCCAACGGTTCGATCACGCCCTTGGCTACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTCGACCTCAAGGTCCAgggc                                                                                                                                    <  1:201354/150‑2 (MQ=255)
                gcggcgCCGTCCAACGGTTCGATCACGCCCTTGGCTACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTc                                                                                                                                   <  2:779990/150‑1 (MQ=255)
                        gTCCAACGGTTCGATCACGCCCTTGGCTACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCtt                                                                                                                           <  2:1498468/150‑1 (MQ=255)
                          ccAACGGTTCGATCACGCCCTTGGCTACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGt                                                                                                                         <  1:996410/150‑1 (MQ=255)
                           cAACGGTTCGATCACGCCCTTGGCTACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGa                                                                                                                                        <  1:468273/134‑1 (MQ=255)
                           cAACGGTTCGATCACGCCCTTGGCTACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGa                                                                                                                                        >  2:468273/1‑134 (MQ=255)
                            aaCGGTTCGATCACGCCCTTGGCTACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGa                                                                                                                                              >  1:592385/1‑127 (MQ=255)
                            aaCGGTTCGATCACGCCCTTGGCTACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGa                                                                                                                                              <  2:592385/127‑1 (MQ=255)
                            aaCGGTTCGATCACGCCCTTGGCTACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCg                                                                                                                       >  2:977563/1‑150 (MQ=255)
                              cGGTTCGATCACGCCCTTGGCTACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGtt                                                                                                                     <  2:112743/150‑1 (MQ=255)
                              cGGTTCGATCACGCCCTTGGCTACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGtt                                                                                                                     <  1:977563/150‑1 (MQ=255)
                                    gATCACGCCCTTGGCTACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTg                                                                                                               <  1:1201187/150‑1 (MQ=255)
                                      tCACGCCCTTGGCTACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGc                                                                                                             >  2:1288887/1‑150 (MQ=255)
                                        aCGCCCTTGGCTACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCCTTGCCTTGGCgt                                                                                                           >  1:743280/1‑150 (MQ=255)
                                            ccTTGGCTACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCCGCAAACGGGCACCATCCACGATGGGCGTCAGGTCCTGGACCTTGAGGTCGAGGGTCCCCGCCTTTTCGTTGCCTGGGCGTGTcg                                                                                                       >  2:289359/1‑149 (MQ=255)
                                              ttGGCTACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTCGACCTCAAGGTCCAggg                                                                                                                                     <  2:673213/118‑1 (MQ=255)
                                              ttGGCTACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTCGACCTCAAGGTCCAggg                                                                                                                                     >  1:673213/1‑118 (MQ=255)
                                                   tACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCgtgt                                                                                                         <  2:1537595/141‑1 (MQ=255)
                                                   tACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCgtgt                                                                                                         >  1:1537595/1‑141 (MQ=255)
                                                    aCCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAAc                                                                                               >  1:1125424/1‑150 (MQ=255)
                                                        gTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTcgc                                                                                           >  2:915061/1‑150 (MQ=255)
                                                         ttGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTcgcg                                                                                          >  1:349455/1‑150 (MQ=255)
                                                                aGGCGCTCCCAGGGAGTGCATGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCgg                                                                                   <  2:349455/150‑1 (MQ=255)
                                                                aGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTcgccgc                                                                                                                              >  2:1537581/1‑107 (MQ=255)
                                                                aGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTcgccgc                                                                                                                              <  1:1537581/107‑1 (MQ=255)
                                                                      tCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCCCGACCTTTCCCTACCGGCCCAcc                                                                             >  2:456553/1‑150 (MQ=255)
                                                                      tCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTCGACCTCAAGGTccag                                                                                                                                       >  2:614534/1‑89 (MQ=255)
                                                                      tCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTCGACCTCAAGGTCCAg                                                                                                                                       <  1:614534/92‑1 (MQ=255)
                                                                        cAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACcgg                                                                           >  2:163546/1‑150 (MQ=255)
                                                                                       agCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAg                                                           >  1:1398306/1‑150 (MQ=255)
                                                                                        gCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAgg                                                          <  2:684777/150‑1 (MQ=255)
                                                                                               cATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAGGGCGTTGt                                                   >  1:338838/1‑150 (MQ=255)
                                                                                                     ccagggccagCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAGGGCGTTGTCAGCca                                             <  2:743280/150‑1 (MQ=255)
                                                                                                     ccagggccagCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAGGGCGTTGTCAGCca                                             >  2:411719/1‑150 (MQ=255)
                                                                                                      cagggccagCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTcgcg                                                                                          <  2:1019838/105‑1 (MQ=255)
                                                                                                      cagggccagCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTAGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTcgcg                                                                                          >  1:1019838/1‑105 (MQ=255)
                                                                                                       agggccagCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAGGGCGTTGTCAGCcatc                                           <  2:54464/150‑1 (MQ=255)
                                                                                                                 aaCGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAGGGCGTTGTCAGCCATCATGCGCTGGa                                 <  2:1398306/150‑1 (MQ=255)
                                                                                                                  aCGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAGGGCGTTGTCAGCCATCATGCGCTGGaa                                <  1:411719/150‑1 (MQ=255)
                                                                                                                        aCCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAGGGCGTTGTCAGCCATCATGCGCTGGAAGATGCg                          <  2:323301/150‑1 (MQ=255)
                                                                                                                           aTCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAGGGCGTTGTCAGCCATCATGCGCTGGAAGATGCGGtt                       >  1:1324331/1‑150 (MQ=255)
                                                                                                                           aTCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAGGGCGTTGTCAGCCATCATGCGCTGGAAGATGCGGtt                       >  1:597373/1‑150 (MQ=255)
                                                                                                                            tCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGCGGGTCGCCGCCTTGTCGTTGCCTTGGCGTCACAGCGCGAACTCGCGCAACCGGCGCACCGGCGGGCGTTGGCGCAGTGCGTTGTCAGCTATCCTGCGGTGGAAGATGCGGTTg                      >  2:810412/1‑150 (MQ=255)
                                                                                                                              cACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAGGGCGTTGTCAGCCATCATGCGCTGGAAGATGCGGTTGTc                    >  2:1047947/1‑150 (MQ=255)
                                                                                                                              cACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAGGGCGTTGTCAGCCATCATGCGCTGGAAGATGCGGTTGTc                    <  1:1047947/150‑1 (MQ=255)
                                                                                                                               aCGAAGGGGGTCAGGCCCTGGACCTTGAGGTGGAGGGTCGCCGCCTGGTCGTTGCCTTGGCGGGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAGGGCGTTGTCAGCCATCATGCTCTGGAAGATGCGGTTGTCg                   >  2:918383/1‑150 (MQ=255)
                                                                                                                               aCGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAGGGCGTTGTCAGCCATCATGCGCTGGAAGATGCGGTTGTCg                   <  1:172724/150‑1 (MQ=255)
                                                                                                                                 gAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTTGCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAGGGCGTTGTCAGCCATCATGCGCTGGAAGATGCGGTTGTCGGc                 >  2:723899/1‑150 (MQ=255)
                                                                                                                                      gCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTgg                                                                >  2:1110719/1‑98 (MQ=255)
                                                                                                                                      gCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTgg                                                                <  1:1110719/98‑1 (MQ=255)
                                                                                                                                      gCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAGGGCGTTGTCAGCCATCATGCGCTGGAAGATGCGGTTGTCGGCTACCt            >  2:171732/1‑150 (MQ=255)
                                                                                                                                          cAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAGGGCGTTGTCAGCCATCATGCGCTGGAAGATGCGGTTGTCGGCTACCTGTGc        <  1:915061/150‑1 (MQ=255)
                                                                                                                                          cAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAGGGCGTTGTCAGCCATCATGCGCTGGAAGATGCGGTTGTCGGCTACCTGTGc        <  1:1288887/150‑1 (MQ=255)
                                                                                                                                           aGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGGCGCGAAACGGCTGGGCCATCATGCGCTGGAAGATGCGGTTGTCGGCTACCTGTGcc       >  2:861846/1‑150 (MQ=255)
                                                                                                                                            ggCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAGGGCGTTGTCAGCCATCATGCGCTGGAAGATGCGGTTGTCGGCTACCTGTGCca      >  1:1240184/1‑150 (MQ=255)
                                                                                                                                                cTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAGGGCGTTGTCAGCCATCATGCGCTGGAAGATGCGGTTGTCGGCTACCTGTGCcagcag  <  2:1125424/150‑1 (MQ=255)
                                                                                                                                                cTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAGGGCGTTGGCGTCCATCATGCGCTGGAAGATGCGGTTGTCGGCTACCTGTGcccacag  <  1:456553/150‑6 (MQ=255)
                                                                                                                                                  ggACCTTGAGGTCGAGCGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAGGGCGTTGTCAGCCATCATGCGCTGGAAgc                              >  1:415392/1‑119 (MQ=255)
                                                                                                                                                  ggACCTTGAGGTCGAGCGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAGGGCGTTGTCAGCCATCATGCGCTGGAAGa                              <  2:415392/120‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GTAGGCCTCTTCATAGGCGGCGCCGTCCAACGGTTCGATCACGCCCTTGGCTACCAGTTGGCGCAGGCGCTCCAGGGTATTGCAGGCAGCTATGCCATTGGCCAGGGCCAGCAAACGGGCACCATCCACGAAGGGCGTCAGGCCCTGGACCTTGAGGTCGAGGGTCGCCGCCTTGTCGTTGCCTTGGCGTGTCAGCACGAACTCGCGCAACCGGCCCACCGGCGGGCGTTGGCGCAGGGCGTTGTCAGCCATCATGCGCTGGAAGATGCGGTTGTCGGCTACCTGTGCCAGCAG  >  NC_002947/430556‑430849

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