Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,119,304 T→C N169S (AAC→AGC)  solA ← sarcosine oxidase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,119,3040TC85.0% 571.8 / 88.9 246N169S (AAC→AGC) solAsarcosine oxidase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (18/19);  new base C (102/107);  total (120/126)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.62e-02

Only 100 of 248 total aligned reads displayed.

CGCAGACAATCGCTTTTTTCGCCTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGTTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCTGCGCAAAAAACCGGAATCAGTTTCAA  >  NC_000913/1119212‑1119394
                                                                                            |                                                                                          
cgcAGACAATCGCTTTTTTCGCCTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGt                                                                                            <  2:1849190/93‑1 (MQ=255)
 gcAGACAATCGCTTTTTTCGCCTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGtt                                                                                           >  2:1326808/1‑93 (MQ=255)
   aGACAATCGCTTTTTTCGCCTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGa                                                                                         <  1:6058248/93‑1 (MQ=255)
    gACAATCGCTTTTTTCGCCTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGaa                                                                                        >  1:5995572/1‑93 (MQ=255)
    gACAATCGCTTTTTTCGCCTGATACTCACCGACAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGaa                                                                                        >  2:1820646/1‑93 (MQ=255)
     aCAATCGCTTTTTTCGCCTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGTTGAAc                                                                                       <  1:3547506/93‑1 (MQ=255)
     aCAATCGCTTTTTTCGCCTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAAc                                                                                       >  2:5152215/1‑93 (MQ=255)
     aCAATCGCTTTTTTCGCCTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAAc                                                                                       >  2:6909524/1‑93 (MQ=255)
     aCAATCGCTTTTTTCGCCTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAAc                                                                                       <  2:4878064/93‑1 (MQ=255)
     aCAATCGCTTTTTTCGCCTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTAACCGGGCAGCTGAAc                                                                                       <  1:2851539/93‑1 (MQ=255)
        aTCGCTTTTTTCGCCTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGt                                                                                    <  2:6556993/93‑1 (MQ=255)
         tCGCTTTTTTCGCCTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGtt                                                                                   >  1:4805196/1‑93 (MQ=255)
         tCGCTTTTTTCGCCTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGtt                                                                                   >  2:3651658/1‑93 (MQ=255)
           gCTTTTTTCGCCTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTgc                                                                                 >  2:7377174/1‑93 (MQ=255)
            cTTTTTTCGCCTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTgcg                                                                                >  2:349381/1‑93 (MQ=255)
             ttttttCGCCTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTgcgc                                                                               >  2:3879121/1‑93 (MQ=255)
                tttCGCCTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGcaca                                                                            <  2:1701856/93‑1 (MQ=255)
                 ttCGCCTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAg                                                                           <  2:5814494/93‑1 (MQ=255)
                 ttCGCCTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAg                                                                           <  2:7075129/93‑1 (MQ=255)
                  tCGCCTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGc                                                                          <  2:214178/93‑1 (MQ=255)
                    gCCTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGTTGAACAGTTGCGCACAGccc                                                                        >  2:3286397/1‑93 (MQ=255)
                    gCCTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGTTGAACAGTTGCGCACAGccc                                                                        >  2:2066076/1‑93 (MQ=255)
                    gCCTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGccc                                                                        >  1:2238234/1‑93 (MQ=255)
                     ccTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCg                                                                       >  1:362560/1‑93 (MQ=255)
                     ccTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCg                                                                       >  1:6674814/1‑93 (MQ=255)
                      cTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGTTGAACAGTTGCGCACAGCCCGc                                                                      <  1:7219827/93‑1 (MQ=255)
                      cTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGc                                                                      >  2:2795302/1‑93 (MQ=255)
                      cTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGc                                                                      <  2:5078717/93‑1 (MQ=255)
                        gATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGctt                                                                    >  2:6380111/1‑93 (MQ=255)
                          tACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGcttcc                                                                  >  2:4155390/1‑93 (MQ=255)
                          tACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGcttcc                                                                  <  2:1454468/93‑1 (MQ=255)
                          tACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGcttcc                                                                  <  2:5607017/93‑1 (MQ=255)
                          tACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGcttcc                                                                  >  1:3795354/1‑93 (MQ=255)
                             tCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGcttccttc                                                               <  1:6548703/93‑1 (MQ=255)
                             tCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGcttccttc                                                               <  1:1370259/93‑1 (MQ=255)
                              cACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCg                                                              >  1:3912251/1‑93 (MQ=255)
                                ccGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGcc                                                            <  2:5910988/93‑1 (MQ=255)
                                ccGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGcc                                                            >  1:7177038/1‑93 (MQ=255)
                                 cGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCa                                                           <  1:2562491/93‑1 (MQ=255)
                                  gTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAg                                                          >  1:5777781/1‑93 (MQ=255)
                                   tCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGt                                                         >  1:96323/1‑93 (MQ=255)
                                    cAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGtt                                                        <  1:4874503/93‑1 (MQ=255)
                                     aGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTg                                                       >  2:4335055/1‑93 (MQ=255)
                                     aGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTg                                                       <  1:3127498/93‑1 (MQ=255)
                                      gCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTgg                                                      >  1:7466263/1‑93 (MQ=255)
                                      gCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTgg                                                      >  2:6658176/1‑93 (MQ=255)
                                       ccGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGa                                                     <  1:4681890/93‑1 (MQ=255)
                                        cGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGAt                                                    <  2:6192730/93‑1 (MQ=255)
                                        cGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGAt                                                    >  1:1987421/1‑93 (MQ=255)
                                         gTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATc                                                   <  2:1800121/93‑1 (MQ=255)
                                         gTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATc                                                   <  2:2529027/93‑1 (MQ=255)
                                         gTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATc                                                   <  2:2737479/93‑1 (MQ=255)
                                         gTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATc                                                   >  1:5863552/1‑93 (MQ=255)
                                          tttCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATcc                                                  <  2:4273887/93‑1 (MQ=255)
                                           ttCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCa                                                 <  1:6119735/93‑1 (MQ=255)
                                            tCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAg                                                <  2:7112434/93‑1 (MQ=255)
                                             cAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGTTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAgg                                               >  1:7310281/1‑93 (MQ=255)
                                               aTAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGtt                                             <  1:1047686/93‑1 (MQ=255)
                                                  gTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTa                                          <  1:1622827/93‑1 (MQ=255)
                                                   ttACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTaa                                         <  2:382164/93‑1 (MQ=255)
                                                    tACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGTTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAAt                                        <  1:1124517/93‑1 (MQ=255)
                                                     aCGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATc                                       <  2:40016/93‑1 (MQ=255)
                                                     aCGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATc                                       >  2:4450358/1‑93 (MQ=255)
                                                         catcatCGTCATGACGAATTGCGGTGACCGGGCAGTTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCa                                   <  2:5422280/93‑1 (MQ=255)
                                                         catcatCGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCa                                   <  1:4529479/93‑1 (MQ=255)
                                                           tcatcGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGt                                 >  1:3070072/1‑93 (MQ=255)
                                                           tcatcGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGt                                 <  1:7045278/93‑1 (MQ=255)
                                                           tcatcGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGt                                 >  1:458387/1‑93 (MQ=255)
                                                             atcGTCATGACGAATTGCGGTGACCGGGCAGTTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTc                               >  2:699441/1‑93 (MQ=255)
                                                             atcGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTc                               >  1:1275162/1‑93 (MQ=255)
                                                             atcGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTc                               >  1:4008188/1‑93 (MQ=255)
                                                             atcGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTc                               >  1:1561585/1‑93 (MQ=255)
                                                              tcGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCg                              >  2:6746818/1‑93 (MQ=255)
                                                              tcGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCg                              >  2:344387/1‑93 (MQ=255)
                                                               cGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGc                             <  2:1635082/93‑1 (MQ=255)
                                                               cGTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGc                             <  2:1139480/93‑1 (MQ=255)
                                                                gTCATGACGAATTGCGGTGACCGGGCAGTTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCt                            <  1:5398636/93‑1 (MQ=255)
                                                                gTCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCt                            <  1:6344237/93‑1 (MQ=255)
                                                                 tCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCTg                           <  2:3929450/93‑1 (MQ=255)
                                                                 tCATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCTg                           <  1:1404623/93‑1 (MQ=255)
                                                                  cATGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCTgc                          <  1:6844471/93‑1 (MQ=255)
                                                                    tGACGAATTGCGGTGACCGGGCAGTTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCTgcgc                        >  2:1221037/1‑93 (MQ=255)
                                                                    tGACGAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCTgcgc                        >  1:5824607/1‑93 (MQ=255)
                                                                     gACGAATTGCGGTGACCGGGCAGTTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCTGCGCa                       <  2:4955893/93‑1 (MQ=255)
                                                                        gAATTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCTGCGCaaaa                    <  2:4784364/93‑1 (MQ=255)
                                                                          aTTGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCTGCGCaaaaaa                  <  2:1050249/93‑1 (MQ=255)
                                                                            tGCGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCTGCGCAAAAAAcc                <  2:2330405/93‑1 (MQ=255)
                                                                             gCGGTGACCGGGCAGTTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCTGCGCAAAAAACCg               >  1:6381795/1‑93 (MQ=255)
                                                                              cGGTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCTGCGCAAAAAACCgg              >  2:4550408/1‑93 (MQ=255)
                                                                               ggTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCTGCGCAAAAAACCGGa             >  1:548912/1‑93 (MQ=255)
                                                                                gTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCTGCGCAAAAAACCGGaa            <  2:817758/93‑1 (MQ=255)
                                                                                gTGACCGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCTGCGCAAAAAACCGGaa            <  2:2488546/93‑1 (MQ=255)
                                                                                    ccGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCTGCGCAAAAAACCGGAATCAg        <  2:3210241/93‑1 (MQ=255)
                                                                                     cGGGCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCTGCGCAAAAAACCGGAATCAGt       <  1:5644308/93‑1 (MQ=255)
                                                                                      gggCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCTGCGCAAAAAACCGGAATCAGtt      <  1:3259160/93‑1 (MQ=255)
                                                                                      gggCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCTGCGCAAAAAACCGGAATCAGtt      <  1:1810377/93‑1 (MQ=255)
                                                                                      gggCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCTGCGCAAAAAACCGGAATCAGtt      <  1:4481746/93‑1 (MQ=255)
                                                                                        gCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCTGCGCAAAAAACCGGAATCAGTTTc    >  1:4330947/1‑93 (MQ=255)
                                                                                        gCAGCTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCTGCGCAAAAAACCGGAATCAGTTTc    >  2:7200330/1‑93 (MQ=255)
                                                                                          aGTTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCTGCGCAAAAAACCGGAATCAGTTTCaa  <  2:2919379/93‑1 (MQ=255)
                                                                                            |                                                                                          
CGCAGACAATCGCTTTTTTCGCCTGATACTCACCGTCAGCCGTTTCAATAGTTACGCCATCATCGTCATGACGAATTGCGGTGACCGGGCAGTTGAACAGTTGCGCACAGCCCGCTTCCTTCGCCAGTTGGATCCAGGTTTTAATCGCCAGTTCGCTGCGCAAAAAACCGGAATCAGTTTCAA  >  NC_000913/1119212‑1119394

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