Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,096,825 T→C F71S (TTC→TCC)  hisI → phosphoribosyl‑AMP cyclohydrolase and phosphoribosyl‑ATP pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,096,8250TC96.5% 295.1 / 1.1 85F71S (TTC→TCC) hisIphosphoribosyl‑AMP cyclohydrolase and phosphoribosyl‑ATP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base C (43/39);  total (46/39)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.70e-01

CAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTTCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGACACCGCTCACCAGTGGCTGTTCCTGTATCAACTGGAAC  >  NC_000913/2096683‑2096971
                                                                                                                                              |                                                                                                                                                  
cAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACg                                                                                                                                            <  2:1311383/151‑1 (MQ=255)
 aaCACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGt                                                                                                                                           <  2:522762/151‑1 (MQ=255)
     cgcgGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAg                                                                                                                                         >  1:796962/1‑149 (MQ=255)
     cgcgGCATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGt                                                                                                                                        >  1:338274/1‑150 (MQ=255)
       cgGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAg                                                                                                                                     <  1:135834/151‑1 (MQ=255)
        ggTATCCGGCGAAGCGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAg                                                                                                                                     >  1:631896/1‑150 (MQ=255)
         gTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGt                                                                                                                                    >  2:31020/1‑150 (MQ=255)
         gTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGt                                                                                                                                    >  1:230134/1‑150 (MQ=255)
         gTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTa                                                                                                                                   <  2:294409/151‑1 (MQ=255)
             ccGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTg                                                                                                                                <  1:753823/150‑1 (MQ=255)
             ccGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGc                                                                                                                               <  1:1192324/151‑1 (MQ=255)
             ccGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGc                                                                                                                               <  2:748987/151‑1 (MQ=255)
                 cGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCg                                                                                                                           <  1:1173924/151‑1 (MQ=255)
                    aGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGAc                                                                                                                        <  1:813236/151‑1 (MQ=255)
                       gCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTg                                                                                                                      <  2:362408/150‑1 (MQ=255)
                            tGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACa                                                                                                                 >  2:1419091/1‑150 (MQ=255)
                             gCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAAc                                                                                                               >  1:1398373/1‑151 (MQ=255)
                                gggCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGac                                                                                                            <  1:74734/151‑1 (MQ=255)
                                  gCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGacac                                                                                                          <  1:1125631/151‑1 (MQ=255)
                                  gCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGacac                                                                                                          >  1:745705/1‑151 (MQ=255)
                                     atatGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACg                                                                                                         <  2:1141105/149‑1 (MQ=255)
                                            cccGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTg                                                                                                <  2:654398/151‑1 (MQ=255)
                                                   gCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCg                                                                                          >  1:289562/1‑150 (MQ=255)
                                                   gCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGa                                                                                         <  2:1117804/151‑1 (MQ=255)
                                                    ccTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTACGACAACGACACGTTACTGGTGCTGGCGa                                                                                         >  2:798039/1‑150 (MQ=255)
                                                       tAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAAt                                                                                       <  2:1119036/149‑1 (MQ=255)
                                                       tAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATcc                                                                                     >  1:1023933/1‑151 (MQ=255)
                                                        aGACAAAATCCTCGAAAGCGGTATAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAATGCGAAACGTCGGGCAATTCCCTCAACGTAGAGAGTATTGACCAGGACTGCGAAAACGACATGTTACTGGAGCTGGCGAATccc                                                                                    >  1:1366483/1‑151 (MQ=255)
                                                           cAAAACCCGCGAAAGCGGCAAAGTCACCTACTTCTCGCGCACTAAACAGCGACTGTCGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATc                                                                                 <  1:145801/151‑1 (MQ=255)
                                                            aaaaCCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCACTCTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATc                                                                                 <  2:1181563/150‑1 (MQ=255)
                                                               aCCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGcc                                                                             <  1:1387222/151‑1 (MQ=255)
                                                                cccTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGcc                                                                             <  1:116084/150‑1 (MQ=255)
                                                                 ccTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCg                                                                           >  2:1357647/1‑151 (MQ=255)
                                                                    cGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGa                                                                          <  1:1420197/149‑1 (MQ=255)
                                                                    cGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGAcc                                                                        >  2:632798/1‑151 (MQ=255)
                                                                     gAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGAAGCGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGAc                                                                         >  1:318744/1‑149 (MQ=255)
                                                                       aaGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTTCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGc                                                                     >  2:1303267/1‑151 (MQ=255)
                                                                        aGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTTCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGcc                                                                    >  2:348564/1‑151 (MQ=255)
                                                                        aGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGcc                                                                    >  1:1403327/1‑151 (MQ=255)
                                                                        aGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGcc                                                                    >  2:1283409/1‑151 (MQ=255)
                                                                          cGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTTCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCcac                                                                  >  1:440334/1‑151 (MQ=255)
                                                                          cGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGgcccgccc                                                                        >  2:205512/1‑145 (MQ=255)
                                                                           ggCAAAGTCCCCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCcaca                                                                 <  2:1191151/151‑1 (MQ=255)
                                                                           ggCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCcac                                                                  >  1:283083/1‑150 (MQ=255)
                                                                           ggCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCcac                                                                  >  2:124700/1‑150 (MQ=255)
                                                                             cAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACaa                                                                >  1:107670/1‑150 (MQ=255)
                                                                              aaaGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCcaca                                                                 >  2:952613/1‑148 (MQ=255)
                                                                              aaaGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCcac                                                                  >  1:539396/1‑147 (MQ=255)
                                                                               aaGTCCCCTTCTTCTCGCCCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTAAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGATTCCCATCGGCCCGACCTGCCACAAAgg                                                             <  1:905338/151‑1 (MQ=255)
                                                                                aGTCACCTTCTTCTCGCGCACAAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAg                                                              <  2:905767/149‑1 (MQ=255)
                                                                                 gTCACCTTCTTCTCGCGCACTAAACAGCGCCTGGGGCCCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCa                                                           <  1:1390885/151‑1 (MQ=255)
                                                                                 gTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGc                                                            >  1:1266631/1‑150 (MQ=255)
                                                                                    aCCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCAcc                                                         >  2:10459/1‑150 (MQ=255)
                                                                                      cttcttCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACcagc                                                      >  1:12602/1‑151 (MQ=255)
                                                                                         cttctCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGt                                                   >  2:1163025/1‑151 (MQ=255)
                                                                                         cttctCGCGCACTAAACAGCGACTGTCGACCAAAGGCGAAACGTCGGCCAATTCCCTCAACGTAGTGAGGATTGCCCCGGACTGCGACAACGACCCGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACcagcag                                                    >  1:56864/1‑150 (MQ=255)
                                                                                          ttctCGCGCACAAAACAGCGACTGTGGGCCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATGGGCCCGACCTGCCACAAAGGCACcagc                                                      <  1:519952/147‑1 (MQ=255)
                                                                                            ctcGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGtt                                                  <  2:785145/149‑1 (MQ=255)
                                                                                                cgcACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTTCCACAAAGGCACCAGCAGCTGCtt                                              >  1:240240/1‑149 (MQ=255)
                                                                                                 gcACTAAACAGCGATTGCGGCCCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTTGTAAGTATTCCTCCGGACCGCCACAACGACACTTTTCGGGTGTTGGTGAATCCCAACGACCAACCCTGCCACAAAGGCGCCAGCAGTTGCtac                                             <  1:554668/149‑3 (MQ=255)
                                                                                                  cACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCg                                            <  2:1157063/149‑1 (MQ=255)
                                                                                                   aCTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGc                                          >  2:1268139/1‑150 (MQ=255)
                                                                                                    cTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACTTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGc                                          <  1:258124/149‑1 (MQ=255)
                                                                                                      aaccAGCGACTGTGGACCAAAGTCAAAACGTCGGGCAATCCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAAGGACACGTTCCTGGTGCTGGCGAATCCCATCGGCCCAACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCg                                         <  1:840588/145‑1 (MQ=255)
                                                                                                       aaCAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGacac                                     >  2:3218/1‑151 (MQ=255)
                                                                                                        aCAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGACAcc                                    <  2:296610/151‑1 (MQ=255)
                                                                                                        aCAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGACAcc                                    <  1:1159266/151‑1 (MQ=255)
                                                                                                             gACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGACACCGCTc                                >  1:988134/1‑150 (MQ=255)
                                                                                                             gACTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGACACCGCTCa                               <  2:382974/151‑1 (MQ=255)
                                                                                                              aCTGTGGACCAAAGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGACACCGCCCAc                              >  2:281149/1‑151 (MQ=255)
                                                                                                                   ggACCAAAGGCGAAACGTCGGGCAATCCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACCCGTTACTGGTGCTGGCGAACCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGACACCGCTCACCAGTg                         <  1:1113399/151‑1 (MQ=255)
                                                                                                                   ggACCAAAGGCGAAACGACGGGCAATTCCCTCATCGTAGTAAGTATTGCCCCGGACTGCGACACCGACACGTTACTGGTGCTGGCGAATCCCATCGCCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGACACCGCTCACCAGt                          <  2:296009/150‑1 (MQ=255)
                                                                                                                        aaaGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGACACCGCTCACCAGTGGCTg                     >  1:751365/1‑150 (MQ=255)
                                                                                                                         aaGGCGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGACACCGCTCACCAGTGGCTGtt                   >  1:25731/1‑151 (MQ=255)
                                                                                                                             cGAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGACACCGCTCACCAGTGGCTGTTCCt                >  1:44153/1‑150 (MQ=255)
                                                                                                                              gAAACGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGACACCGCTCACCAGTGGCTGTTCCTGt              <  2:456967/151‑1 (MQ=255)
                                                                                                                              gAAACGTCGGGCAATTCCCTCAACGGAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCAACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGAGACCGCTCACCAGTGGCTGTTCCTGt              >  2:826622/1‑151 (MQ=255)
                                                                                                                                  cGTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGACACCGCTCACCAGTGGCTGTTCCTGTATc           >  1:952961/1‑150 (MQ=255)
                                                                                                                                   gTCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGACACCGCTCACCAGTGGCTGTTCCTGTATc           >  2:1285138/1‑149 (MQ=255)
                                                                                                                                    tCGGGCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGATTCGGCGACACCGCTCACCAGTGGCTGTTCCTGTATCaa         >  2:1088383/1‑150 (MQ=255)
                                                                                                                                      gggCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGCCCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGACACCGCTCACCAGTGGCTGTTCCTGTATCAACTg      <  1:811107/151‑1 (MQ=255)
                                                                                                                                      gggCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGACACCGCTCACCAGTGGCTGTTCCTGTATCAACTg      >  2:226610/1‑151 (MQ=255)
                                                                                                                                        gCAATTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGACACCGCTCACCAGTGGCTGTTCCTGTATCAACTgg     >  1:806850/1‑150 (MQ=255)
                                                                                                                                          aaTTCCCTCAACGTAGTGAGTTTTGCCCCGGCCTGCGACAACGACTCGTTACTGGTGCTGGGGAATCCCATCGGCCCATCCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGACACCGCTCACCAGTGGCTGTTCCTGTATCAACTGGaac  <  2:853312/151‑1 (MQ=255)
                                                                                                                                           aTTCCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGACACCGCTCACCAGTGGCTGTTCCTGTATCAACTGGaac  >  2:646125/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
CAACACGCGGTATCCGGCGAAGTGCTAATGCTGGGCTATATGAACCCGGAAGCCTTAGACAAAACCCTCGAAAGCGGCAAAGTCACCTTCTTCTCGCGCACTAAACAGCGACTGTGGACCAAAGGCGAAACGTCGGGCAATTTCCTCAACGTAGTGAGTATTGCCCCGGACTGCGACAACGACACGTTACTGGTGCTGGCGAATCCCATCGGCCCGACCTGCCACAAAGGCACCAGCAGTTGCTTCGGCGACACCGCTCACCAGTGGCTGTTCCTGTATCAACTGGAAC  >  NC_000913/2096683‑2096971

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