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% 305.4 / ‑1.3 89D445N (GAT→AAT) arnAfused UDP‑L‑Ara4N formyltransferase/UDP‑GlcA C‑4'‑decarboxylase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (49/38);  total (49/40)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.93e-01

TGAATATACCCGCAACCCGCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCGATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCCTCTTCCGCCCGTTTAACTG  >  NC_000913/2369226‑2369516
                                                                                                                                                 |                                                                                                                                                 
tGAATATACCCGCAACCCGCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGa                                                                                                                                               <  2:242959/150‑1 (MQ=255)
  aaTATACCCGCAACCCGCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGa                                                                                                                                            <  2:733910/151‑1 (MQ=255)
       aCCCGCAACCCGCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCAtt                                                                                                                                       >  1:1157243/1‑151 (MQ=255)
        cccGCAACCCGCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCAtt                                                                                                                                       <  2:27242/150‑1 (MQ=255)
               cccGCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAAttt                                                                                                                                >  2:60979/1‑150 (MQ=255)
                ccGCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTa                                                                                                                               >  2:594656/1‑150 (MQ=255)
                 cGCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAAt                                                                                                                             >  1:834471/1‑151 (MQ=255)
                  gCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCACATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCGATGAGGACCATTCTAATTTAAtc                                                                                                                            <  1:469402/151‑1 (MQ=255)
                     gcgcGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAAtcgtc                                                                                                                         >  2:233296/1‑151 (MQ=255)
                       gcgTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCgg                                                                                                                       <  2:183150/151‑1 (MQ=255)
                       gcgTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGGACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCgg                                                                                                                       >  2:291897/1‑151 (MQ=255)
                        cgTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGc                                                                                                                      <  2:400808/151‑1 (MQ=255)
                              tGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTg                                                                                                                >  2:1103803/1‑151 (MQ=255)
                              tGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTg                                                                                                                >  1:357221/1‑151 (MQ=255)
                               gAACTTGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTg                                                                                                                >  2:1039053/1‑150 (MQ=255)
                                aaCTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGaa                                                                                                              <  1:402262/151‑1 (MQ=255)
                                 acccGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAAt                                                                                                             <  2:180505/148‑1 (MQ=255)
                                    cGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGTCCGGTGAATaaa                                                                                                          <  1:777200/151‑1 (MQ=255)
                                    cGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATaaa                                                                                                          <  1:267363/151‑1 (MQ=255)
                                    cGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATaa                                                                                                           >  2:565973/1‑150 (MQ=255)
                                     gATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAAc                                                                                                         >  2:532328/1‑151 (MQ=255)
                                     gATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAAc                                                                                                         >  2:210178/1‑151 (MQ=255)
                                       ttttGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAAcc                                                                                                        >  2:597747/1‑150 (MQ=255)
                                        tttGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCa                                                                                                       <  1:761898/150‑1 (MQ=255)
                                        tttGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCAc                                                                                                      >  1:329749/1‑151 (MQ=255)
                                        tttGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCAc                                                                                                      >  1:287288/1‑151 (MQ=255)
                                          tGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACg                                                                                                     <  2:954847/150‑1 (MQ=255)
                                           gAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCt                                                                                                   >  1:418931/1‑151 (MQ=255)
                                             agagAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTgg                                                                                                 >  2:398001/1‑151 (MQ=255)
                                             agagAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTgg                                                                                                 <  1:1206137/151‑1 (MQ=255)
                                               agaATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGAt                                                                                               <  2:397628/151‑1 (MQ=255)
                                               agaATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGAt                                                                                               <  1:175050/151‑1 (MQ=255)
                                               agaATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGAt                                                                                               <  2:1021826/151‑1 (MQ=255)
                                                      gcgcATTATCCGCTACTGCGTGAAGTACCTTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTGTATGGGATGTGTAGCGAGATATTCTACAATGAGGACCATTCTAATTTAATCGTCGAAATGGTGAATAAACCGCGCTGGATTTATTc                                                                                         <  1:876650/150‑1 (MQ=255)
                                                      gcgcATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTAt                                                                                           <  2:1207958/148‑1 (MQ=255)
                                                          aTTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTAt                                                                                    >  1:909917/1‑151 (MQ=255)
                                                            tATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATc                                                                                   >  1:886025/1‑150 (MQ=255)
                                                             aTCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCa                                                                                  <  1:1117529/150‑1 (MQ=255)
                                                                cGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAAc                                                                              <  2:19915/151‑1 (MQ=255)
                                                                cGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAAc                                                                              <  2:327440/151‑1 (MQ=255)
                                                                cGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAAc                                                                              <  1:791473/151‑1 (MQ=255)
                                                                     cTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAAtt                                                                          >  2:99954/1‑150 (MQ=255)
                                                                     cTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAAtt                                                                          >  1:662939/1‑150 (MQ=255)
                                                                        cGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACtt                                                                      >  1:109019/1‑151 (MQ=255)
                                                                         gTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTg                                                                     <  1:1140698/151‑1 (MQ=255)
                                                                         gTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTg                                                                     <  1:892472/151‑1 (MQ=255)
                                                                         gTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGATGGATTTATTCGGTATCAAAACAATTACTTg                                                                     >  1:565180/1‑151 (MQ=255)
                                                                           gAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGa                                                                    <  2:980737/150‑1 (MQ=255)
                                                                           gAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGAt                                                                   >  2:844716/1‑151 (MQ=255)
                                                                           gAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGAt                                                                   <  1:484049/151‑1 (MQ=255)
                                                                           gAAGCACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCGATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGa                                                                    <  2:1076468/150‑1 (MQ=255)
                                                                             aGTACCGTAAGCGAATCATCTTCCCGTCAACTTCTGAAGTTTATGGGATGGGTAGCGATAAATACTTCACTGATGACCATTCTAATTTAATCTTCGGCCCGGTGAATAAACCACGCTGGATTCATTCGGCATCAACACAATTACTTGATCg                                                                 >  1:1172847/1‑151 (MQ=255)
                                                                                aCCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCggg                                                               <  2:158772/150‑1 (MQ=255)
                                                                                aCCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGt                                                              <  2:978557/151‑1 (MQ=255)
                                                                                 ccGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTg                                                             <  2:316715/151‑1 (MQ=255)
                                                                                  cGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGa                                                            <  2:1104301/151‑1 (MQ=255)
                                                                                   gTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGAt                                                           >  2:61782/1‑151 (MQ=255)
                                                                                       gCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTgg                                                       <  2:672853/151‑1 (MQ=255)
                                                                                          aaTCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGcc                                                    >  1:275114/1‑151 (MQ=255)
                                                                                            tcatcTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCt                                                   >  1:814569/1‑150 (MQ=255)
                                                                                             catcTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTAt                                                 >  1:721490/1‑151 (MQ=255)
                                                                                                 ttCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCg                                             >  1:485462/1‑151 (MQ=255)
                                                                                                  tCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGc                                              <  2:287296/149‑1 (MQ=255)
                                                                                                  tCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCg                                             >  1:155289/1‑150 (MQ=255)
                                                                                                      gTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGaaaa                                         >  1:399995/1‑150 (MQ=255)
                                                                                                        cAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGaaaa                                         >  2:721430/1‑148 (MQ=255)
                                                                                                         aaCTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAag                                       >  2:977892/1‑149 (MQ=255)
                                                                                                             tCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGtt                                 >  2:828984/1‑151 (MQ=255)
                                                                                                               aGAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGttt                                >  1:471856/1‑150 (MQ=255)
                                                                                                                gAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGCGGGTTTAc                              >  1:976077/1‑151 (MQ=255)
                                                                                                                gAAGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTAc                              >  1:589425/1‑151 (MQ=255)
                                                                                                                 aaGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTAc                              >  2:1017350/1‑150 (MQ=255)
                                                                                                                  aGTTTATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTAc                              >  1:570886/1‑149 (MQ=255)
                                                                                                                      tATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTc                         <  2:1118574/150‑1 (MQ=255)
                                                                                                                      tATGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTc                         >  1:1160673/1‑150 (MQ=255)
                                                                                                                       aTGGGATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCa                        <  2:730503/150‑1 (MQ=255)
                                                                                                                         gggATGTGTAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCAccc                     >  2:994923/1‑151 (MQ=255)
                                                                                                                               tgtAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCCTCTTcc               >  1:565427/1‑151 (MQ=255)
                                                                                                                                gtAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCCTCTTcc               <  2:160691/150‑1 (MQ=255)
                                                                                                                                 tAGCGATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCCTCTTCCGc             >  2:290524/1‑151 (MQ=255)
                                                                                                                                     gATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCCTCTTCCGCCCGt         <  2:447775/151‑1 (MQ=255)
                                                                                                                                     gATAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCACGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCCTCTTCCGCCCGt         <  2:1006154/151‑1 (MQ=255)
                                                                                                                                       tAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCCTCTTCCGCCCGttt       >  1:910548/1‑151 (MQ=255)
                                                                                                                                       tAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCCTCTTCCGCCCGttt       >  1:417920/1‑151 (MQ=255)
                                                                                                                                       tAAATACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCCTCTTCCGCCCGttt       >  1:357998/1‑151 (MQ=255)
                                                                                                                                        aaaTACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCCTCTTCCGCCCGTTTa      <  1:692999/151‑1 (MQ=255)
                                                                                                                                         aaTACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCCTCTTCCGCCCGTTTaa     <  2:672543/151‑1 (MQ=255)
                                                                                                                                           tACTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCCTCTTCCGCCCGTTTAACt   >  2:467687/1‑151 (MQ=255)
                                                                                                                                            aCTTCAATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCCTCTTCCGCCCGTTTAACTg  >  1:176533/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TGAATATACCCGCAACCCGCTGCGCGTATTTGAACTCGATTTTGAAGAGAATCTGCGCATTATCCGCTACTGCGTGAAGTACCGTAAGCGAATCATCTTCCCGTCAACTTCAGAAGTTTATGGGATGTGTAGCGATAAATACTTCGATGAGGACCATTCTAATTTAATCGTCGGCCCGGTGAATAAACCACGCTGGATTTATTCGGTATCAAAACAATTACTTGATCGGGTGATCTGGGCCTATGGCGAAAAAGAGGGTTTACAGTTCACCCTCTTCCGCCCGTTTAACTG  >  NC_000913/2369226‑2369516

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: