Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 52,787 C→T 5.1% P210P (CCG→CCA)  pdxA ← 4‑hydroxy‑L‑threonine phosphate dehydrogenase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091352,7870CT 5.1% 254.1 / 3.2 78P210P (CCG→CCA) pdxA4‑hydroxy‑L‑threonine phosphate dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base C (40/34);  new base T (2/2);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.99e-01

GTTATCAAGATATTTCGGCTGAAACAGGGTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCA  >  NC_000913/52646‑52932
                                                                                                                                             |                                                                                                                                                 
gtTATCAAGATATTTCGGCTGAAACAGGGTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTc                                                                                                                                             <  2:895824/148‑1 (MQ=255)
    tCAAGATATTTCGGCTGAAACAGGGTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGTGGATTCAGcc                                                                                                                                         >  2:663186/1‑148 (MQ=255)
    tCAAGATATTTCGGCTGAAACAGGGTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGTATTCAGccc                                                                                                                                        <  4:486488/149‑1 (MQ=255)
    tCAAGATATTTCGGCTGAAACAGGGTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGccc                                                                                                                                        >  3:512612/1‑149 (MQ=255)
         atatTTCGGCTGAAACAGGGTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGa                                                                                                                                   >  4:331004/1‑149 (MQ=255)
           atTTCGGCTGAAACAGGGTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGAcc                                                                                                                                 >  4:596750/1‑149 (MQ=255)
           atTTCGGCTGAAACAGGGTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGAc                                                                                                                                  >  2:840308/1‑148 (MQ=255)
                   tGAAACAGGGTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATgc                                                                                                                          >  1:240963/1‑148 (MQ=255)
                       aCAGGGTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTc                                                                                                                                             >  2:105728/1‑125 (MQ=255)
                       aCAGGGTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTc                                                                                                                                             <  1:105728/125‑1 (MQ=255)
                        cAGGGTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGtt                                                                                                                    >  3:596363/1‑149 (MQ=255)
                        cAGGGTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGtt                                                                                                                    >  3:852307/1‑149 (MQ=255)
                         aGGGTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTc                                                                                                                   <  3:476452/149‑1 (MQ=255)
                         aGGGTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTc                                                                                                                   >  4:621944/1‑149 (MQ=255)
                          gggTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGAc                                                                                                                                  >  3:836860/1‑133 (MQ=255)
                          gggTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGAc                                                                                                                                  <  4:836860/133‑1 (MQ=255)
                           ggTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCg                                                                                                                  <  1:145262/148‑1 (MQ=255)
                              aTCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGTGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCa                                                                                                               <  1:663186/148‑1 (MQ=255)
                              aTCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCa                                                                                                               <  2:732076/148‑1 (MQ=255)
                                     ggcagCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCaa                                                                                                        <  2:285149/148‑1 (MQ=255)
                                     ggcagCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCaa                                                                                                        <  2:519824/148‑1 (MQ=255)
                                     ggcagCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCaa                                                                                                        <  1:302581/148‑1 (MQ=255)
                                         gCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTg                                                                                                   >  3:684001/1‑149 (MQ=255)
                                            gCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGTGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATa                                                                                                            <  4:301876/137‑1 (MQ=255)
                                            gCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGTGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATa                                                                                                            >  3:301876/1‑137 (MQ=255)
                                            gCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTc                                                                                                <  3:393027/149‑1 (MQ=255)
                                             cccGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTgg                                                                                                  >  2:365678/1‑146 (MQ=255)
                                             cccGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTgg                                                                                                  <  1:365678/146‑1 (MQ=255)
                                               cGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCg                                                                                              >  1:696628/1‑148 (MQ=255)
                                                  tGAGTTTCATCCCCTGCGCGCGGATCTCATTGAGCACCGGAATAATGGTGTTTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGTAGACCAGAATGCTCGCTTCGGCAATACCAAATTTGGTCCGCaaa                                                                                          >  4:311612/1‑149 (MQ=255)
                                                  tGAGTATCACCCCCTGCGCCCGCAGCTCATTGAGCAGCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGCCAATACCAAATTTGGTCCGCaa                                                                                           >  2:879831/1‑148 (MQ=255)
                                                      tttCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGt                                                                                      >  4:102493/1‑149 (MQ=255)
                                                         cATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGa                                                                                    <  2:571052/148‑1 (MQ=255)
                                                              ccTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCaaa                                                                              >  3:434880/1‑149 (MQ=255)
                                                               cTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCaaa                                                                              >  1:624344/1‑148 (MQ=255)
                                                                 gcgcCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGaa                                                                                                                             <  1:161514/99‑1 (MQ=255)
                                                                 gcgcCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGaa                                                                                                                             >  2:161514/1‑99 (MQ=255)
                                                                   gcCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAg                                                                          <  2:624344/148‑1 (MQ=255)
                                                                        cAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATc                                                                     >  1:281055/1‑148 (MQ=255)
                                                                              aTTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCg                                                               <  1:490410/148‑1 (MQ=255)
                                                                              aTTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGt                                                              <  3:596750/149‑1 (MQ=255)
                                                                               ttGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTg                                                             <  4:611547/149‑1 (MQ=255)
                                                                                    cACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCaaa                                                         <  2:388261/148‑1 (MQ=255)
                                                                                       cGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGt                                                      >  2:629638/1‑148 (MQ=255)
                                                                                       cGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCGTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTg                                                     >  3:562590/1‑149 (MQ=255)
                                                                                        ggAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGt                                                                                                 >  1:303742/1‑104 (MQ=255)
                                                                                        ggAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGt                                                                                                 <  2:303742/104‑1 (MQ=255)
                                                                                        ggAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTc                                                                                                >  2:53360/1‑105 (MQ=255)
                                                                                        ggAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTc                                                                                                <  1:53360/105‑1 (MQ=255)
                                                                                        ggAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTg                                                     >  1:703560/1‑148 (MQ=255)
                                                                                            taatGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAgg                                                 <  1:840308/148‑1 (MQ=255)
                                                                                            taatGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAgg                                                 <  2:12200/148‑1 (MQ=255)
                                                                                              atGGTGTCTATCTGTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGt                                              <  4:684170/149‑1 (MQ=255)
                                                                                                ggTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTg                                             <  2:703560/148‑1 (MQ=255)
                                                                                                  tgtCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGAt                                           >  1:640564/1‑148 (MQ=255)
                                                                                                     cTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGt                                                                                      <  1:853116/102‑1 (MQ=255)
                                                                                                     cTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGt                                                                                      >  2:853116/1‑102 (MQ=255)
                                                                                                      tATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCg                                       >  2:864664/1‑148 (MQ=255)
                                                                                                       aTCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTc                                     >  4:748257/1‑149 (MQ=255)
                                                                                                          tctTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGt                                                              <  3:230410/121‑1 (MQ=255)
                                                                                                          tctTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGt                                                              >  4:230410/1‑121 (MQ=255)
                                                                                                          tctTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGc                                  >  4:645107/1‑149 (MQ=255)
                                                                                                          tctTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGc                                  >  3:103434/1‑149 (MQ=255)
                                                                                                          tctTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGc                                  >  4:798993/1‑149 (MQ=255)
                                                                                                            ttCCGGATCCATATGACGGCCTTCGCCTGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCg                                 >  2:655996/1‑148 (MQ=255)
                                                                                                                gTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATc                            >  3:545985/1‑149 (MQ=255)
                                                                                                                      atatGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCaaaa                                                                             <  4:893580/94‑1 (MQ=255)
                                                                                                                      atatGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCaaaa                                                                             >  3:893580/1‑94 (MQ=255)
                                                                                                                       tatGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGc                     <  3:737051/149‑1 (MQ=255)
                                                                                                                                cTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGTGTCTGCGATATCGCGCAGCGGTAAATg             >  2:128609/1‑148 (MQ=255)
                                                                                                                                cTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATg             >  2:872304/1‑148 (MQ=255)
                                                                                                                                   cGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGcgt          <  2:90964/148‑1 (MQ=255)
                                                                                                                                   cGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGcgt          <  2:240963/148‑1 (MQ=255)
                                                                                                                                      ccGCGTGCGTATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGt                                      <  3:817384/117‑1 (MQ=255)
                                                                                                                                      ccGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGt                                      >  4:817384/1‑117 (MQ=255)
                                                                                                                                      ccGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGcgtcgt       <  2:706530/148‑1 (MQ=255)
                                                                                                                                       cgcgTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTg     >  4:777217/1‑149 (MQ=255)
                                                                                                                                          gTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCa  <  4:852307/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
GTTATCAAGATATTTCGGCTGAAACAGGGTATCGGCAGGCAGCGGCCCGTTGAGTTTCATCCCCTGCGCCCGCAGCTCATTGAGCACCGGAATAATGGTGTCTATCTCTTCCGTACCCATATGACCGCCTTCGCCCGCGTGCGGATTCAGCCCGCAGACCAGAATGCGCGGTTCGGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCA  >  NC_000913/52646‑52932

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑