Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,869,012 T→C I116V (ATT→GTT)  b3691 ← D‑galactonate transport

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,869,0120TC100.0% 1232.6 / NA 387I116V (ATT→GTT) b3691D‑galactonate transport
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (154/233);  total (154/233)

Only 100 of 402 total aligned reads displayed.

TATAAAAACCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAATACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAA  >  NC_000913/3868922‑3869101
                                                                                          |                                                                                         
tataAAAACCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAacac                                                                                         >  1:807074/1‑93 (MQ=255)
tataAAAACCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAacac                                                                                         <  2:4856595/93‑1 (MQ=255)
 ataAAAACCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACAcc                                                                                        >  1:6102641/1‑93 (MQ=255)
 ataAAAACCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACAcc                                                                                        <  1:6301070/93‑1 (MQ=255)
 ataAAAACCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACAcc                                                                                        <  1:3900566/93‑1 (MQ=255)
  taAAAACCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCg                                                                                       >  1:5843246/1‑93 (MQ=255)
  taAAAACCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCg                                                                                       <  1:1488084/93‑1 (MQ=255)
  taAAAACCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCg                                                                                       <  1:6412636/93‑1 (MQ=255)
  taAAAACCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCg                                                                                       >  2:7780515/1‑93 (MQ=255)
      aaCCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTa                                                                                   <  1:521307/93‑1 (MQ=255)
      aaCCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTa                                                                                   <  2:5300801/93‑1 (MQ=255)
      aaCCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTa                                                                                   <  2:7334224/93‑1 (MQ=255)
      aaCCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTa                                                                                   <  2:7423410/93‑1 (MQ=255)
      aaCCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTa                                                                                   >  2:4607467/1‑93 (MQ=255)
      aaCCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTa                                                                                   >  2:7839245/1‑93 (MQ=255)
      aaCCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTa                                                                                   >  2:1768489/1‑93 (MQ=255)
        ccAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATc                                                                                 <  1:3348015/93‑1 (MQ=255)
        ccAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATc                                                                                 >  2:7084260/1‑93 (MQ=255)
        ccAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATc                                                                                 <  2:2490818/93‑1 (MQ=255)
         cAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATcg                                                                                >  2:6991093/1‑93 (MQ=255)
          aaCGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATcgc                                                                               >  1:3835070/1‑93 (MQ=255)
           aCGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATcgcg                                                                              <  1:739719/93‑1 (MQ=255)
               cAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCa                                                                          <  2:5622946/93‑1 (MQ=255)
               cAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCa                                                                          >  1:2611862/1‑93 (MQ=255)
               cAGAAGCGCGTTCATGTGCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCa                                                                          <  2:3423639/93‑1 (MQ=255)
                 gAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCGGGCGCTTCGAAAACACCGGTTATCGCGCGCAGa                                                                        >  2:4074798/1‑93 (MQ=255)
                 gAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGa                                                                        <  2:6711573/93‑1 (MQ=255)
                 gAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGa                                                                        >  1:6810319/1‑93 (MQ=255)
                   aGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGAcc                                                                      <  1:7716354/93‑1 (MQ=255)
                        gTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTa                                                                 >  1:5120662/1‑93 (MQ=255)
                        gTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTa                                                                 <  2:4395759/93‑1 (MQ=255)
                         ttCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTaa                                                                <  2:1045550/93‑1 (MQ=255)
                         ttCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTaa                                                                >  1:2292958/1‑93 (MQ=255)
                         ttCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTaa                                                                <  2:2913125/93‑1 (MQ=255)
                         ttCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTaa                                                                <  2:6329419/93‑1 (MQ=255)
                         ttCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTaa                                                                >  1:3977263/1‑93 (MQ=255)
                          tCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAAt                                                               <  2:135137/93‑1 (MQ=255)
                          tCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAAt                                                               <  1:3285065/93‑1 (MQ=255)
                            aTGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGa                                                             <  2:2235819/93‑1 (MQ=255)
                            aTGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGa                                                             <  1:959271/93‑1 (MQ=255)
                             tGTTCCGGGAACCAGCTGGTAACCATCCGGTTATTGGTCGGGAAAGCAGGAGCTTCGAAAAAACAGGTTATCGCGCGCAGCACAATTAATGaa                                                            >  2:5061592/1‑92 (MQ=255)
                              gTTCCGGGAACCAGCTGGTCGCCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACa                                                           <  1:3765855/93‑1 (MQ=255)
                               ttCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACAt                                                          <  2:1490554/93‑1 (MQ=255)
                               ttCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACAt                                                          >  2:1966925/1‑93 (MQ=255)
                                 ccGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTa                                                        >  2:7726917/1‑93 (MQ=255)
                                  cGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTaa                                                       <  1:2347884/93‑1 (MQ=255)
                                       aCCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGccct                                                  >  1:5641986/1‑92 (MQ=255)
                                       aCCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCg                                                  >  1:2764355/1‑93 (MQ=255)
                                         cAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTg                                                <  1:5562108/93‑1 (MQ=255)
                                         cAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTg                                                >  2:597496/1‑93 (MQ=255)
                                           gCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGc                                              <  1:4381648/93‑1 (MQ=255)
                                           gCTGGGCCCCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCCGACCAATTAATGACATTAAGCCCGTGGc                                              <  1:2204354/93‑1 (MQ=255)
                                              ggTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCaaa                                           <  1:5937217/93‑1 (MQ=255)
                                                tCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGc                                         >  2:5673218/1‑93 (MQ=255)
                                                tCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGc                                         <  1:3131297/93‑1 (MQ=255)
                                                 cACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGATTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGcc                                        <  1:125956/93‑1 (MQ=255)
                                                  aCCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGccc                                       <  1:4482150/93‑1 (MQ=255)
                                                   ccATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGAGCAATTAATGACATTAAGCCCGTGGCAAAGCCCt                                      >  1:3033172/1‑93 (MQ=255)
                                                   ccATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCt                                      <  2:6241316/93‑1 (MQ=255)
                                                    cATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTg                                     <  1:1931976/93‑1 (MQ=255)
                                                     aTCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTgg                                    <  2:794090/93‑1 (MQ=255)
                                                      tCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGa                                   >  2:6882871/1‑93 (MQ=255)
                                                       ccGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGaa                                  <  1:2509700/93‑1 (MQ=255)
                                                         ggTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATa                                >  2:537172/1‑93 (MQ=255)
                                                          gTTATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATaa                               <  1:6750661/93‑1 (MQ=255)
                                                           ttATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATaaa                              >  2:2181297/1‑93 (MQ=255)
                                                           ttATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATaaa                              >  1:6939479/1‑93 (MQ=255)
                                                           ttATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATaaa                              >  1:7532187/1‑93 (MQ=255)
                                                           gtATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATaaa                              <  1:7259469/92‑1 (MQ=255)
                                                            tATTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAg                             <  2:2403973/93‑1 (MQ=255)
                                                             aTTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGt                            <  1:1202783/93‑1 (MQ=255)
                                                             aTTGGTCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGt                            <  2:3710191/93‑1 (MQ=255)
                                                                 ggcGGGAAACCGGCCACTTCAAAACCCCCGTTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGcc                        <  2:1575882/91‑1 (MQ=255)
                                                                  tCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCa                       <  1:3421616/93‑1 (MQ=255)
                                                                  tCGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCa                       <  1:5274148/93‑1 (MQ=255)
                                                                   cGGGAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCAc                      <  2:4244358/93‑1 (MQ=255)
                                                                    gggAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACt                     <  1:1025394/93‑1 (MQ=255)
                                                                     ggAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTg                    >  2:2082194/1‑93 (MQ=255)
                                                                     ggAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTg                    >  1:49782/1‑93 (MQ=255)
                                                                      gAAAGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGa                   <  1:4589764/93‑1 (MQ=255)
                                                                        aaGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGAcc                 <  2:5771614/93‑1 (MQ=255)
                                                                        aaGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGAcc                 >  1:7335900/1‑93 (MQ=255)
                                                                         aGCAGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCa                >  2:3956020/1‑93 (MQ=255)
                                                                            aGGCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGcc             >  2:1150496/1‑93 (MQ=255)
                                                                             ggCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCa            >  1:841339/1‑93 (MQ=255)
                                                                             ggCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCa            >  2:2750718/1‑93 (MQ=255)
                                                                             ggCGCTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCa            >  1:3235571/1‑93 (MQ=255)
                                                                              gcgcTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCaa           >  2:301346/1‑93 (MQ=255)
                                                                               cgcTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCaaa          <  1:6629420/93‑1 (MQ=255)
                                                                               cgcTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCaaa          <  1:4892956/93‑1 (MQ=255)
                                                                                gcTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAat         <  1:4257674/93‑1 (MQ=255)
                                                                                gcTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAat         <  2:3075162/93‑1 (MQ=255)
                                                                                 cTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAata        >  1:2928825/1‑93 (MQ=255)
                                                                                 cTTCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAata        <  2:2292957/93‑1 (MQ=255)
                                                                                  ttCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAatat       >  2:6454402/1‑93 (MQ=255)
                                                                                  ttCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAatat       <  1:4110959/93‑1 (MQ=255)
                                                                                  ttCGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAatat       >  1:5402422/1‑93 (MQ=255)
                                                                                    cGAAAACACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCg     <  2:3247390/93‑1 (MQ=255)
                                                                                       aaaCACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATCAAGTCGCCACTGACCAGCCAAATATCGCaa  >  1:3581570/1‑93 (MQ=255)
                                                                                       aaaCACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCaa  >  2:3558461/1‑93 (MQ=255)
                                                                                          |                                                                                         
TATAAAAACCAACGGCAGAAGCGCGTTCATGTTCCGGGAACCAGCTGGTCACCATCCGGTTATTGGTCGGGAAAGCAGGCGCTTCGAAAATACCGGTTATCGCGCGCAGACCAATTAATGACATTAAGCCCGTGGCAAAGCCCTGGAATAAAGTCGCCACTGACCAGCCAAATATCGCAA  >  NC_000913/3868922‑3869101

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