Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,858,723 C→T A212T (GCG→ACG)  b1776 ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,858,7230CT85.3% 603.4 / 67.6 226A212T (GCG→ACG) b1776putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (17/16);  new base T (80/112);  total (98/128)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.52e-01

Only 100 of 226 total aligned reads displayed.

AAATTTGCTGGCAGCGCGCCACCACATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGCCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGTGGGCCAGTGCCG  >  NC_000913/1858633‑1858814
                                                                                          |                                                                                           
aaaTTTGCTGGCAGCGCGCCACCACATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTca                                                                                           >  1:6542/1‑93 (MQ=255)
  aTTTGCTGGCAGCGCGCCACCACATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTcaca                                                                                         <  1:4247820/93‑1 (MQ=255)
   tttGCTGGCAGCGCGCCACCACATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACAt                                                                                        >  2:6096270/1‑93 (MQ=255)
      gCTGGCAGCGCGCCACCACATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACAtctt                                                                                     <  2:7036070/93‑1 (MQ=255)
      gCTGGCAGCGCGCCACCACATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACAtctt                                                                                     <  1:2721143/93‑1 (MQ=255)
       cTGGCAGCGCGCCACCACATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACAtcttc                                                                                    <  2:5847120/93‑1 (MQ=255)
       cTGGCAGCGCGCCACCACATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACAtcttc                                                                                    <  1:4949746/93‑1 (MQ=255)
         ggCAGCGCGCCACCACATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATcttctt                                                                                  <  2:6862679/93‑1 (MQ=255)
         ggCAGCGCGCCACCACATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATcttctt                                                                                  <  2:5012491/93‑1 (MQ=255)
           cAGCGCGCCACCACATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGc                                                                                <  1:5799897/93‑1 (MQ=255)
            aGCGCGCCACCACATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCa                                                                               <  2:1566090/93‑1 (MQ=255)
            aGCGCGCCACCACATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCa                                                                               >  1:5489237/1‑93 (MQ=255)
             gcgcgcCACCACATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAg                                                                              >  1:5405939/1‑93 (MQ=255)
             gcgcgcCACCACATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAg                                                                              >  2:4000869/1‑93 (MQ=255)
              cgcgcCACCACATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGa                                                                             >  1:5888638/1‑93 (MQ=255)
              cgcgcCACCACATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGa                                                                             <  2:6567740/93‑1 (MQ=255)
                 gcCACCACATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCa                                                                          >  1:3114860/1‑93 (MQ=255)
                 gcCACCACATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGCCACATCTTCTTGCAGACCa                                                                          >  1:5614319/1‑93 (MQ=255)
                  ccaccaCATCTTCGGTAGAACCATTTACTACTGCCGTCGCCCCCAGTTCTTTTGCAACCGGGAAGCGGACCGCCACATCTTCTTGCAGACCaa                                                                         >  1:1315963/1‑93 (MQ=255)
                  ccaccaCATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCaa                                                                         <  1:4957800/93‑1 (MQ=255)
                  ccaccaCATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCaa                                                                         <  2:513849/93‑1 (MQ=255)
                   caccacATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAac                                                                        <  2:1816338/93‑1 (MQ=255)
                   caccacATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAac                                                                        >  1:7046700/1‑93 (MQ=255)
                   caccacATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAac                                                                        <  1:3058689/93‑1 (MQ=255)
                   caccacATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAac                                                                        >  2:6368654/1‑93 (MQ=255)
                    accacATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAacg                                                                       >  2:5850305/1‑93 (MQ=255)
                     ccacATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAacga                                                                      >  2:3692870/1‑93 (MQ=255)
                     ccacATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAacga                                                                      >  2:6504222/1‑93 (MQ=255)
                     ccacATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAacga                                                                      >  2:1361453/1‑93 (MQ=255)
                     ccacATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAacga                                                                      >  1:1586384/1‑93 (MQ=255)
                     ccacATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAacga                                                                      <  2:4638570/93‑1 (MQ=255)
                      cacaTCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAacgac                                                                     >  1:3588369/1‑93 (MQ=255)
                      cacaTCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAacgac                                                                     <  2:3866975/93‑1 (MQ=255)
                        caTCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAacgacga                                                                   >  2:1223382/1‑93 (MQ=255)
                        caTCTTCGGTAGAACCATGTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAacgacga                                                                   <  2:565386/93‑1 (MQ=255)
                          tCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAacgacgacg                                                                 >  1:3810305/1‑93 (MQ=255)
                          tCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAacgacgacg                                                                 >  2:737859/1‑93 (MQ=255)
                           cTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAAcgacgacga                                                                <  2:2269491/93‑1 (MQ=255)
                           cTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAAcgacgacga                                                                <  2:2171062/93‑1 (MQ=255)
                           cTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAAcgacgacga                                                                <  1:1723404/93‑1 (MQ=255)
                            ttCGGTAGACCCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGCCACATCTTCTTGCAGACCAACGACGACGat                                                               <  1:6623152/93‑1 (MQ=255)
                            ttCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGat                                                               <  1:4842974/93‑1 (MQ=255)
                               ggTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATAtt                                                            <  1:4147361/93‑1 (MQ=255)
                               ggTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATAtt                                                            <  1:1191896/93‑1 (MQ=255)
                                 tAGAACCATTTACTACTGCCGTCGCCCCCAGTTCTTTTGCAACCGGGAAGCGGACCGCCACATCTTCTTGCAGACCAACGACGACGATATTTa                                                          <  1:6260221/93‑1 (MQ=255)
                                  aGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTAc                                                         >  1:5318238/1‑93 (MQ=255)
                                  aGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTAc                                                         <  1:1655831/93‑1 (MQ=255)
                                   gAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGCCACATCTTCTTGCAGACCAACGACGACGATATTTAcc                                                        <  2:2839295/93‑1 (MQ=255)
                                   gAACAATTTGCTACTGCCGTCGCCCCTATTTCTTTTGCAACCGGGAAGCGGACCGCCACATCTTCTTGCAGACCAACGACGACGATATTTAcc                                                        <  1:2915542/93‑1 (MQ=255)
                                     aCCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCATCACATCTTCTTGCAGACCAACGACGACGATATTTACCGc                                                      >  2:5001888/1‑93 (MQ=255)
                                       cATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGccc                                                    <  2:4459096/93‑1 (MQ=255)
                                        aTTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGCCACATCTTCTTGCAGACCAACGACGACGATATATACCGccca                                                   >  2:5778941/1‑92 (MQ=255)
                                         tttACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGccccc                                                  >  2:3234664/1‑93 (MQ=255)
                                          ttACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCa                                                 <  2:7217552/93‑1 (MQ=255)
                                           tactacTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCAt                                                <  2:2499423/93‑1 (MQ=255)
                                             ctactGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATaa                                              <  2:765392/93‑1 (MQ=255)
                                             ctactGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATaa                                              <  1:1431931/93‑1 (MQ=255)
                                             ctactGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATaa                                              <  1:3843999/93‑1 (MQ=255)
                                              tactGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAAt                                             >  1:1072464/1‑93 (MQ=255)
                                              tactGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAAt                                             >  1:777060/1‑93 (MQ=255)
                                              tactGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAAt                                             >  1:2479493/1‑93 (MQ=255)
                                               actGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAAtt                                            >  2:3174060/1‑93 (MQ=255)
                                                ctGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTATCGCCCCCATAATTc                                           <  2:5268022/93‑1 (MQ=255)
                                                ctGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTc                                           <  2:5159910/93‑1 (MQ=255)
                                                ctGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTc                                           >  1:1942536/1‑93 (MQ=255)
                                                 tGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTcg                                          >  2:5505793/1‑93 (MQ=255)
                                                 tGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTcg                                          <  2:3398492/93‑1 (MQ=255)
                                                      tCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCa                                     >  2:6476660/1‑93 (MQ=255)
                                                       cGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCAt                                    <  1:1987802/93‑1 (MQ=255)
                                                        gCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGCCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCAtt                                   <  1:2526548/93‑1 (MQ=255)
                                                        gCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGCCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCAtt                                   >  1:3417803/1‑93 (MQ=255)
                                                          aCCCAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTg                                 <  1:1536091/93‑1 (MQ=255)
                                                           cccAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGt                                >  1:5239099/1‑93 (MQ=255)
                                                           cccAGTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACAACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGt                                >  1:3519704/1‑93 (MQ=255)
                                                            ccAGTTCTTTTGCAACCGGGAAGCGGACCGCCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTa                               <  1:6509355/93‑1 (MQ=255)
                                                             cAGTTCTTTTGCAACCGGGAAGCGGACCGCCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTAc                              >  1:633316/1‑93 (MQ=255)
                                                               gTTCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACgg                            >  2:6169629/1‑93 (MQ=255)
                                                                ttCTTTTGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGga                           >  1:1203249/1‑93 (MQ=255)
                                                                    tttGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCATCGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAAc                       <  2:7257373/93‑1 (MQ=255)
                                                                    tttGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAAc                       <  2:469693/93‑1 (MQ=255)
                                                                    tttGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAAc                       <  2:4793859/93‑1 (MQ=255)
                                                                    tttGCAACCGGGAAGCGGACCGCCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAAc                       <  1:7402676/93‑1 (MQ=255)
                                                                     ttGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGcaca                      <  1:6974557/93‑1 (MQ=255)
                                                                     ttGCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACa                      <  1:927700/93‑1 (MQ=255)
                                                                       gCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGc                    <  1:4390195/93‑1 (MQ=255)
                                                                       gCAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGc                    >  1:3915949/1‑93 (MQ=255)
                                                                        cAACCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGcc                   <  1:7015165/93‑1 (MQ=255)
                                                                        cAACCGGGAAGCGGACCGCCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGcc                   <  2:5938667/93‑1 (MQ=255)
                                                                         aaCCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGccc                  <  2:1845705/93‑1 (MQ=255)
                                                                          aCCGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCg                 >  1:4997893/1‑93 (MQ=255)
                                                                            cGGGAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAg               >  2:4210522/1‑93 (MQ=255)
                                                                             gggAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAgc              <  1:7299128/93‑2 (MQ=255)
                                                                              ggAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAgcg             >  2:6256211/1‑91 (MQ=255)
                                                                              ggAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAgcg             <  1:211026/93‑3 (MQ=255)
                                                                               gAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCgg            >  2:4887900/1‑93 (MQ=255)
                                                                                 aGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGc          <  2:2899051/93‑1 (MQ=255)
                                                                                    ggACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAg       <  2:2737879/93‑1 (MQ=255)
                                                                                    ggACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAg       <  2:7329815/93‑1 (MQ=255)
                                                                                     gACCGCCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGTGGGCCAGt      <  2:4347365/93‑1 (MQ=255)
                                                                                        cGCCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGTGGGCCAGTGcc   <  2:1179133/93‑1 (MQ=255)
                                                                                          |                                                                                           
AAATTTGCTGGCAGCGCGCCACCACATCTTCGGTAGAACCATTTACTACTGCCGTCGCACCCAGTTCTTTTGCAACCGGGAAGCGGACCGCCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGTGGGCCAGTGCCG  >  NC_000913/1858633‑1858814

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