Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,941,055 C→T 21.1% H154H (CAC→CAT)  PP_4348 → cystathionine beta‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,941,0550CT21.1% 152.9 / 34.8 80H154H (CAC→CAT) PP_4348cystathionine beta‑lyase
Reads supporting (aligned to +/- strand):  ref base C (36/27);  new base T (10/7);  total (46/34)
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.98e-01

TAAAAATTGAATACTTCACACCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACCACGGGGCTGTAGAGGGCATCTCTCTCGGAGCACTAATAACCAAGCGAAA  >  NC_002947/4940907‑4941204
                                                                                                                                                    |                                                                                                                                                     
tAAAAATTGAATACTTCACACCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGaccacc                                                                                                                                                      >  2:2257166/1‑150 (MQ=255)
  aaaaTTGAATACTTCACACCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCa                                                                                                                                                    >  1:1318573/1‑150 (MQ=255)
         aaTACTTCACACCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGa                                                                                                                                             >  1:997344/1‑150 (MQ=255)
          aTACTTCAAACCTATCGAAGTCGCGACTCTAAAAAATCAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATCGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGAt                                                                                                                                            >  2:699440/1‑150 (MQ=255)
           tACTTCACACCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCATCCAGCGCTGAtt                                                                                                                                           >  1:1904141/1‑150 (MQ=255)
           tACTTCACACCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGAtt                                                                                                                                           <  2:1028887/150‑1 (MQ=255)
               tCACACCTAGCGAAGCCGCGACTCTCAAAAATAAACTAACGCTAAAGACAAAATTAATCATAGTTGAAAACCCAGACGCATTCACGTACGAGGTCCTGGATATTAGTTCGATCGTCGAATCAAGACACGGCCACCGAGCGCTGATTGTgg                                                                                                                                       >  1:1058731/1‑150 (MQ=255)
               tCACACCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTgg                                                                                                                                       <  2:997344/150‑1 (MQ=255)
                cacaCCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGc                                                                                                                                      <  2:1318573/150‑1 (MQ=255)
                 acacCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATGTTAGTGCGATCGTCGAATCAAGACACGACCATCCAGCGCTGATTGTGGcc                                                                                                                                     <  2:946349/150‑1 (MQ=255)
                 acacCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGcc                                                                                                                                     >  2:1519786/1‑150 (MQ=255)
                     cTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGaca                                                                                                                                 >  2:818033/1‑150 (MQ=255)
                       aGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGacaac                                                                                                                               <  1:2106698/150‑1 (MQ=255)
                            aGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATg                                                                                                                          <  2:1044782/150‑1 (MQ=255)
                              ccGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCATCCAGCGCTGATtg                                                                                                                                          >  1:647548/1‑132 (MQ=255)
                              ccGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCATCCAGCGCTGATtg                                                                                                                                          <  2:647548/132‑1 (MQ=255)
                              ccGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATggg                                                                                                                        >  1:1969850/1‑150 (MQ=255)
                                gcgACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGgca                                                                                                                      >  1:2011383/1‑150 (MQ=255)
                                     tctAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGACTTAACGTACGAGGTGATATGGATTAGGGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATc                                                                                                                 >  1:1451573/1‑150 (MQ=255)
                                        aaaaaaTAAACTAACTCTAAAGACAAAATTAATCATATTTGAATACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACt                                                                                                              >  2:1695795/1‑150 (MQ=255)
                                             aTAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTc                                                                                                         <  1:1883260/150‑1 (MQ=255)
                                                 actaactCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTcacca                                                                                                     >  1:511413/1‑150 (MQ=255)
                                                  ctaactCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGAtt                                                                                                                                           <  1:1216519/111‑1 (MQ=255)
                                                  ctaactCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGAtt                                                                                                                                           >  2:1216519/1‑111 (MQ=255)
                                                   taactCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCATCCAGCGCTGAt                                                                                                                                            <  2:1858151/109‑1 (MQ=255)
                                                   taactCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCATCCAGCGCTGAt                                                                                                                                            >  1:1858151/1‑109 (MQ=255)
                                                     actCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCa                                                                                                                                                    <  2:1775436/99‑1 (MQ=255)
                                                     actCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCa                                                                                                                                                    >  1:1775436/1‑99 (MQ=255)
                                                         tAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTAt                                                                                             >  1:534205/1‑150 (MQ=255)
                                                              aCAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCATCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTa                                                                                              <  1:1185752/144‑1 (MQ=255)
                                                              aCAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCATCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTa                                                                                              >  2:1185752/1‑144 (MQ=255)
                                                              aCAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTc                                                                                        >  1:1298234/1‑150 (MQ=255)
                                                                aaaaTTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCgg                                                                                      >  1:2217337/1‑150 (MQ=255)
                                                                aaaaTTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCg                                                                                       >  2:423262/1‑149 (MQ=255)
                                                                aaaaTTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCg                                                                                       <  1:423262/149‑1 (MQ=255)
                                                                   aTTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCATCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGgcgc                                                                                   >  2:1070378/1‑150 (MQ=255)
                                                                   aTTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTAt                                                                                             <  2:374045/140‑1 (MQ=255)
                                                                   aTTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTAt                                                                                             >  1:374045/1‑140 (MQ=255)
                                                                      aaTCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGAGACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGa                                                                                <  2:2190457/150‑1 (MQ=255)
                                                                         cATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGcc                                                                                                                                     <  2:2150224/94‑1 (MQ=255)
                                                                         cATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGcc                                                                                                                                     >  1:2150224/1‑94 (MQ=255)
                                                                           tatTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGa                                                                                                                                   <  2:1752666/94‑1 (MQ=255)
                                                                           tatTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGa                                                                                                                                   >  1:1752666/1‑94 (MQ=255)
                                                                             tttGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGc                                                                                                                                      <  1:1020417/89‑1 (MQ=255)
                                                                             tttGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGc                                                                                                                                      >  2:1020417/1‑89 (MQ=255)
                                                                                      ccAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATc                                                                                            >  2:765805/1‑122 (MQ=255)
                                                                                      ccAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATc                                                                                            <  1:765805/122‑1 (MQ=255)
                                                                                        aGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAAt                                                              >  1:1479603/1‑150 (MQ=255)
                                                                                          gcgcATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATtgtg                                                                                                                                        >  1:33162/1‑74 (MQ=255)
                                                                                          gcgcATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATtgtg                                                                                                                                        <  2:33162/74‑1 (MQ=255)
                                                                                               ttCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCATCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCaaa                                                       <  2:1904141/150‑1 (MQ=255)
                                                                                                  aCGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCATCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGgcgc                                                                                   >  2:604158/1‑119 (MQ=255)
                                                                                                  aCGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCATCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGgcgc                                                                                   <  1:604158/119‑1 (MQ=255)
                                                                                                      aCGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACCACg                                                <  2:1451573/150‑1 (MQ=255)
                                                                                                         aGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACTACgggg                                             <  1:818033/150‑1 (MQ=255)
                                                                                                          ggTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACCACGGGGc                                            <  1:1519786/150‑1 (MQ=255)
                                                                                                            tCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACCACGGGGCTg                                          >  2:2271901/1‑150 (MQ=255)
                                                                                                              aTGGATATTAGTGCGATCGTCGAATCAAGACACGACCATCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACCACg                                                >  1:13250/1‑142 (MQ=255)
                                                                                                              aTGGATATTAGTGCGATCGTCGAATCAAGACACGACCATCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACCACg                                                <  2:13250/142‑1 (MQ=255)
                                                                                                               gtgATATGAGTGCGATCGGCGAATCATGACACGACCACCCAGCGCTGATTTTCGCCTACAACACGTGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACCACGGGGCTGTag                                       <  2:1058731/148‑1 (MQ=255)
                                                                                                                ggATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACCACGGGGCTGTaga                                      <  2:1479603/150‑1 (MQ=255)
                                                                                                                 gATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACCACGGGGCTGTagag                                     >  2:1426537/1‑150 (MQ=255)
                                                                                                                  atatTAGTGCGATCGTCGAATCAAGACACGACCATCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACCACGGGGCTGTAGAgg                                    >  2:1681464/1‑150 (MQ=255)
                                                                                                                      tAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACCACGGGGCTGTAGAGGGCAt                                <  2:1969850/150‑1 (MQ=255)
                                                                                                                        gTGCGATCGTGGAATCAAGACACGACCACCCAGCGCTGATTGTGGCGGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACCACGGGGCTGTAGAGGGCAtct                              >  2:1064136/1‑150 (MQ=255)
                                                                                                                            gATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATggg                                                                                                                        <  2:462115/56‑1 (MQ=255)
                                                                                                                            gATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATggg                                                                                                                        >  1:462115/1‑56 (MQ=255)
                                                                                                                            gATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATAc                                                    >  2:1358514/1‑124 (MQ=255)
                                                                                                                            gATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATAc                                                    <  1:1358514/124‑1 (MQ=255)
                                                                                                                                gtcgAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACCACGGGGCTg                                          >  2:1596841/1‑130 (MQ=255)
                                                                                                                                gtcgAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACCACGGGGCTg                                          <  1:1596841/130‑1 (MQ=255)
                                                                                                                                 tcgAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACCACGGGGCTGTAGAGGGCATCTCTCTCGGAg                     >  2:949986/1‑150 (MQ=255)
                                                                                                                                    aaTCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACCACGGGGCTGTAGAGGGCATCTCTCTCGGAGCAc                  >  1:2218055/1‑150 (MQ=255)
                                                                                                                                    aaTCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACCACGGGGCTGTAGAGGGCATCTCTCTCGGAGCAc                  <  2:2011383/150‑1 (MQ=255)
                                                                                                                                      tCAAGACACGACCATCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCATTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACCACGGGGCTGTAGAGGGCATCTCTCTCGGAGCACta                >  1:982145/1‑150 (MQ=255)
                                                                                                                                          gACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAAt                                                              >  2:1609211/1‑100 (MQ=255)
                                                                                                                                          gACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAAt                                                              <  1:1609211/100‑1 (MQ=255)
                                                                                                                                           acacGACCATCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACCACGGGGCTGTAGAGGGCATCTCTCTCGGAGCACTAATAAc           >  2:545082/1‑150 (MQ=255)
                                                                                                                                           acacGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACCACGGGGCTGTAGAGGGCATCTCTCTCGGAGCACTAATAAc           >  1:1826566/1‑150 (MQ=255)
                                                                                                                                           acacGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGACAAAGTACCACGGGGCTGTAGAGGGCATCTCTCTCGGAGCACTAATAAc           >  1:1421742/1‑150 (MQ=255)
                                                                                                                                             acGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACCACGGGGCTGTAGAGGGCATCTCTCTCGGAGCACTAATAACCa         >  1:1978030/1‑150 (MQ=255)
                                                                                                                                                accaccCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACCACGGGGCTGTAGAGGGCATCTCTCTCGGAGCACTAATAACCAAGc      >  1:2146651/1‑150 (MQ=255)
                                                                                                                                                    cccAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACCACGGGGCTGTAGAGGGCATCTCTCTCGGAGCACTAATAACCAAGCGaaa  >  2:1953004/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TAAAAATTGAATACTTCACACCTAGCGAAGCCGCGACTCTAAAAAATAAACTAACTCTAAAGACAAAATTAATCATATTTGAAAACCCAGGCGCATTCACGTACGAGGTCATGGATATTAGTGCGATCGTCGAATCAAGACACGACCACCCAGCGCTGATTGTGGCCGACAACACATGGGCAGCATCACTATTTCACCACCCCTTATCACTCGGCGCGGACATCAGCTTGATTTCAATGAGCAAATACCACGGGGCTGTAGAGGGCATCTCTCTCGGAGCACTAATAACCAAGCGAAA  >  NC_002947/4940907‑4941204

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