Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 432,018 C→T 5.5% C100Y (TGC→TAC) ‡ PP_0354 ← CBS domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947432,0180CT 5.5% 265.3 / 2.4 90C100Y (TGC→TAC) ‡PP_0354CBS domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (41/43);  new base T (3/2);  total (45/45)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ACCGCGCAGGGCAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCACCGCAGGCGT  >  NC_002947/431870‑432162
                                                                                                                                                    |                                                                                                                                                
acCGCGCAGGGCAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCa                                                                                                                                                 >  1:1093141/1‑150 (MQ=255)
acCGCGCAGGGCAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCa                                                                                                                                                 <  2:1638266/150‑1 (MQ=255)
       aGGGCAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCg                                                                                                                                          >  2:718442/1‑150 (MQ=255)
                  ccGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGAtctc                                                                                                                               <  1:1324668/150‑1 (MQ=255)
                   cGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCg                                                                                                                              >  2:196784/1‑150 (MQ=255)
                           aTCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGt                                                                                                                      >  1:1599350/1‑150 (MQ=255)
                           aTCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACATGGt                                                                                                                      <  2:1805390/150‑1 (MQ=255)
                            tCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGtt                                                                                                                     >  1:1064132/1‑150 (MQ=255)
                            tCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGtt                                                                                                                     <  2:719770/150‑1 (MQ=255)
                              aGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTc                                                                                                                   <  1:419936/150‑1 (MQ=255)
                                aCATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGa                                                                                                                 <  1:2240201/150‑1 (MQ=255)
                                  aTCGAGAACAGGCGGTTCCCTGCGACCTTGCTCCGTTTCACGCGGAATGTGGCGTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTTCCAGCTGGAAGCCCTCCCCGCTGGTGATCTCGAACGTGGTTTCGaca                                                                                                               <  2:260318/150‑1 (MQ=255)
                                  aTCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGGAGTGCTTCCCGCTggcagcc                                                                                                                                 <  2:1309726/132‑6 (MQ=255)
                                  aTCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGaca                                                                                                               >  1:897079/1‑150 (MQ=255)
                                  aTCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGAGAGTGCTTCCCGCTggcagcc                                                                                                                                 >  1:1309726/1‑127 (MQ=255)
                                     gagaACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACAccc                                                                                                            <  1:1531637/150‑1 (MQ=255)
                                     gagaACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACAccc                                                                                                            <  2:2083892/150‑1 (MQ=255)
                                     gagaACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCAACCAGCGGGATTAGCAGTGCGATGCACACCCCGCTAGGTAGCTCTAACGTGGTTTCGactaca                                                                                                            >  2:23878/1‑146 (MQ=255)
                                                gATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAg                                                                                                 <  2:1064132/150‑1 (MQ=255)
                                                gATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAg                                                                                                 <  2:806289/150‑1 (MQ=255)
                                                       gCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGTGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTgg                                                                                          <  1:846202/150‑1 (MQ=255)
                                                       gCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTgg                                                                                          <  1:293471/150‑1 (MQ=255)
                                                         ggCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGgcg                                                                                        >  1:512450/1‑150 (MQ=255)
                                                         ggCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGgcg                                                                                        >  2:482767/1‑150 (MQ=255)
                                                          gCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGgcgc                                                                                       >  2:616924/1‑150 (MQ=255)
                                                           ccTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTg                                                                                           <  1:699097/145‑1 (MQ=255)
                                                           ccTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTg                                                                                           >  2:699097/1‑145 (MQ=255)
                                                               gTTCAGTTTCAGGCAGTAGGTGTGGTCCCCGGCCAGGTACTCGGCGCGGGTTGCCCGCTCACCGCGCAGGGCTGCCAGCGGGAAGGACTCCCCGCTGGGGATCTCGAACGTGGTTTCGACACCCTCTTAGGCCAGGTGATGGCGCTCggc                                                                                  >  2:634051/1‑148 (MQ=255)
                                                                 tCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCaa                                                                                <  1:482767/150‑1 (MQ=255)
                                                                 tCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCaa                                                                                >  2:1757053/1‑150 (MQ=255)
                                                                     ttGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTg                                                                                                                                    >  1:269669/1‑94 (MQ=255)
                                                                     ttGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTg                                                                                                                                    <  2:269669/94‑1 (MQ=255)
                                                                     ttGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCAc                                                                            <  1:1469148/150‑1 (MQ=255)
                                                                         aggcagAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCgg                                                                        >  2:2050838/1‑150 (MQ=255)
                                                                          ggcagAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGAAATGCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGc                                                                       <  1:23878/150‑1 (MQ=255)
                                                                           gcagAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGcc                                                                      >  2:601216/1‑150 (MQ=255)
                                                                           gcagAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGcc                                                                      >  2:1090544/1‑150 (MQ=255)
                                                                            cagAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGccc                                                                     >  2:2377961/1‑150 (MQ=255)
                                                                                     tCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGGAGTGCTTCCCGCTgg                                                                                                                                      <  2:466819/76‑1 (MQ=37)
                                                                                     tCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGGAGTGCTTCCCGCTgg                                                                                                                                      >  1:466819/1‑76 (MQ=37)
                                                                                     tCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACg                                                            >  1:186270/1‑150 (MQ=255)
                                                                                      cTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGa                                                           <  1:2029955/150‑1 (MQ=255)
                                                                                              gCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGAt                                                          <  2:1763477/143‑1 (MQ=255)
                                                                                              gCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGAt                                                          >  1:1763477/1‑143 (MQ=255)
                                                                                               ccAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAAACACGCGGCCCCGCTTGACGTTGTAGAAGt                                                  <  2:149615/150‑1 (MQ=255)
                                                                                                     gCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGaa                                                    >  2:1834697/1‑142 (MQ=255)
                                                                                                     gCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGaa                                                    <  1:1834697/142‑1 (MQ=255)
                                                                                                      cTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGttt                                                                                                       >  1:1029174/1‑90 (MQ=255)
                                                                                                      cTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGttt                                                                                                       >  1:1020543/1‑90 (MQ=255)
                                                                                                      cTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGttt                                                                                                       <  2:1029174/90‑1 (MQ=255)
                                                                                                      cTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGttt                                                                                                       <  2:1020543/90‑1 (MQ=255)
                                                                                                      cTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGgcg                                                                                        >  1:2177650/1‑105 (MQ=255)
                                                                                                      cTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGgcg                                                                                        <  2:2177650/105‑1 (MQ=255)
                                                                                                      cTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCaa                                                                                <  2:1658162/113‑1 (MQ=255)
                                                                                                      cTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCaa                                                                                >  1:1658162/1‑113 (MQ=255)
                                                                                                          gTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGc                                       <  1:196784/150‑1 (MQ=255)
                                                                                                            gCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGACGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGccc                                     <  1:501584/150‑1 (MQ=255)
                                                                                                            gCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGaca                                                                                                               <  1:1971090/76‑1 (MQ=255)
                                                                                                            gCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGaca                                                                                                               >  2:1971090/1‑76 (MQ=255)
                                                                                                               ggTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTTGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTc                                  <  1:61577/150‑1 (MQ=255)
                                                                                                               ggTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGc                                       >  2:925293/1‑145 (MQ=255)
                                                                                                               ggTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGc                                       <  1:925293/145‑1 (MQ=255)
                                                                                                                    cccGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAg                             <  1:1273434/150‑1 (MQ=255)
                                                                                                                     ccGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCAcgc                                                                          <  1:2309952/104‑1 (MQ=255)
                                                                                                                     ccGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCAcgc                                                                          >  2:2309952/1‑104 (MQ=255)
                                                                                                                     ccGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGc                                                                       >  1:2403164/1‑107 (MQ=255)
                                                                                                                     ccGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGc                                                                       <  2:2403164/107‑1 (MQ=255)
                                                                                                                     ccGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAgg                            >  2:813959/1‑150 (MQ=255)
                                                                                                                     ccGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAgg                            >  1:321908/1‑150 (MQ=255)
                                                                                                                     ccGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAgg                            >  2:18781/1‑150 (MQ=255)
                                                                                                                      cGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGt                           >  1:836672/1‑150 (MQ=255)
                                                                                                                       gCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTg                                                                                           <  1:380962/85‑1 (MQ=255)
                                                                                                                       gCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTg                                                                                           >  2:380962/1‑85 (MQ=255)
                                                                                                                           aCCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGAc                      >  2:819419/1‑150 (MQ=255)
                                                                                                                             ccAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGa                    >  1:1396291/1‑150 (MQ=255)
                                                                                                                             ccAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGa                    >  2:690426/1‑150 (MQ=255)
                                                                                                                              ctgcagGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGAt                   <  2:2087585/148‑1 (MQ=255)
                                                                                                                              cagcagGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGAt                   <  1:718442/150‑1 (MQ=255)
                                                                                                                                     gCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCTCGTCGGCAGGTTTGACGATGCTTTCa            <  2:1599350/150‑1 (MQ=255)
                                                                                                                                         ccAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCACCGc        >  2:892195/1‑150 (MQ=255)
                                                                                                                                         ccAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCACCGc        >  1:548100/1‑150 (MQ=255)
                                                                                                                                          cAGCGGGGAGTGCTTCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCACCGCa       >  1:932336/1‑150 (MQ=255)
                                                                                                                                          cAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCACCGCa       <  2:548100/150‑1 (MQ=255)
                                                                                                                                          cAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCACCGCa       <  1:616924/150‑1 (MQ=255)
                                                                                                                                           aGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCACCGCAg      <  1:819419/150‑1 (MQ=255)
                                                                                                                                            gCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCACCGCAgg     >  1:132327/1‑150 (MQ=255)
                                                                                                                                             cGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCACCGCAGGc    <  1:335446/150‑1 (MQ=255)
                                                                                                                                               ggAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGc                                       <  2:535864/113‑1 (MQ=255)
                                                                                                                                               ggAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGc                                       >  1:535864/1‑113 (MQ=255)
                                                                                                                                               ggAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCACCGCAGGCGt  >  1:2295545/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
ACCGCGCAGGGCAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGATGCTTTCACCGCAGGCGT  >  NC_002947/431870‑432162

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