Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,758,650 A→G 5.2% I404I (ATT→ATC selB ← selenocysteyl‑tRNA‑specific translation elongation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,758,6500AG 5.2% 301.5 / 4.7 96I404I (ATT→ATCselBselenocysteyl‑tRNA‑specific translation elongation factor
Reads supporting (aligned to +/- strand):  ref base A (49/42);  new base G (2/3);  total (51/45)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

CTTCCATTGGCAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACGCCAGCCA  >  NC_000913/3758502‑3758799
                                                                                                                                                    |                                                                                                                                                     
cTTCCATTGGCAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGaa                                                                                                                                                      <  1:1218256/150‑1 (MQ=255)
cTTCAATTGGCAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGaa                                                                                                                                                      <  1:1144173/150‑1 (MQ=255)
 ttCCATTGGCAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAat                                                                                                                                                     <  1:857291/150‑1 (MQ=255)
  tCCATTGGCAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAata                                                                                                                                                    >  2:1372613/1‑150 (MQ=255)
  tCCATTGGCAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAata                                                                                                                                                    >  2:778430/1‑150 (MQ=255)
   ccATTGGCAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAatat                                                                                                                                                   >  1:1166158/1‑150 (MQ=255)
     aTTGGCAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATTTGTCGAGAATTTTCCGCTTCCAGCGGGCGGCAACGGCGCATTCAACAAGCTATAACCAGCCTGAATATAAc                                                                                                                                                <  1:1428775/150‑1 (MQ=255)
     aTTGGCAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATaa                                                                                                                                                 >  1:1130154/1‑150 (MQ=255)
        ggCAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCa                                                                                                                                              >  2:119911/1‑150 (MQ=255)
             cGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTg                                                                                                                                         <  1:307257/150‑1 (MQ=255)
             cGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATAGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTg                                                                                                                                         >  2:394941/1‑150 (MQ=255)
               ccATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGtt                                                                                                                                       <  1:1520549/150‑1 (MQ=255)
                 aTACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGGTTATATTCAGGCTggtta                                                                                                                                     <  2:650016/150‑5 (MQ=255)
                  tACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTgca                                                                                                                                    <  2:923543/150‑1 (MQ=255)
                  tACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGACAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGGTTATATTCAGGCTggttat                                                                                                                                    <  2:979921/150‑6 (MQ=255)
                   tcGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTGATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCTGCAACCTGCGCATTCATCAAGCTATAACCAGCCTGAATATAACCAGGCTGTTgcag                                                                                                                                   <  2:637512/149‑1 (MQ=255)
                    cGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGGTTATATTCAGGCTggttatag                                                                                                                                  <  1:45492/150‑8 (MQ=255)
                        gCAGACGTTCGCGCCCAGGGCAAGGTTCATCGCTATGAAGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTACGCTGCCAGAGGGCGGAAACCGGCGGATACAACAAGCTATAAACAGAATGAATATAACCAGGCTGTTGCAGCAAtt                                                                                                                              >  2:1382859/1‑150 (MQ=255)
                             cGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCgcgc                                                                                                                         <  1:12580/150‑1 (MQ=255)
                              gTTCGCGCCCAGGGCCAGGTTCATCGCGATCTTGCTCATGATAAGTCGCTAATGTCTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACAGGCGCATTCAACAAGCTATAACCAGCCTGAAAATAACCAGGCTGTTGAAGCAATTCGCGCa                                                                                                                        >  2:29734/1‑150 (MQ=255)
                                      ccAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCg                                                                                                                <  1:871144/150‑1 (MQ=255)
                                        aGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGcc                                                                                                              >  2:942762/1‑150 (MQ=255)
                                              aGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAg                                                                                                        >  2:997080/1‑150 (MQ=255)
                                                   cATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGgc                                                                                                   <  1:119911/150‑1 (MQ=255)
                                                      cgcgATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGgcgcg                                                                                                <  1:1361834/150‑1 (MQ=255)
                                                            gTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAgg                                                                                          >  1:236954/1‑150 (MQ=255)
                                                            gTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATCCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAgg                                                                                          >  2:630924/1‑150 (MQ=255)
                                                              tGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCg                                                                                        <  1:778430/150‑1 (MQ=255)
                                                               gCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGa                                                                                       >  1:395360/1‑150 (MQ=255)
                                                                cTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGTAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGaa                                                                                      <  1:1246185/150‑1 (MQ=255)
                                                                cTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGaa                                                                                      >  1:655266/1‑150 (MQ=255)
                                                                 tCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAATCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGaaa                                                                                     <  2:1532404/150‑1 (MQ=255)
                                                                    tGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATcc                                                                                  >  2:1296398/1‑150 (MQ=255)
                                                                      aTAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGc                                                                                <  1:518227/150‑1 (MQ=255)
                                                                       tAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCa                                                                               >  2:70559/1‑150 (MQ=255)
                                                                         aGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAg                                                                             >  1:1375090/1‑150 (MQ=255)
                                                                          gTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCa                                                                                                                        <  2:724348/106‑1 (MQ=255)
                                                                          gTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCa                                                                                                                        >  1:724348/1‑106 (MQ=255)
                                                                            cGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGtt                                                                          >  2:152461/1‑150 (MQ=255)
                                                                            cGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGtt                                                                          >  1:883216/1‑150 (MQ=255)
                                                                            cGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGtt                                                                          >  1:608115/1‑150 (MQ=255)
                                                                             gCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTa                                                                         >  1:1046800/1‑150 (MQ=255)
                                                                              cTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTaa                                                                        >  1:675035/1‑150 (MQ=255)
                                                                                     gtCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAAccgcgcc                                                                 >  1:270580/1‑150 (MQ=255)
                                                                                     gtCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAAccgcgcc                                                                 >  2:1152032/1‑150 (MQ=255)
                                                                                      tCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAAccgcgccg                                                                <  1:845476/150‑1 (MQ=255)
                                                                                       cGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAAccgcgccgc                                                               >  1:101308/1‑150 (MQ=255)
                                                                                         agaATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGt                                                             >  2:989760/1‑150 (MQ=255)
                                                                                          gaATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGtt                                                            >  1:52202/1‑150 (MQ=255)
                                                                                              tttCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAg                                                        <  1:394941/150‑1 (MQ=255)
                                                                                              tttCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAg                                                        <  1:1382859/150‑1 (MQ=255)
                                                                                              tttCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAg                                                        <  1:1427280/150‑1 (MQ=255)
                                                                                              tttCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAg                                                        <  2:1283004/150‑1 (MQ=255)
                                                                                              tttCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTAAAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAg                                                        <  1:1429790/150‑1 (MQ=255)
                                                                                                 ccGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATg                                                     <  2:627778/150‑1 (MQ=255)
                                                                                                 ccGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATg                                                     <  2:1130154/150‑1 (MQ=255)
                                                                                                  cGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGa                                                    >  2:500735/1‑150 (MQ=255)
                                                                                                   gCTGCCAGGGGGCGGCAACCGGCTCATTCAACAAGCTATAACCAGCCTGAATATATCCAGGATGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCGGGCGCGCGCAGGCTAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGaa                                                   <  2:1378967/150‑1 (MQ=255)
                                                                                                   gCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGaa                                                   >  2:136173/1‑150 (MQ=255)
                                                                                                    cTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAAc                                                  >  1:708556/1‑150 (MQ=255)
                                                                                                     tGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACa                                                 >  2:1166355/1‑150 (MQ=255)
                                                                                                     tGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACa                                                 <  2:240097/150‑1 (MQ=255)
                                                                                                      gCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAg                                                <  1:201100/150‑1 (MQ=255)
                                                                                                      cccAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACAAGCCTGAATATACCCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAg                                                <  2:293891/149‑1 (MQ=255)
                                                                                                        cAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGACTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGAt                                              >  2:343224/1‑150 (MQ=255)
                                                                                                          gcgggcggCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATaa                                            >  2:47523/1‑150 (MQ=255)
                                                                                                              gaggcAACCGGCGCATTCAACAAGCTATAACCAGCCTTAATATAACCAGTCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACgca                                        <  1:566009/148‑1 (MQ=255)
                                                                                                               cggcAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACgcat                                       <  1:1397424/150‑1 (MQ=255)
                                                                                                                   aaCCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGCTTAACCGCGCCGCGTTCCAGATGAACAGATAACgcatcggc                                   >  1:249630/1‑150 (MQ=255)
                                                                                                                    aCCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACgcatcggca                                  <  1:989760/150‑1 (MQ=255)
                                                                                                                    aCCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACgcatcggca                                  <  1:792156/150‑1 (MQ=255)
                                                                                                                        gcgcATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGc                              <  2:1166158/150‑1 (MQ=255)
                                                                                                                          gcATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGctc                            <  2:1273951/150‑1 (MQ=255)
                                                                                                                          gcATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGctc                            >  2:240141/1‑150 (MQ=255)
                                                                                                                          gcATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGctc                            >  2:1212631/1‑150 (MQ=255)
                                                                                                                             ttCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTgc                         >  1:520366/1‑150 (MQ=255)
                                                                                                                             ttCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATCCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTgc                         <  1:630924/150‑1 (MQ=255)
                                                                                                                                    aGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTg                  <  1:1372613/150‑1 (MQ=255)
                                                                                                                                     gCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGc                 <  1:466747/150‑1 (MQ=255)
                                                                                                                                     gCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGc                 <  2:599733/150‑1 (MQ=255)
                                                                                                                                      cTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCa                <  1:1132781/150‑1 (MQ=255)
                                                                                                                                       tataACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCaa               >  1:934996/1‑150 (MQ=255)
                                                                                                                                       tataACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCaa               >  2:914393/1‑150 (MQ=255)
                                                                                                                                         taACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAaga             >  2:213012/1‑150 (MQ=255)
                                                                                                                                         taACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAaga             <  1:1163935/150‑1 (MQ=255)
                                                                                                                                         gcaacAGCCTGGTTATATTCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAaga             >  2:1500575/5‑150 (MQ=255)
                                                                                                                                          aaCCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAagag            <  2:1141678/150‑1 (MQ=255)
                                                                                                                                          aaCCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAagag            >  1:1163156/1‑150 (MQ=255)
                                                                                                                                          aaCCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAagag            >  2:744968/1‑150 (MQ=255)
                                                                                                                                           aCCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAagag            <  2:1138381/149‑1 (MQ=255)
                                                                                                                                           aCCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAagag            >  1:1138381/1‑149 (MQ=255)
                                                                                                                                           aCCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAgaga           >  1:396721/1‑150 (MQ=255)
                                                                                                                                            ccAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGAc          >  1:772740/1‑150 (MQ=255)
                                                                                                                                            ccAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGAc          <  2:608115/150‑1 (MQ=255)
                                                                                                                                            acAGCCTGGTTATATTCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGAc          >  2:1267374/2‑150 (MQ=255)
                                                                                                                                             cAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACg         >  1:379008/1‑150 (MQ=255)
                                                                                                                                                  tGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGATCTGGCTAGCCAATTCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACgccagc    >  2:1252194/1‑150 (MQ=255)
                                                                                                                                                  tGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACgccagc    >  2:1072882/1‑150 (MQ=255)
                                                                                                                                                    aaTATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACgccagcca  <  2:381403/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTTCCATTGGCAACGCCATACGTCGCAGACGTTCGCGCCCAGGGCCAGGTTCATCGCGATGTTGCTCATGATAAGTCGCTAATGTGTCGAGAATTTTCCGCTGCCAGCGGGCGGCAACCGGCGCATTCAACAAGCTATAACCAGCCTGAATATAACCAGGCTGTTGCAGCAATTCGCGCATCCCTTCGCCGTTGAGCTGGCGCGCCCAGGCGAAATCCGCAAGGTTAACCGCGCCGCGTTCCAGATGAACAGATAACGCATCGGCATCGCTCTGCGCCCGTGCAAGAGACGCCAGCCA  >  NC_000913/3758502‑3758799

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