Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 531,584 C→T 7.8% R213Q (CGA→CAA)  allS ← DNA‑binding transcriptional activator AllS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913531,5840CT 7.8% 271.8 / 9.3 89R213Q (CGA→CAA) allSDNA‑binding transcriptional activator AllS
Reads supporting (aligned to +/- strand):  ref base C (39/43);  new base T (2/5);  total (41/48)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

GTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCC  >  NC_000913/531437‑531729
                                                                                                                                                   |                                                                                                                                                 
gTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGc                                                                                                                                                 <  2:358743/150‑1 (MQ=255)
gTTGGGATTACCCGGCTGACCAGATGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGc                                                                                                                                                 >  2:185266/1‑150 (MQ=255)
 ttGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCGATTGcc                                                                                                                                                <  1:1047580/150‑1 (MQ=255)
 ttGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGcc                                                                                                                                                >  1:186388/1‑150 (MQ=255)
      aTTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCGATTGCCAgggc                                                                                                                                           <  2:1260002/150‑3 (MQ=255)
      aTTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCGATTGCCAgggc                                                                                                                                           <  1:941320/150‑3 (MQ=255)
       ttACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCGATTGCCAgggca                                                                                                                                          <  2:2588/150‑4 (MQ=255)
        tACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAACCCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGAc                                                                                                                                         >  2:1362466/1‑150 (MQ=255)
         aCCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGAcg                                                                                                                                        <  1:649886/150‑1 (MQ=255)
         aCCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGAcg                                                                                                                                        >  2:540657/1‑150 (MQ=255)
         aCCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGAcg                                                                                                                                        <  2:406007/150‑1 (MQ=255)
          cccGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGAcgc                                                                                                                                       >  1:1311341/1‑150 (MQ=255)
          cccGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGAcgc                                                                                                                                       >  2:1346494/1‑150 (MQ=255)
             ggCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGtt                                                                                                                                    >  1:880828/1‑150 (MQ=255)
              gCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGttt                                                                                                                                   <  2:102775/150‑1 (MQ=255)
               cTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTg                                                                                                                                  <  2:1038475/150‑1 (MQ=255)
                tGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTgg                                                                                                                                 <  2:1297536/150‑1 (MQ=255)
                tGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTgg                                                                                                                                 <  2:1324747/150‑1 (MQ=255)
                 gACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGt                                                                                                                                >  2:177384/1‑150 (MQ=255)
                  aCCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGtt                                                                                                                               <  2:497328/150‑1 (MQ=255)
                      gttgttGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAgg                                                                                                                           >  1:476321/1‑150 (MQ=255)
                      gttgttGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAgg                                                                                                                           >  1:481560/1‑150 (MQ=255)
                      gttgttGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAgg                                                                                                                           <  2:501307/150‑1 (MQ=255)
                        tgttgATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCATATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTg                                                                                                                         >  2:741356/1‑150 (MQ=255)
                         gttgATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGc                                                                                                                        >  2:1301311/1‑150 (MQ=255)
                           tgATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCgg                                                                                                                      >  1:127423/1‑150 (MQ=255)
                            gATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCggg                                                                                                                     <  2:176033/150‑1 (MQ=255)
                                 tCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTg                                                                                                                <  2:309819/150‑1 (MQ=255)
                                 tCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTg                                                                                                                <  2:751669/150‑1 (MQ=255)
                                   gATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTc                                                                                                              >  2:364285/1‑150 (MQ=255)
                                    aTCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCt                                                                                                             <  1:736339/150‑1 (MQ=255)
                                      cATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTc                                                                                                           <  2:583734/150‑1 (MQ=255)
                                      cATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTc                                                                                                           <  2:254815/150‑1 (MQ=255)
                                         tGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAat                                                                                                        <  1:1198921/150‑1 (MQ=255)
                                          gACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAata                                                                                                       >  1:1456051/1‑150 (MQ=255)
                                           aCTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAatat                                                                                                      <  2:73051/150‑1 (MQ=255)
                                               gCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGAc                                                                                                  >  2:497067/1‑150 (MQ=255)
                                                cAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGAcc                                                                                                 >  1:882027/1‑150 (MQ=255)
                                                           tGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATTGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAgcg                                                                                      <  1:1006598/150‑1 (MQ=255)
                                                           tGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATTGCCAGGCGACGCGGTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAgcg                                                                                      <  1:832693/150‑1 (MQ=255)
                                                             gCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGTCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGcgcc                                                                                    >  2:1443603/1‑150 (MQ=255)
                                                             gCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTg                                                                                                    >  2:191284/1‑134 (MQ=255)
                                                             gCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTg                                                                                                    <  1:191284/134‑1 (MQ=255)
                                                             gCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGcgcc                                                                                    >  1:1232228/1‑150 (MQ=255)
                                                              cAAAAACCCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGCCGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGCCCGCCGGAAAGcgccg                                                                                   >  1:1041651/1‑150 (MQ=255)
                                                              cAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGcgccg                                                                                   >  1:1313589/1‑150 (MQ=255)
                                                              cAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGcgccg                                                                                   >  2:281648/1‑150 (MQ=255)
                                                              cAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACAGCCGGAAAGcgccg                                                                                   >  1:1016800/1‑150 (MQ=255)
                                                                    aCCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGCTTCCATCGCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGGGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTg                                                                             <  1:624854/150‑1 (MQ=255)
                                                                        aTGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCt                                                                         <  2:1048343/150‑1 (MQ=255)
                                                                         tGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCtt                                                                        <  2:253827/150‑1 (MQ=255)
                                                                            cAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTg                                                                     <  2:348232/150‑1 (MQ=255)
                                                                            cAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTg                                                                     >  1:708601/1‑150 (MQ=255)
                                                                               cgcccgccAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTa                                                                  <  2:131068/150‑1 (MQ=255)
                                                                                 cccgccAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCCGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGc                                                                <  1:221711/150‑1 (MQ=255)
                                                                                  ccgccAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCg                                                               >  2:1366194/1‑150 (MQ=255)
                                                                                    gccAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGc                                                             >  2:1524139/1‑150 (MQ=255)
                                                                                       aGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGctct                                                          <  1:497067/150‑1 (MQ=255)
                                                                                           gggCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCaa                                                      <  1:603174/150‑1 (MQ=255)
                                                                                            ggcggcGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTTCGCTTCTGTTAGCGGCTCTTCAAc                                                     <  1:741356/150‑1 (MQ=255)
                                                                                              cggcgATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGt                                                   >  1:405817/1‑150 (MQ=255)
                                                                                               ggcgATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGtt                                                  >  2:454023/1‑150 (MQ=255)
                                                                                                  gATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGc                                               <  1:1366194/150‑1 (MQ=255)
                                                                                                   aTTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAATGTGCGGGCGCGGTCTTCAATATTGACCGCCGGAAAGCGCAGCAACTGCGATTCTGTTAGCGGCTCTTCAACGTTCGcc                                              <  2:267910/150‑1 (MQ=255)
                                                                                                   aTTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCACTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGcc                                              <  2:186388/150‑1 (MQ=255)
                                                                                                    ttttCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCa                                             >  2:833948/1‑150 (MQ=255)
                                                                                                            tCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGgatga                                     >  1:1302012/1‑150 (MQ=255)
                                                                                                            tCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGgatga                                     >  2:482505/1‑150 (MQ=255)
                                                                                                            tCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGgatga                                     >  1:840122/1‑150 (MQ=255)
                                                                                                              cATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATc                                   >  2:1035487/1‑150 (MQ=255)
                                                                                                               atatCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGt                                                                    <  1:1467694/116‑1 (MQ=255)
                                                                                                               atatCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGt                                                                    >  2:1467694/1‑116 (MQ=255)
                                                                                                               atatCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATcc                                  <  2:1456051/150‑1 (MQ=255)
                                                                                                                tatCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCCACGTTCGCCAGCGGATGCTccg                                 >  2:1022391/1‑150 (MQ=255)
                                                                                                                tatCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATccg                                 <  2:73146/150‑1 (MQ=255)
                                                                                                                        aCAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATg                         <  2:1311341/150‑1 (MQ=255)
                                                                                                                          aataatCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGAc                       <  2:882027/150‑1 (MQ=255)
                                                                                                                          aataatCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGAc                       <  1:1301311/150‑1 (MQ=255)
                                                                                                                           ataatCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACa                      <  2:1232228/150‑1 (MQ=255)
                                                                                                                               tctcTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAgc                  <  2:1061537/150‑1 (MQ=255)
                                                                                                                                ctctTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAgcg                 >  2:1094852/1‑150 (MQ=255)
                                                                                                                                   tttttGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCa              <  1:405910/150‑1 (MQ=255)
                                                                                                                                      ttGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTg           <  1:1094852/150‑1 (MQ=255)
                                                                                                                                      gtcgccTGGCGATTGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTg           >  1:791110/5‑150 (MQ=255)
                                                                                                                                       tcgccTGGCGATTGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGc          >  1:73746/4‑150 (MQ=255)
                                                                                                                                          ccTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAACGCGCCATTGCAcc       >  1:632715/1‑150 (MQ=255)
                                                                                                                                           cTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCg      <  1:177384/150‑1 (MQ=255)
                                                                                                                                              gCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATc   >  1:970010/1‑150 (MQ=255)
                                                                                                                                               cAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATcc  <  2:842515/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GTTGGGATTACCCGGCTGACCAGTTGTTGATTATCGATCATTGACTGGCAAAGCGATTTTGGCAAAAAACCAATGCCAACGCCCGCCAGATGGGCGGCGATTTTCGTTTCCATATCGGGAACAATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCC  >  NC_000913/531437‑531729

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