Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,573,889 C→T 5.4% A139A (GCG→GCA eutT ← cobalamin adenosyltransferase involved in ethanolamine utilization

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,573,8890CT 5.4% 236.8 / 2.6 74A139A (GCG→GCAeutTcobalamin adenosyltransferase involved in ethanolamine utilization
Reads supporting (aligned to +/- strand):  ref base C (35/35);  new base T (2/2);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

CGTGGTCGAGATAGCGCAGCGGCTGATGAGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATC  >  NC_000913/2573753‑2574032
                                                                                                                                        |                                                                                                                                               
cGTGGTCGAGATAGCGCAGCGGCTGATGAGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACgg                                                                                                                                      <  4:178894/148‑1 (MQ=255)
 gTGGTCGAGATAGCGCAGCGGCTGATGAGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGc                                                                                                                                     <  3:354468/148‑1 (MQ=255)
  tGGTCGAGATAGCGCAGCGGCTGATGAGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTg                                                                                                                                   >  6:74574/1‑149 (MQ=255)
     tcgttaTAGCGCATCTGCTGATGATAAAGAATATGCATATCTTCGTCAGAAATCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCa                                                                                                                                <  8:307290/144‑1 (MQ=255)
     tCGAGATAGCGCAGCGGCTGATGAGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCa                                                                                                                                >  1:349104/1‑149 (MQ=255)
       gagaTAGCGCAGCGGCTGATGAGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACg                                                                                                                              >  1:429425/1‑149 (MQ=255)
             gcgcAGCGGCTGATGAGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTcc                                                                                                                         <  2:291958/148‑1 (MQ=255)
               gcAGCGGCTGATGAGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGc                                                                                                                       >  3:30622/1‑148 (MQ=255)
                        gatgaGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCTGCCAGCCACGGCTGCCACGGTTCCGCCAGtt                                                                                                                  >  6:151129/1‑144 (MQ=255)
                        gatgaGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCTGCCAGCCACGGCTGCCACGGTTCCGCCAGtt                                                                                                                  <  5:151129/144‑1 (MQ=255)
                         atgaGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAAtt                                                                                                             <  1:106670/148‑1 (MQ=255)
                         atgaGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAAt                                                                                                              <  7:288723/147‑1 (MQ=255)
                          tgaGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTg                                                                                                           >  4:75982/1‑149 (MQ=255)
                             gaAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAg                                                                                                        <  2:429425/149‑1 (MQ=255)
                              aaaGCCGATGCAGATCTTCGTCAGAAAGCCCGGCAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCCGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGc                                                                                                       >  7:87051/1‑149 (MQ=255)
                              aaaGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCCCGGTTCCGCCAGTTCAATTTGCAGc                                                                                                       >  6:289763/1‑149 (MQ=255)
                              aaaGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCa                                                                                                         <  3:71697/147‑1 (MQ=255)
                              aaaGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCa                                                                                                         >  4:71697/1‑147 (MQ=255)
                                 gCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCca                                                                                                     <  5:222441/148‑1 (MQ=255)
                                 gCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCca                                                                                                     >  6:222441/1‑148 (MQ=255)
                                  ccGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCcaca                                                                                                   <  2:349104/149‑1 (MQ=255)
                                  ccGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCcaca                                                                                                   >  7:216388/1‑149 (MQ=255)
                                      tGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGc                                                                                               <  3:75982/149‑1 (MQ=255)
                                         aGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGGATAATGTTGCCAAGACGCGAACGGATATCTGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAg                                                                                            >  7:290867/1‑149 (MQ=255)
                                           aTCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAgcg                                                                                          >  1:147442/1‑149 (MQ=255)
                                                  tCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGcc                                                                                              <  6:182688/138‑1 (MQ=255)
                                                  tCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGcc                                                                                              >  5:182688/1‑138 (MQ=255)
                                                     gAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTTCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCt                                                                                >  1:300017/1‑149 (MQ=255)
                                                      aaaGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTg                                                                               <  3:268317/149‑1 (MQ=255)
                                                          cccGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCga                                                                           <  1:375342/149‑1 (MQ=255)
                                                           ccGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTc                                                                             >  2:2474/1‑146 (MQ=255)
                                                           ccGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTc                                                                             >  4:61420/1‑146 (MQ=255)
                                                               cAATCGCCTGGCAGCCCAGCGGTTCCCACAGCGAATCGGCGCGCATAATGTTGCTAAGACGCGAACGTATATCCGCCAGCCACGGCTGCCACGGTTACGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAAcc                                                                      <  5:289763/149‑1 (MQ=255)
                                                                        ggCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGaaaa                                                             >  8:223296/1‑149 (MQ=255)
                                                                         gCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAAc                                                            >  8:2671/1‑149 (MQ=255)
                                                                            gCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCa                                                         <  7:2671/149‑1 (MQ=255)
                                                                             cccAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACCCTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAg                                                        >  5:405489/1‑149 (MQ=255)
                                                                              ccAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCTAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCGGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGa                                                       >  1:186466/1‑149 (MQ=255)
                                                                              ccAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGGTGCCAAGACGCGAACGGATATCTGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGa                                                       <  8:290867/149‑1 (MQ=255)
                                                                                    gTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGa                                                 >  4:4230/1‑149 (MQ=255)
                                                                                     ttCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGGGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGAt                                                <  8:148442/149‑1 (MQ=255)
                                                                                          ccAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAACCCCAGACGCGGATCGCtt                                           >  4:20565/1‑149 (MQ=255)
                                                                                          ccAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCtt                                           >  5:191055/1‑149 (MQ=255)
                                                                                          ccAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCtt                                           >  3:302966/1‑149 (MQ=255)
                                                                                          ccAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCtt                                           >  1:246002/1‑149 (MQ=255)
                                                                                          ccAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCtt                                           >  7:313299/1‑149 (MQ=255)
                                                                                           cAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCa                                                                                                                     <  1:422925/74‑1 (MQ=255)
                                                                                           cAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGACa                                                                                                                     >  2:422925/1‑74 (MQ=255)
                                                                                            aGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCtttt                                         <  4:30622/149‑1 (MQ=255)
                                                                                             gcgcATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTg                                        <  3:61420/149‑1 (MQ=255)
                                                                                              cgcATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTgg                                       <  8:216388/149‑1 (MQ=255)
                                                                                                cATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGc                                      >  2:267462/1‑148 (MQ=255)
                                                                                                 aTCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGa                                    <  1:2474/149‑1 (MQ=255)
                                                                                                  tCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGAc                                   <  2:147442/149‑1 (MQ=255)
                                                                                                       gcgcATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGcc                                                                                                                      >  6:96755/1‑61 (MQ=255)
                                                                                                       gcgcATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGcc                                                                                                                      <  5:96755/61‑1 (MQ=255)
                                                                                                         gcATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCga                                                                           >  8:422181/1‑102 (MQ=255)
                                                                                                         gcATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCAttttt                            >  5:121703/1‑149 (MQ=255)
                                                                                                           aTAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCAttttt                            <  8:276991/147‑1 (MQ=255)
                                                                                                           aTAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCAttttt                            >  7:276991/1‑147 (MQ=255)
                                                                                                           aTAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTc                          <  8:87051/149‑1 (MQ=255)
                                                                                                           aTAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTc                          >  6:199636/1‑149 (MQ=255)
                                                                                                            tAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTcc                         <  6:191055/149‑1 (MQ=255)
                                                                                                             aaTGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTccg                        >  6:359352/1‑149 (MQ=255)
                                                                                                                 ttGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGa                    >  8:408756/1‑149 (MQ=255)
                                                                                                                   gCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGt                                                                              >  6:316801/1‑89 (MQ=255)
                                                                                                                   gCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGt                                                                              <  5:316801/89‑1 (MQ=255)
                                                                                                                    ccAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTGCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAg                 <  6:360527/149‑1 (MQ=255)
                                                                                                                      aaGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTTCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGt               <  2:300017/149‑1 (MQ=255)
                                                                                                                      aaGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGt               <  7:408756/149‑1 (MQ=255)
                                                                                                                        gACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGa             <  3:20565/149‑1 (MQ=255)
                                                                                                                              aaCGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGc       <  5:74574/149‑1 (MQ=255)
                                                                                                                              aaCGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGc       <  7:223296/149‑1 (MQ=255)
                                                                                                                                   atatCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATc  <  5:199636/149‑1 (MQ=255)
                                                                                                                                        cGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTa    <  7:430142/142‑1 (MQ=255)
                                                                                                                                        cGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTa    >  8:243722/1‑142 (MQ=255)
                                                                                                                                        cGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTa    <  7:243722/142‑1 (MQ=255)
                                                                                                                                        cGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTa    >  8:430142/1‑142 (MQ=255)
                                                                                                                                        |                                                                                                                                               
CGTGGTCGAGATAGCGCAGCGGCTGATGAGAAAGCCGATGCAGATCTTCGTCAGAAAGCCCGACAATCGCCTGGCAGCCCAGCGGTTCCCCCAGCGCATCGGCGCGCATAATGTTGCCAAGACGCGAACGGATATCCGCCAGCCACGGCTGCCACGGTTCCGCCAGTTCAATTTGCAGCCACACTGCCAGCGCAATGGTGCTGTCGAGAACCGCGCGAAAACCCAGACGCGGATCGCTTTTGGCGACCATTTTTTCCGCCGACAGGTGAGTCAGCGTATC  >  NC_000913/2573753‑2574032

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