Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,160,113 G→T H226Q (CAC→CAA sthA ← pyridine nucleotide transhydrogenase, soluble

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,160,1130GT98.7% 222.7 / NA 76H226Q (CAC→CAAsthApyridine nucleotide transhydrogenase, soluble
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (29/46);  minor base C (1/0);  total (30/46)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATA  >  NC_000913/4159976‑4160251
                                                                                                                                         |                                                                                                                                          
gttcgcccGTTGGCATCGAGCATGCAGTTAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATa                                                                                                                                       <  2:669574/137‑1 (MQ=255)
ggtgggacGTTGGCATAGAGCAGGCAGGCAGCTTGCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATa                                                                                                                                       <  2:714683/136‑1 (MQ=255)
 tGCGACCGGTGGCATAGAGCAGGCAGTCAGCTTCCAGTTTTTTACCCGGCTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAg                                                                                                                                      <  2:1026283/143‑1 (MQ=255)
   cGATCGTTGGCATAGAGCAGGCAGTCAGCTTTGAGTTTTTTACCCGACTTCAGATGCAGGATCACACCAGCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGga                                                                                                                                    <  1:1029867/143‑1 (MQ=255)
     aCCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGACTCACTACGCCACTGTTCTAGAACTGATAGgaga                                                                                                                                  >  1:820438/1‑143 (MQ=255)
       cGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCGTCGTCACAGCCTTCGCTCTTCTCGTACTCTTCGTTGTGACGCATCACTACGTCACTGTTCCAGAGTTGATCGgagaga                                                                                                                                >  1:179704/1‑143 (MQ=255)
       cGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTTGGACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGgagaga                                                                                                                                >  2:965807/1‑143 (MQ=255)
       cGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGgagaga                                                                                                                                >  1:829155/1‑143 (MQ=255)
         ttGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGgagagcga                                                                                                                              >  2:209372/1‑143 (MQ=255)
             cATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCt                                                                                                                          <  2:196868/143‑1 (MQ=255)
              aTAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACCTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTg                                                                                                                         >  1:919817/1‑143 (MQ=255)
               tAGAGCAGGGAGTCAGCTTTGCGTCTTTTACCCGACTGCAGAGGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGa                                                                                                                        <  2:356567/143‑1 (MQ=255)
               tAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGa                                                                                                                        <  1:582499/143‑1 (MQ=255)
                agagCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGAc                                                                                                                       <  1:1010275/143‑1 (MQ=255)
                  agCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGACTCTGACAt                                                                                                                     >  1:5016/1‑143 (MQ=255)
                  agCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACAt                                                                                                                     >  2:808979/1‑143 (MQ=255)
                     aggcagTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACAtctc                                                                                                                  <  2:931227/143‑1 (MQ=255)
                         agTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGACTTGATAGGAGAGAGAATCTGACATCTCTTGa                                                                                                              >  1:644437/1‑143 (MQ=255)
                          gTCAGCTTTCAGTTTTTTACCCGACTCTAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGAt                                                                                                             <  2:658913/143‑1 (MQ=255)
                           tCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACATCCTTCGATCTTCTCGTATTCTTCGTTGTGACGAATTACTCCGCCACTGTTCCAGAATTGTTAGGAGAGACAATCTGCCATCACTTGCTc                                                                                                            >  2:814208/1‑143 (MQ=255)
                               cTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCACAATTGATACGAGAGAGAATCTGACATCTCTTGACCgaga                                                                                                        >  2:402842/1‑143 (MQ=255)
                               cTTTCAGTTCTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCgaga                                                                                                        <  1:9946/143‑1 (MQ=255)
                                 ttCAGTTTTTTACCCGACTTCAGATGCATGATCACACCAGCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGaaa                                                                                                      <  2:410144/143‑1 (MQ=255)
                                    aGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGc                                                                                                   <  2:748077/143‑1 (MQ=255)
                                     gTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTCTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGACTCTGACATCTCTTGAGCGAGAAATGcc                                                                                                  >  2:707095/1‑143 (MQ=255)
                                       tttttACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCcag                                                                                                <  2:575321/143‑1 (MQ=255)
                                         tttACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAAGCTGACATCTCTTGATCGAGTACTGCcagca                                                                                              >  2:895748/1‑143 (MQ=255)
                                          ttACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCcagcag                                                                                             <  1:520078/143‑1 (MQ=255)
                                              ccGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCCCTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGCACTTGCAGCAGGCgg                                                                                         >  2:658785/1‑142 (MQ=255)
                                              ccGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTAACGAATAAGTGCCAGCAGGCGa                                                                                         >  2:266918/1‑143 (MQ=255)
                                                 aCTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATcg                                                                                      <  1:375885/143‑1 (MQ=255)
                                                    tCAGCTGCATCATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCgg                                                                                   <  2:1001118/143‑1 (MQ=255)
                                                    tCAGATGCATGATCACACCATCGTCACAGTCTTTGTTTTTCTTGGACTGTACGTTGTGACGAAGCACTACGCCACGGGTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCgg                                                                                   <  1:962268/143‑1 (MQ=255)
                                                    tCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCgg                                                                                   <  1:191544/143‑1 (MQ=255)
                                                      aGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTATGACGAATCAGTCAGCCACTGCTTCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCggat                                                                                 >  2:132250/1‑141 (MQ=255)
                                                            aTGATCACACCATCGTCACAGCATTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGAt                                                                           <  1:706604/143‑1 (MQ=255)
                                                             tGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATc                                                                          <  1:592084/143‑1 (MQ=255)
                                                                tCACACCCTCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGa                                                                       <  2:874266/143‑1 (MQ=255)
                                                                 caccCCATCGTCCCAGCCTTCGATCTTCTCGTACTCTTCTTTGTGACGAAGCACTACGCCACTATTCCAGAATTGATAGGAGCGAGATTCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGAt                                                                      <  2:1063336/143‑1 (MQ=255)
                                                                    acCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCCGCACGCGATCGCGGGTGTTGATCAGATACa                                                                   >  1:507212/1‑143 (MQ=255)
                                                                      cagcGTCCCAGCCTTCGATCTTCTCGTACTCCTCGTTGTGACGAATCACTCCGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACt                                                                 <  2:669248/140‑1 (MQ=255)
                                                                         cGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACtttt                                                              >  2:91713/1‑143 (MQ=255)
                                                                          gTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTa                                                             <  1:931938/143‑1 (MQ=255)
                                                                          gTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTa                                                             <  2:691142/143‑1 (MQ=255)
                                                                            cacgGCCTTCGATCTCCTCGTACTCTTCGTTGTGACGAATCACTGCGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACa                                                           <  2:1037017/143‑1 (MQ=255)
                                                                             acaGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGAGTTGATCAGAACCACTTTTACag                                                          >  2:58892/1‑142 (MQ=255)
                                                                               aGCCTTCGATCTTCTCGTACTCATCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATcc                                                        <  2:937311/143‑1 (MQ=255)
                                                                                 ccTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGACAGAGAATCTGACATCACTTGATCGAGAAATGACAACAGGAGATCGCGGGTGTTGATCAGATCCACTTTTACATCCAt                                                      >  1:640074/1‑143 (MQ=255)
                                                                                   ttCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATAc                                                    >  1:525008/1‑143 (MQ=255)
                                                                                       aTCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACcgcg                                                <  2:300600/143‑1 (MQ=255)
                                                                                          ttctCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTGCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGaa                                             <  2:595743/143‑1 (MQ=255)
                                                                                            ctcGTTCTCTTCGTTGTGACGAATCACTACTCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGa                                           <  2:466930/143‑1 (MQ=255)
                                                                                            ctcGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGa                                           >  1:1031531/1‑143 (MQ=255)
                                                                                             tcGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTCGATCGAGAAATGCCAGCAGGCGCTCGCGGGTCTCGATCAGCTACACGTCTCCACCCATCCCGCgggacac                                          >  2:1007017/1‑138 (MQ=255)
                                                                                              cGTACTCTTCGTTGTGACGAAGCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATc                                         <  1:808979/143‑1 (MQ=255)
                                                                                                 gctctGCGTTGTGAGGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGAc                                      <  2:258692/142‑1 (MQ=255)
                                                                                                 aCTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGAc                                      <  2:919817/143‑1 (MQ=255)
                                                                                                   tctTCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGc                                    <  1:437010/143‑1 (MQ=255)
                                                                                                   tctTCGTTGTGACGAATCACTACGCCACGGTGGCAGAATGGAAAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGc                                    <  1:314374/143‑1 (MQ=255)
                                                                                                     ttCGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCat                                  >  2:836455/1‑143 (MQ=255)
                                                                                                       cGTTGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCatat                                <  2:469291/143‑1 (MQ=255)
                                                                                                         ttGTGACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCCGATTc                              >  1:155700/1‑143 (MQ=255)
                                                                                                          cgtgACGAAGCGCGCCGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTca                             <  1:90652/142‑1 (MQ=255)
                                                                                                             gACGAATCACTAGGCCACTGTTGCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAg                          <  1:332073/143‑1 (MQ=255)
                                                                                                             gACGAATCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAg                          <  1:91713/143‑1 (MQ=255)
                                                                                                                   tCACTACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATc                    >  1:453341/1‑143 (MQ=255)
                                                                                                                       tACGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGTCCATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTc                >  1:495480/1‑143 (MQ=255)
                                                                                                                       tACGCCACTGTTCCAGAATTGAGAGGAGAGAGAATCTGACATCTCTTGAGCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTc                <  1:757032/143‑1 (MQ=255)
                                                                                                                         cgccgcTGTTCCAGAGTTGCTAGGAGAGGGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCa              <  2:399427/143‑1 (MQ=255)
                                                                                                                         cGCCACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGATGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCa              <  1:361066/143‑1 (MQ=255)
                                                                                                                            cACTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGAGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATGGACGCATATTAACAGCGGACAACTcacgca           >  2:398342/1‑138 (MQ=255)
                                                                                                                             aCTGTTCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCAc          <  1:965807/143‑1 (MQ=255)
                                                                                                                                 ttCCAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTAACTCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATa      >  1:508652/1‑143 (MQ=255)
                                                                                                                                 ttCCAGAATGGATAGGAGAGCGAATCTTACATCTCTTGATCTAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATa      <  2:894214/143‑1 (MQ=255)
                                                                                                                                    cAGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACACATATTCACAGCCGATCACTCCAGCACAATAGAt   >  1:899221/1‑143 (MQ=255)
                                                                                                                                     aGAATTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATa  <  1:58892/143‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                          
GTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATA  >  NC_000913/4159976‑4160251

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