Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,488,611 T→C 7.4% D40G (GAC→GGC) ‡ PP_3094 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,488,6110TC 7.4% 186.0 / 3.0 66D40G (GAC→GGC) ‡PP_3094hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/30);  new base C (2/3);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCATCGGTAGCCGCGAAGCGCTCGAGGTCGGTCAGCTGATACACGCGCATCACCACCGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGATGTGCTTGTGTTGCATGACGGCAAAACCTCCGGGTACTGGTCAG  >  NC_002947/3488469‑3488757
                                                                                                                                              |                                                                                                                                                  
ggCATCGGTAGCCGCGAAGCGCTCGAGGTCGGTCAGCTGATACACGCGCATCACCACCGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGgccg                                                                                                                                              >  2:232293/1‑149 (MQ=255)
ggCATCGGTAGCCGCGAAGCGCTCGAGGTCGGTCAGCTGATACACGCGCATCACCACCGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGgccg                                                                                                                                              >  2:104824/1‑149 (MQ=255)
ggCATCGGTAGCCGCGAAGCGCTCGAGGTCGGTCAGCTGATACACGCGCATCACCACCGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGgc                                                                                                                                                <  3:200975/147‑1 (MQ=255)
  cATCGGTAGCCGCGAAGCGCTCGAGGTCGGTCAGCTGATACACGCGCATCACCACCGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGgccgcc                                                                                                                                            <  4:85134/149‑1 (MQ=255)
   aTCGGTAGCCGCGAAGCGCTCGAGGTCGGTCAGCTGATACACGCGCATCACCACCGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGgccgcc                                                                                                                                            >  1:172297/1‑148 (MQ=255)
       gTAGCCGCGAAGCGCTCGAGGTCGGTCAGCTGATACACGCGCATCACCACCGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCt                                                                                                                                       <  1:232293/149‑1 (MQ=255)
            cgcgAAGCGCTCGAGGTCGGTCAGCTGATACACGCGCATCACCACCGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCGACATGCTGgccgaca                                                                                                                                           <  2:28972/140‑4 (MQ=255)
            cgcgAAGCGCTCGAGGTCGGTCAGCTGATACACGCGCATCACCACCGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCGACATGCTGgccgaca                                                                                                                                           >  1:28972/1‑137 (MQ=255)
            cgcgAAGCGCTCGAGGTCGGTCAGCTGATACACGCGCATCACCACCGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCtt                                                                                                                                  <  3:215032/149‑1 (MQ=255)
             gcgAAGCGCTCGAGGTCGGTCAGCTGATACACGCGCATCACCACCGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTc                                                                                                                                 >  8:197190/1‑149 (MQ=255)
              cgAAGCGCTCGAGGTCGGTCAGCTGATACACGCGCATCACCACCGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTc                                                                                                                                 >  1:332461/1‑148 (MQ=255)
                aaGCGCTCGAGGTCGGTCAGCTGATACACGCGCATCACCACCGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCAt                                                                                                                              <  8:87844/149‑1 (MQ=255)
                  gcgcTCGAGGGCGGTCAGCTGATACACGCGCATCACCACCGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATgc                                                                                                                            <  3:35100/149‑1 (MQ=255)
                               gtcAGCTGATACACGCGCATCACCACCGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATAcg                                                                                                               >  3:154174/1‑149 (MQ=255)
                                     tGATACACGCGCATCACCACCGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATAcggtcggt                                                                                                         >  4:100019/1‑149 (MQ=255)
                                       aTACACGCGCATCACCACCGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATAcggtcggt                                                                                                         <  5:42467/147‑1 (MQ=255)
                                          cacGCGCATCACCACCGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGAcc                                                                                                     <  1:104824/148‑1 (MQ=255)
                                           acGCGCATCACCACCGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCGACATGCTGgccgaca                                                                                                                                           <  8:60482/109‑4 (MQ=255)
                                           acGCGCATCACCACCGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCGACATGCTGgccgaca                                                                                                                                           >  7:60482/1‑106 (MQ=255)
                                           acGCGCATCACCACCGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGAcccc                                                                                                   >  8:118519/1‑149 (MQ=255)
                                           acGCGCATCACCACCGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGAcccc                                                                                                   >  3:152416/1‑149 (MQ=255)
                                                    accaccGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCgg                                                                                           <  2:332461/148‑1 (MQ=255)
                                                    accaccGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGc                                                                                          >  2:367823/1‑149 (MQ=255)
                                                    accaccGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACAccgccgccgcc                                                                                          >  8:272588/1‑147 (MQ=255)
                                                       cccGACAGCGGCATGCCCTCGGCGTCGGGGGTTTTCTGGTTTCCGCCGTCGGTCGGGATGATTTCCTTTTCGCTGTGGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAg                                                                                       <  4:108652/148‑1 (MQ=255)
                                                          gACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATAcggtcggt                                                                                                         >  2:21492/1‑128 (MQ=255)
                                                          gACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATAcggtcggt                                                                                                         <  1:21492/128‑1 (MQ=255)
                                                           aCAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCAt                                                                                                                              <  6:351074/106‑1 (MQ=255)
                                                           aCAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCAt                                                                                                                              >  5:351074/1‑106 (MQ=255)
                                                                         tCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCCGCAAggcggc                                                                     >  8:354048/1‑149 (MQ=255)
                                                                            gCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGc                                                                  <  1:367823/149‑1 (MQ=255)
                                                                                 ggggTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGt                                                               >  1:179331/1‑147 (MQ=255)
                                                                                 ggggTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCa                                                             <  6:311322/149‑1 (MQ=255)
                                                                                  gggTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGTCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAg                                                            >  5:108185/1‑149 (MQ=255)
                                                                                            ggTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAg                                                   >  3:201131/1‑148 (MQ=255)
                                                                                             gTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGa                                                 <  7:118519/149‑1 (MQ=255)
                                                                                             gTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGa                                                 >  8:213560/1‑149 (MQ=255)
                                                                                             gTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGa                                                 <  5:109890/149‑1 (MQ=255)
                                                                                                gcccccGGCGGAGGTCAGCCTGACCTTTGCGCGGGCCGCCAGCCTGCCGGCCGCCCAGCTGTCTTCCCTGCGCTTGCCGCTACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGa                                              <  4:255279/149‑1 (MQ=255)
                                                                                                  cgccgTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGttt                                                      <  2:290803/139‑1 (MQ=255)
                                                                                                  cgccgTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGttt                                                      >  1:290803/1‑139 (MQ=255)
                                                                                                  cgccgTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGAtg                                            >  6:212643/1‑149 (MQ=255)
                                                                                                      gTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGATGTGCt                                        >  2:108721/1‑149 (MQ=255)
                                                                                                       tCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGATGTGCtt                                       <  4:152416/149‑1 (MQ=255)
                                                                                                         ggAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGATGTGCTtgt                                     <  7:197190/149‑1 (MQ=255)
                                                                                                          gAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGATGTGCTtgtg                                    >  8:169667/1‑149 (MQ=255)
                                                                                                            ggTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGATGTGCTTGTGtt                                  <  4:201131/149‑1 (MQ=255)
                                                                                                             gTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGATGTGCTTGTGTTg                                 >  1:143892/1‑149 (MQ=255)
                                                                                                             gTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGATGTGCTTGTGTTg                                 >  6:289193/1‑149 (MQ=255)
                                                                                                             gTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGATGTGCTTGTGTTg                                 >  3:163979/1‑149 (MQ=255)
                                                                                                              tCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGATGTGCTTGTGTTGc                                <  1:322135/149‑1 (MQ=255)
                                                                                                                aGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGATGTGCTTGTGTTGCa                               <  5:289193/148‑1 (MQ=255)
                                                                                                                 ggATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGGTTGAGGACGATGTGCTTGTGTTGCAt                              >  8:301611/1‑148 (MQ=255)
                                                                                                                  gATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGATGTGCTTGTGTTGCATGa                            >  4:174694/1‑149 (MQ=255)
                                                                                                                   aTGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGTCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGATGTGCTTGTGTTGCATGAc                           <  7:199535/149‑1 (MQ=255)
                                                                                                                        cTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGATGTGCTTGTGTTg                                 <  1:136092/138‑1 (MQ=255)
                                                                                                                        cTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGATGTGCTTGTGTTg                                 >  2:136092/1‑138 (MQ=255)
                                                                                                                             cGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGATGTGCTTGTGTTGCATGACGGCAAAACCt                 >  5:225181/1‑149 (MQ=255)
                                                                                                                                  tCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGttt                                                      <  4:123711/107‑1 (MQ=255)
                                                                                                                                  tCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGttt                                                      >  3:123711/1‑107 (MQ=255)
                                                                                                                                  tCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGATGTGCTTGTGTTGCATGACGGCAAAACCTCCgg             <  4:154174/148‑1 (MQ=255)
                                                                                                                                     gCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGATGTGCTTGTGTTGCATGACGGCAAAACCTCCGGGTAc         <  6:225181/149‑1 (MQ=255)
                                                                                                                                     gCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCAGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGATGTGCTTGTGTTGCATGACGGCAAAACCTCCGGGTAc         >  3:58846/1‑149 (MQ=255)
                                                                                                                                      ccAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGATGTGCTTGTGTTGCATGACGGCAAAACCTCCGGGTACt        <  7:169667/149‑1 (MQ=255)
                                                                                                                                      ccAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGATGTGCTTGTGTTGCATGACGGCAAAACCTCCGGGTACt        <  2:179331/149‑1 (MQ=255)
                                                                                                                                            tGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGATGTGCTTGTGTTGCATGACGGCAAAACCTCCGGGTACTGGTcag  <  7:213560/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GGCATCGGTAGCCGCGAAGCGCTCGAGGTCGGTCAGCTGATACACGCGCATCACCACCGACAGCGGCTTGCCCTCGGCGTCGGGGTTGAGCTGGTTGCCGCCGTCGGAGGTCAGGATGACCTTTTCGCTGTCGGCCAGCATGTCGGCCGCCCAGCTGTCTTCCATGCGCTTGCCGATACGGTCGGTGACCCCGCAGCCGGCCAGCGCTACCAGCAAGGCGGCTGCGGTCAGGCGTTTGAGGACGATGTGCTTGTGTTGCATGACGGCAAAACCTCCGGGTACTGGTCAG  >  NC_002947/3488469‑3488757

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: ‑