Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 624,133 T→C V9A (GTA→GCA)  entC → isochorismate hydroxymutase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913624,1330TC99.7% 1241.4 / ‑5.2 358V9A (GTA→GCA) entCisochorismate hydroxymutase 2
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (185/171);  total (187/171)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

Only 100 of 373 total aligned reads displayed.

GGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGTACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGTTTCGCCCGCT  >  NC_000913/624043‑624222
                                                                                          |                                                                                         
 gTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGcaca                                                                                        >  2:2701780/1‑93 (MQ=255)
 gTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGcaca                                                                                        <  2:1757467/93‑1 (MQ=255)
 gTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGcaca                                                                                        >  1:2497164/1‑93 (MQ=255)
  ttAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCAcag                                                                                       >  2:3899450/1‑93 (MQ=255)
  ttAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGAATACGTCACTGGCTGAGGAAGCAcag                                                                                       <  1:6754227/93‑1 (MQ=255)
   tAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCAcagc                                                                                      <  1:4436069/93‑1 (MQ=255)
    aGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCAcagca                                                                                     <  1:2058938/93‑1 (MQ=255)
    aGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCAcagca                                                                                     >  2:3741383/1‑93 (MQ=255)
    aGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCAcagca                                                                                     >  2:1578905/1‑93 (MQ=255)
     gCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCAcagcag                                                                                    >  1:8312606/1‑93 (MQ=255)
        aCCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGATGAAGCACAGCAGAcc                                                                                 >  2:2299332/1‑93 (MQ=255)
        aCCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGAAAAGCAGAcc                                                                                 >  1:3816243/1‑93 (MQ=255)
         ccGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCa                                                                                <  1:2788402/93‑1 (MQ=255)
          cGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCAt                                                                               <  1:3434933/93‑1 (MQ=255)
           gAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATg                                                                              <  2:4550829/93‑1 (MQ=255)
           gAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATg                                                                              >  2:834122/1‑93 (MQ=255)
           gAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATg                                                                              >  1:1373439/1‑93 (MQ=255)
            aaaaTATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATgg                                                                             >  2:8242128/1‑93 (MQ=255)
            aaaaTATAAATAATAATCATTATTAAAGCCTTTATCATTTGTGGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATgg                                                                             <  1:4690517/93‑1 (MQ=255)
             aaaTATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGc                                                                            <  1:5287198/93‑1 (MQ=255)
             aaaTATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGc                                                                            <  1:6542191/93‑1 (MQ=255)
                tataAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAac                                                                         >  1:5223157/1‑93 (MQ=255)
                tataAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAac                                                                         >  1:2692783/1‑93 (MQ=255)
                  taAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAacac                                                                       <  1:3269949/93‑1 (MQ=255)
                  taAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAacac                                                                       >  1:2110970/1‑93 (MQ=255)
                    aaTGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACtt                                                                     <  1:5909457/93‑1 (MQ=255)
                     aTGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTg                                                                    <  1:878238/93‑1 (MQ=255)
                      tGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCTCAGCAGACCATGGCAACACTTgc                                                                   >  2:4733259/1‑93 (MQ=255)
                      tGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTgc                                                                   >  1:6962358/1‑93 (MQ=255)
                      tGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTgc                                                                   >  2:1292025/1‑93 (MQ=255)
                         tAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGcc                                                                <  1:5576040/93‑1 (MQ=255)
                          aaTCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGccc                                                               >  2:7027886/1‑93 (MQ=255)
                          aaTCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGccc                                                               >  1:7455308/1‑93 (MQ=255)
                          aaTCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGcca                                                               <  1:7089771/93‑2 (MQ=255)
                            tCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCCCAGCAGACCATGGCAACACTTGCGCCCat                                                             >  1:4859517/1‑92 (MQ=255)
                            tCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCaa                                                             >  1:1707802/1‑93 (MQ=255)
                            tCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCaa                                                             <  2:3410236/93‑1 (MQ=255)
                             cATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAAt                                                            <  1:2657121/93‑1 (MQ=255)
                             cATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAAt                                                            <  2:3240444/93‑1 (MQ=255)
                              attattAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATc                                                           >  2:2920461/1‑93 (MQ=255)
                               ttattaAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCg                                                          >  2:182643/1‑93 (MQ=255)
                                tattaAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGc                                                         <  2:6635547/93‑1 (MQ=255)
                                 attaAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCt                                                        >  1:7868941/1‑93 (MQ=255)
                                  ttaAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCtt                                                       >  2:6168949/1‑93 (MQ=255)
                                  ttaAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCtt                                                       >  2:4498014/1‑93 (MQ=255)
                                   taAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCttt                                                      <  2:7727896/93‑1 (MQ=255)
                                     aaGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCttttt                                                    <  2:7357292/93‑1 (MQ=255)
                                      aGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTc                                                   >  2:278946/1‑93 (MQ=255)
                                      aGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTc                                                   <  1:8413799/93‑1 (MQ=255)
                                       gCCTTTATCATTTTGTGGAGTATGATATGGATACGTCACTGGCTGAGGAAGCGCAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCt                                                  <  2:4975940/93‑1 (MQ=255)
                                        ccTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCtt                                                 <  2:1124096/93‑1 (MQ=255)
                                         cTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCttt                                                >  2:5921261/1‑93 (MQ=255)
                                          tttATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTa                                               >  2:2228230/1‑93 (MQ=255)
                                          tttATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTa                                               >  2:6474257/1‑93 (MQ=255)
                                           ttATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTAt                                              <  2:8548976/93‑1 (MQ=255)
                                           ttATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTAt                                              >  1:503569/1‑93 (MQ=255)
                                            tATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATg                                             >  2:8423163/1‑93 (MQ=255)
                                             aTCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGt                                            >  2:5261509/1‑93 (MQ=255)
                                             aTCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGt                                            >  1:5480121/1‑93 (MQ=255)
                                             aTCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGt                                            <  1:427087/93‑1 (MQ=255)
                                              tCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTc                                           >  2:7550128/1‑93 (MQ=255)
                                                aTTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGc                                         <  2:5350693/93‑1 (MQ=255)
                                                aTTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGc                                         >  2:2812152/1‑93 (MQ=255)
                                                aTTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGc                                         >  1:4036807/1‑93 (MQ=255)
                                                aTTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGc                                         <  1:2232652/93‑1 (MQ=255)
                                                 ttttGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGcc                                        >  1:6487972/1‑93 (MQ=255)
                                                   ttGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGt                                      <  2:7333567/93‑1 (MQ=255)
                                                    tgtgGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTa                                     >  2:7985554/1‑93 (MQ=255)
                                                    tgtgGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTa                                     >  1:5137772/1‑93 (MQ=255)
                                                      tgGAGGATGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTAcc                                   <  2:3635438/93‑1 (MQ=255)
                                                       ggaggaTGATATGGATACGTCACTGGCTGAGGCAGCACAGCAGACCATGGCAACACTTGAGCCCAATCGCTTTTTCTTTATGTCGCCGTACCg                                  >  1:1142717/1‑93 (MQ=255)
                                                       ggaggaTGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCg                                  <  1:609250/93‑1 (MQ=255)
                                                       ggaggaTGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCg                                  <  1:2916471/93‑1 (MQ=255)
                                                         aggaTGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCa                                <  2:1971169/93‑1 (MQ=255)
                                                         aggaTGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCa                                >  1:762843/1‑93 (MQ=255)
                                                         aggaTGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCa                                <  2:5749089/93‑1 (MQ=255)
                                                          ggaTGATATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAg                               >  2:8766447/1‑93 (MQ=255)
                                                              gatATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGtttt                           >  1:1760843/1‑93 (MQ=255)
                                                              gatATGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGtttt                           <  2:6832861/93‑1 (MQ=255)
                                                                tatGGATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTac                         <  1:1538764/93‑1 (MQ=255)
                                                                   ggATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTacgac                      <  2:5968772/93‑1 (MQ=255)
                                                                   ggATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTacgac                      <  2:4272327/93‑1 (MQ=255)
                                                                    gATACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTacgacg                     >  1:3699764/1‑93 (MQ=255)
                                                                      tACGTCACTGGCTGAGGTAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCTCAGATTTACGACGTc                   >  2:7843974/1‑93 (MQ=255)
                                                                      tACGTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTc                   <  2:8628600/93‑1 (MQ=255)
                                                                         gTCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAgg                >  1:2332925/1‑93 (MQ=255)
                                                                          tCACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGa               <  2:621500/93‑1 (MQ=255)
                                                                           cACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTCCGACGTCAGGAt              >  2:3102169/1‑93 (MQ=255)
                                                                           cACTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGCTTTACGACGTCAGGAt              >  2:5550952/1‑93 (MQ=255)
                                                                            aCTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATg             >  1:8669116/1‑93 (MQ=255)
                                                                            aCTGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATg             >  2:246829/1‑93 (MQ=255)
                                                                             cTGGCTGAGGAAGCACAGCAGGCCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGt            <  2:6368494/93‑1 (MQ=255)
                                                                              tGGCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGtt           >  1:7745676/1‑93 (MQ=255)
                                                                                gCTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGTTTc         <  1:6341918/93‑1 (MQ=255)
                                                                                 cTGAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGTTTCg        <  1:5409171/93‑1 (MQ=255)
                                                                                   gAGGAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGTTTCGcc      >  1:902346/1‑93 (MQ=255)
                                                                                    aGGAAGCACAGCAGACCCTGGCAGCGCTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGTTTCGccc     <  2:4841570/93‑1 (MQ=255)
                                                                                      gAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGTTTCGCCCGc   >  1:5522620/1‑93 (MQ=255)
                                                                                      gAAGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGTTTCGCCCGc   <  1:1200152/93‑1 (MQ=255)
                                                                                       aaGCACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGTTTCGCCCGCt  <  2:8061825/93‑1 (MQ=255)
                                                                                          |                                                                                         
GGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGTACAGCAGACCATGGCAACACTTGCGCCCAATCGCTTTTTCTTTATGTCGCCGTACCGCAGTTTTACGACGTCAGGATGTTTCGCCCGCT  >  NC_000913/624043‑624222

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: