Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,369,371 G→A D445N (GAT→AAT)  arnA → fused UDP‑L‑Ara4N formyltransferase/UDP‑GlcA C‑4'‑decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,369,3710GA97.8% 319.3 / ‑1.7 89D445N (GAT→AAT) arnAfused UDP‑L‑Ara4N formyltransferase/UDP‑GlcA C‑4'‑decarboxylase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (37/50);  total (39/50)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.28e-01

AATATACCCGCAACCCGCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCGATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCCTCTTCCGCCCGTTTAACTGG  >  NC_000913/2369228‑2369517
                                                                                                                                               |                                                                                                                                                  
aaTATACCCGCAACCCGCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGa                                                                                                                                             <  2:706005/151‑1 (MQ=255)
 atatACCCGCAACCCGCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGAc                                                                                                                                            <  1:834870/151‑1 (MQ=255)
 atatACCCGCAACCCGCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGAc                                                                                                                                            >  1:572968/1‑151 (MQ=255)
   ataCCCGCAACCCGCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGAcc                                                                                                                                           <  2:706691/150‑1 (MQ=255)
   ataCCCGCAACCCGCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCa                                                                                                                                          <  2:24024/151‑1 (MQ=255)
     aCCCGCAACCCGCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCAtt                                                                                                                                        >  2:980177/1‑151 (MQ=255)
     aCCCGCAACCCGCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCAtt                                                                                                                                        <  2:1312943/151‑1 (MQ=255)
     aCCCGCAACCCGCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCAtt                                                                                                                                        <  1:277045/151‑1 (MQ=255)
       ccGCAACCCGCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCt                                                                                                                                      <  1:325421/151‑1 (MQ=255)
        cGCAACCCGCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTa                                                                                                                                     <  1:587566/151‑1 (MQ=255)
            aCCCGCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAAttt                                                                                                                                 >  2:31545/1‑151 (MQ=255)
             cccGCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAAttt                                                                                                                                 <  1:898393/150‑1 (MQ=255)
             cccGCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTa                                                                                                                                <  1:317187/151‑1 (MQ=255)
                 cTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAAt                                                                                                                              >  1:1311665/1‑149 (MQ=255)
                   gcgcGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAAtcgtc                                                                                                                          >  2:657739/1‑151 (MQ=255)
                     gcgTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCGATGAGGACCATTCTAATTTAATCGTCgg                                                                                                                        >  2:358736/1‑151 (MQ=255)
                         aTTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGccc                                                                                                                     <  1:774112/150‑1 (MQ=255)
                            tGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGt                                                                                                                  <  2:327312/150‑1 (MQ=255)
                                   gATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAAc                                                                                                          >  1:663443/1‑151 (MQ=255)
                                       ttGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCAc                                                                                                       <  2:622857/150‑1 (MQ=255)
                                        tGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACg                                                                                                      <  2:693023/150‑1 (MQ=255)
                                        tGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACg                                                                                                      <  1:1248493/150‑1 (MQ=255)
                                        tGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACg                                                                                                      <  2:619544/150‑1 (MQ=255)
                                              gaATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGTCCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGAtt                                                                                               <  2:820221/151‑1 (MQ=255)
                                               aaTCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGAttt                                                                                              >  1:1324258/1‑151 (MQ=255)
                                               aaTCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGAttt                                                                                              <  2:381384/151‑1 (MQ=255)
                                               aaTCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGAttt                                                                                              <  2:995823/151‑1 (MQ=255)
                                               aaTCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGAatt                                                                                              >  1:1395192/1‑151 (MQ=255)
                                               aaTCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCAGTGAATAAACCACGCTGGAttt                                                                                              <  2:1061856/151‑1 (MQ=255)
                                                aTCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGAttt                                                                                              >  1:925827/1‑150 (MQ=255)
                                                  cTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGCACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTAtt                                                                                           >  1:1081569/1‑151 (MQ=255)
                                                    gctcATGATCCAGTACTGCAAGATTTACCGTAGGCGAATCGGCTTCCTGTCAACTGCTGATGTTTATGGGATGTGTAGCGATAAAGACTTCAATGAGGTCCATTGTGATTTAATCGTCGGCCTGGTTAATAAACCACGCTGGAGTTAATc                                                                                          <  1:1036901/147‑1 (MQ=255)
                                                    gcgcATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTc                                                                                          <  2:787684/150‑1 (MQ=255)
                                                      gcATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCGGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCg                                                                                         <  1:501410/149‑1 (MQ=255)
                                                      gcATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAAGCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGt                                                                                       <  2:1347481/151‑1 (MQ=255)
                                                        aTTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCGATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTAt                                                                                     >  2:805048/1‑151 (MQ=255)
                                                         ttATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATc                                                                                    >  2:1380745/1‑151 (MQ=255)
                                                           aTCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCa                                                                                   >  1:1049166/1‑150 (MQ=255)
                                                           aTCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCa                                                                                   <  1:1251886/150‑1 (MQ=255)
                                                             ccGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCaaa                                                                                 <  2:1226053/150‑1 (MQ=255)
                                                                  aCTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAAt                                                                            <  2:1065320/150‑1 (MQ=255)
                                                                    tGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTCTGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTAc                                                                         >  2:121442/1‑151 (MQ=255)
                                                                     gCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACt                                                                        >  2:340367/1‑151 (MQ=255)
                                                                      cGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACtt                                                                       <  1:955712/151‑1 (MQ=255)
                                                                      cGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACtt                                                                       <  2:1141115/151‑1 (MQ=255)
                                                                       gTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTg                                                                      <  2:3051/151‑1 (MQ=255)
                                                                        tGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTg                                                                      >  2:705371/1‑150 (MQ=255)
                                                                        tGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTg                                                                      >  2:670884/1‑150 (MQ=255)
                                                                         gAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGa                                                                     <  2:162362/150‑1 (MQ=255)
                                                                           aGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCg                                                                  <  2:354774/151‑1 (MQ=255)
                                                                                cGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTg                                                              >  2:995843/1‑150 (MQ=255)
                                                                                cGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTg                                                              >  1:1131630/1‑150 (MQ=255)
                                                                                cGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTg                                                              >  2:804503/1‑150 (MQ=255)
                                                                                cGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTg                                                              <  2:1043952/150‑1 (MQ=255)
                                                                                  tAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGAt                                                            <  1:870917/150‑1 (MQ=255)
                                                                                   aaGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCt                                                          <  2:534315/151‑1 (MQ=255)
                                                                                      cGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTggg                                                       <  2:1425707/151‑1 (MQ=255)
                                                                                       gAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGc                                                      >  2:1103501/1‑151 (MQ=255)
                                                                                       gAATCATCTTCCCGCCAACTTCAGAAGCTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGTTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTggg                                                       >  1:390250/1‑150 (MQ=255)
                                                                                          tcatcTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCt                                                    >  2:804098/1‑150 (MQ=255)
                                                                                           catcTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTAt                                                  >  1:130508/1‑151 (MQ=255)
                                                                                             tcTTCCCGTCAACTTCAGGAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGCTTTATTCGGCATCAAAACAATTACTTGATCGGGAGATCTGGGCCTATg                                                 >  1:607982/1‑150 (MQ=255)
                                                                                               ttCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCg                                              >  2:1156946/1‑151 (MQ=255)
                                                                                                tCCCGTCAACTTCAGAAGTTTATGTGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAACCCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATgg                                                <  1:588474/148‑1 (MQ=255)
                                                                                                tCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATg                                                 <  1:434446/147‑1 (MQ=255)
                                                                                                  ccGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGaaa                                           <  1:588137/151‑1 (MQ=255)
                                                                                                     tCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGaaaaa                                         <  1:838384/150‑1 (MQ=255)
                                                                                                         cTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAggg                                    >  2:1006579/1‑151 (MQ=255)
                                                                                                           tCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAggg                                    <  2:332533/149‑1 (MQ=255)
                                                                                                           tCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGt                                   <  1:678259/150‑1 (MQ=255)
                                                                                                           tCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTCGGCCTATGGCGAAAAAGAGGGt                                   <  1:408728/150‑1 (MQ=255)
                                                                                                            cAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGttt                                 <  1:264548/151‑1 (MQ=255)
                                                                                                            cAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGtt                                  <  2:982652/150‑1 (MQ=255)
                                                                                                               aaGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTAc                               >  2:1399678/1‑150 (MQ=255)
                                                                                                               aaGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACa                              <  2:987624/151‑1 (MQ=255)
                                                                                                                aGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTAc                               >  1:595991/1‑149 (MQ=255)
                                                                                                                aGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAg                             <  2:1059868/151‑1 (MQ=255)
                                                                                                                  tttATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGtt                           >  2:99488/1‑151 (MQ=255)
                                                                                                                   ttATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTc                          >  2:104014/1‑151 (MQ=255)
                                                                                                                         gATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCctc                    >  1:811111/1‑151 (MQ=255)
                                                                                                                          aTGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCctct                   >  1:902815/1‑151 (MQ=255)
                                                                                                                            gtgtAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCCTCtt                  <  1:921910/150‑1 (MQ=255)
                                                                                                                                aGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCCTCTTCCGcc             <  2:1360445/151‑1 (MQ=255)
                                                                                                                                  cGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCCTCTTCCGCCCg           >  1:840015/1‑151 (MQ=255)
                                                                                                                                    aTAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCCTCTTCCGCCCGtt         >  2:894042/1‑151 (MQ=255)
                                                                                                                                     tAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCCTCTTCCGCCCGttt        >  1:441290/1‑151 (MQ=255)
                                                                                                                                       aaTACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCCTCTTCCGCCCGTTTa       >  1:1424628/1‑150 (MQ=255)
                                                                                                                                          aCTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCCTCTTCCGCCCGTTTAACTg   <  2:392613/151‑1 (MQ=255)
                                                                                                                                           cTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCCTCTTCCGCCCGTTTAACTgg  >  1:398685/1‑151 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
AATATACCCGCAACCCGCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCGATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCCTCTTCCGCCCGTTTAACTGG  >  NC_000913/2369228‑2369517

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