Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,965,577 A→T 6.9% M128L (ATG→TTG)  PP_3496 → endoribonuclease L‑PSP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,965,5770AT 6.9% 167.9 / 2.9 58M128L (ATG→TTG) PP_3496endoribonuclease L‑PSP
Reads supporting (aligned to +/- strand):  ref base A (26/28);  new base T (2/2);  total (28/30)
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.34e-01

CACCTGACCATCTACCTGACCGACATGGCCGACCGCGCAGCGTTCAACGAGGTGTACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCG  >  NC_002947/3965430‑3965721
                                                                                                                                                    |                                                                                                                                                 
cACCTGACCATCTACCTGACCGACATGGCCGACCGCGCAGCGTTCAACGAGGTGTACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGa                                                                                                                                                   <  6:60756/148‑1 (MQ=255)
 aCCTGACCATCTACCTGACCGACATGGCCGACCGCGCAGCGTTCAACGAGGTGTACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATg                                                                                                                                                 <  3:118933/149‑1 (MQ=255)
 aCCTGACCATCTACCTGACCGACATGGCCGACCGCGCAGCGTTCAACGAGGTGTACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATg                                                                                                                                                 >  4:360275/1‑149 (MQ=255)
      aCCATCTACCTGACCGACATGGCCGACCGCGCAGCGTTCAACGAGGTGTACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCt                                                                                                                                              <  1:186710/147‑1 (MQ=255)
       ccATCTACCTGACCGACATGGCCGACCGCGCAGCGTTCAACGAGGTGTACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGcc                                                                                                                                           <  7:304817/149‑1 (MQ=255)
        cATCTACCTGACCGACATGGCCGACCGCGCAGCGTTCAACGAGGTGTACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGTTGGCTGcc                                                                                                                                           >  6:151539/1‑148 (MQ=255)
        cATCTACCTGACCGACATGGCCGACCGCGCAGCGTTCAACGAGGTGTACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGcc                                                                                                                                           >  4:287492/1‑148 (MQ=255)
          tCTACCTGACCGACATGGCCGACCGCGCAGCGTTCAACGAGGTGTACCCGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAg                                                                                                                                        >  6:295135/1‑149 (MQ=255)
           cTACCTGACCGACATGGCCGACCGCGCAGCGTTCAGCGAGGTGTACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAgc                                                                                                                                       >  4:93715/1‑149 (MQ=255)
                  aCCGACATGGCCGACCGCGCAGCGTTCAGCGAGGTGTACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGa                                                                                                                                  <  3:93715/147‑1 (MQ=255)
                  aCCGACATGGCCGACCGCGCAGCGTTCAACGAGGTGTACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGa                                                                                                                                  <  4:351306/147‑1 (MQ=255)
                             cGACCGCGCAGCGTTCAACGAGGTGTACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTgg                                                                                                                     >  5:175068/1‑149 (MQ=255)
                             cGACCGCGCAGCGTTCAACGAGGTGTACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTgg                                                                                                                     <  4:113439/149‑1 (MQ=255)
                             cGACCGCGCAGCGTTCAACGAGGTGTACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTgg                                                                                                                     <  5:275265/149‑1 (MQ=255)
                              gACCGCGCAGCGTTCAACGAGGTGTACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTggg                                                                                                                    >  3:142744/1‑149 (MQ=255)
                               aCCGCGCAGCGTTCAACGAGGTGTACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTgggg                                                                                                                   <  3:360275/149‑1 (MQ=255)
                                  gcgcAGCGTTCAACGAGGTGTACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTg                                                                                                                >  1:189518/1‑149 (MQ=255)
                                    gcAGCGTTCAACGAGGTGTACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTgcg                                                                                                              >  3:309834/1‑149 (MQ=255)
                                      aGCGTTCAACGAGGTTTACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCa                                                                                                            <  1:351488/149‑1 (MQ=255)
                                       gCGTTCAACGAGGTGTACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAg                                                                                                           <  2:189518/149‑1 (MQ=255)
                                          ttCAACGAGGTGTACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAgcg                                                                                                                                      <  4:133756/119‑1 (MQ=255)
                                          ttCAACGAGGTGTACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAgcg                                                                                                                                      >  3:133756/1‑119 (MQ=255)
                                                  ggTGTACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAt                                                                                                >  1:379587/1‑149 (MQ=255)
                                                    tgtACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGTTGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑Cg                                                                                             >  1:322202/1‑149 (MQ=255)
                                                       aCCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTtgt                                                                                                                       >  7:222179/1‑121 (MQ=255)
                                                       aCCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTtgt                                                                                                                       <  8:222179/121‑1 (MQ=255)
                                                        ccAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTTCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCg                                                                                                    <  6:208061/139‑1 (MQ=255)
                                                        ccAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCg                                                                                                    >  5:208061/1‑139 (MQ=255)
                                                        ccAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGa                                                                                         >  3:389756/1‑149 (MQ=255)
                                                        ccAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGa                                                                                         <  3:130133/149‑1 (MQ=255)
                                                         cAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGAt                                                                                        >  1:189566/1‑149 (MQ=255)
                                                         cAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGAt                                                                                        >  3:338861/1‑149 (MQ=255)
                                                         cAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGAt                                                                                        <  6:15302/149‑1 (MQ=255)
                                                             gcTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGATTGGt                                                                                    >  1:75690/1‑149 (MQ=255)
                                                              cttcttCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGATTGGTc                                                                                   <  4:175069/149‑1 (MQ=255)
                                                                 cttcAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAATTCCGGGGGATGGATCCAAACCGATGAGCT‑CGGATGTTGGGGTTCCCCACCCCCCCGCATACGCTGATTggagga                                                                                 >  6:15883/1‑144 (MQ=255)
                                                                 cttcAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGTTGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGATTGGTCGc                                                                                 <  5:151539/148‑1 (MQ=255)
                                                                     aGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGATTGGTCGCAAAAt                                                                            >  2:218402/1‑149 (MQ=255)
                                                                      gCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGATTGGTCGCAAAATc                                                                           <  2:379587/149‑1 (MQ=255)
                                                                          gCCCTGGCCGGTGCGTGCGGCGGTGTGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTTGG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGTGTGCGCAGCGCCCCGACAT‑CGCTGATTGGTCGCAAAATCGCCg                                                                       <  5:52305/149‑1 (MQ=255)
                                                                            ccTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGATTGGTCGCAAAATCGCCggg                                                                     >  7:332105/1‑149 (MQ=255)
                                                                                 ccGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGATTGGTCGCAAAATCGCCGGGGCCGc                                                                >  1:339873/1‑149 (MQ=255)
                                                                                         tgcgGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGTTGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGc                                                        <  2:322202/149‑1 (MQ=255)
                                                                                          gcggcgGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGcc                                                       >  2:333648/1‑149 (MQ=255)
                                                                                           cggcggTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTGAGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGcc                                                       >  1:81368/1‑149 (MQ=255)
                                                                                           cggcggTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGcc                                                       <  8:233410/148‑1 (MQ=255)
                                                                                            ggcggTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGGTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCt                                                     <  3:287492/149‑1 (MQ=255)
                                                                                                  gggtgTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGATGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTcgcg                                               <  8:10176/149‑1 (MQ=255)
                                                                                                     tgtgGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACg                                            >  7:4319/1‑149 (MQ=255)
                                                                                                     tgtgGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACg                                            <  8:332105/149‑1 (MQ=255)
                                                                                                              aCTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGc                                   <  4:338861/149‑1 (MQ=255)
                                                                                                               cTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCt                                  <  2:75690/149‑1 (MQ=255)
                                                                                                                     tttGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCca                              >  6:223690/1‑147 (MQ=255)
                                                                                                                         aaGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGtt                         <  3:235938/148‑1 (MQ=255)
                                                                                                                         aaGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGtt                         <  8:4319/148‑1 (MQ=255)
                                                                                                                          aGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGGGGGGGTGCGCAGCGCCCCGACAT‑CGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTgc                       <  4:309834/149‑1 (MQ=255)
                                                                                                                            gCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTgcgc                     >  3:381244/1‑149 (MQ=255)
                                                                                                                                            tGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAAt      <  4:389756/149‑1 (MQ=255)
                                                                                                                                                ggCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCGGCGCCCCGACAT‑CGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCg  >  3:303477/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CACCTGACCATCTACCTGACCGACATGGCCGACCGCGCAGCGTTCAACGAGGTGTACCAGCGCTTCTTCAGCAAGCCCTGGCCGGTGCGTGCGGCGGTGGGTGTGGCTTCACTGGCCTTTGAAGGCATGCGTGTG‑GAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACAT‑CGCTGATTGGTCGCAAAATCGCCGGGGCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCG  >  NC_002947/3965430‑3965721

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