Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 675,023 A→G Y250H (TAT→CAT)  ybeQ ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913675,0230AG85.1% 626.9 / 95.2 242Y250H (TAT→CAT) ybeQhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/15);  new base G (114/92);  total (135/107)
Fisher's exact test for biased strand distribution p-value = 8.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 246 total aligned reads displayed.

ACCAGGCAGTGGCGTCATTGCATTCAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACC  >  NC_000913/674933‑675114
                                                                                          |                                                                                           
aCCAGGCAGTGGCGTCATTGCATTCAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCa                                                                                           >  2:153012/1‑93 (MQ=255)
aCCAGGCAGTGGCGTCATTGCATTCAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCa                                                                                           <  1:3325954/93‑1 (MQ=255)
  cAGGCAGTGGCGTCATTGCATTCAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCata                                                                                         >  1:4973248/1‑93 (MQ=255)
  cAGGCAGTGGCGTCATTGCATTCAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCata                                                                                         >  2:6643037/1‑93 (MQ=255)
  cAGGCAGTGGCGTCATTGCATTCAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCata                                                                                         <  2:5789062/93‑1 (MQ=255)
  cAGGCAGTGGCGTCATTGCATTCAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCata                                                                                         >  2:2037600/1‑93 (MQ=255)
   aGGCAGTGGCGTCATTGCATTCAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCatat                                                                                        <  1:2307018/93‑1 (MQ=255)
    ggCAGTGGCGTCATTGCATTCAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCAtata                                                                                       >  2:887604/1‑93 (MQ=255)
    ggCAGTGGCGTCATTGCATTCAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACAtata                                                                                       <  1:6404709/93‑1 (MQ=255)
     gCAGTGGCGTCATTGCATTCAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATaa                                                                                      >  2:4542715/1‑93 (MQ=255)
          ggCGTCATTGCATTCAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGAt                                                                                 <  1:4834213/93‑1 (MQ=255)
           gCGTCATTGCATTCAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATg                                                                                <  2:2973297/93‑1 (MQ=255)
            cGTCATTGCATTCAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGt                                                                               >  2:4600004/1‑93 (MQ=255)
             gTCATTGCATTCAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGtt                                                                              >  2:3323895/1‑93 (MQ=255)
              tCATTGCATTCAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTg                                                                             >  2:3470250/1‑93 (MQ=255)
              tCATTGCATTCAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTg                                                                             >  2:3969035/1‑93 (MQ=255)
               cATTGCATTCAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGa                                                                            >  1:785823/1‑93 (MQ=255)
                aTTGCATTCAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGAc                                                                           >  2:1607987/1‑93 (MQ=255)
                 ttGCATTCAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACg                                                                          <  2:3274283/93‑1 (MQ=255)
                 ttGCATTCAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACg                                                                          <  2:570986/93‑1 (MQ=255)
                    cATTCAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTaa                                                                       <  1:418540/93‑1 (MQ=255)
                     aTTCAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAg                                                                      >  1:6332817/1‑93 (MQ=255)
                        cAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAgcgc                                                                   >  1:2707067/1‑93 (MQ=255)
                        cAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAgcgc                                                                   <  2:84270/93‑1 (MQ=255)
                         aGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGcc                                                                  <  2:3630143/93‑1 (MQ=255)
                         aGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGcc                                                                  <  2:1665879/93‑1 (MQ=255)
                          gCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCg                                                                 >  1:1078962/1‑93 (MQ=255)
                            ggCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGAc                                                               >  2:5049877/1‑93 (MQ=255)
                            ggCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGAc                                                               >  2:4063776/1‑93 (MQ=255)
                             gCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACg                                                              >  2:1576842/1‑93 (MQ=255)
                             gCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACg                                                              >  1:2508567/1‑93 (MQ=255)
                              cTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGc                                                             <  1:1158772/93‑1 (MQ=255)
                               tttCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCTGACGCt                                                            >  2:6592937/1‑93 (MQ=255)
                               tttCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCt                                                            >  2:6296136/1‑93 (MQ=255)
                               tttCGTAAACCATTCAAAGGAAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGATTATGCATATAACCGATGTTGACGTAAGCGCAGACGCt                                                            >  1:1462021/1‑93 (MQ=255)
                                ttCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTc                                                           >  1:3395129/1‑93 (MQ=255)
                                 tCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTcc                                                          <  2:695096/93‑1 (MQ=255)
                                 tCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTcc                                                          <  2:6377110/93‑1 (MQ=255)
                                 tCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTcc                                                          >  2:1838769/1‑93 (MQ=255)
                                   gTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTcccc                                                        >  2:3434444/1‑93 (MQ=255)
                                    tAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCt                                                       <  1:5818299/93‑1 (MQ=255)
                                     aaaCCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTg                                                      >  1:2244263/1‑93 (MQ=255)
                                      aaCCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGa                                                     >  1:5780858/1‑93 (MQ=255)
                                        ccATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGc                                                   <  1:6635823/93‑1 (MQ=255)
                                        ccATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGc                                                   >  2:1158230/1‑93 (MQ=255)
                                         cATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCg                                                  <  2:5008460/93‑1 (MQ=255)
                                          aTTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCgg                                                 >  1:4613377/1‑93 (MQ=255)
                                           ttCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGc                                                >  1:491602/1‑93 (MQ=255)
                                           ttCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGc                                                >  2:177299/1‑93 (MQ=255)
                                           ttCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGc                                                >  1:5024111/1‑93 (MQ=255)
                                             cAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTCTGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAg                                              >  2:400705/1‑93 (MQ=255)
                                              aaaGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGc                                             >  2:2694400/1‑93 (MQ=255)
                                              aaaGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGc                                             >  2:4312771/1‑93 (MQ=255)
                                                aGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCt                                           <  1:4444403/93‑1 (MQ=255)
                                                aGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCt                                           >  2:5168135/1‑93 (MQ=255)
                                                 ggCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCtt                                          <  2:1898177/93‑1 (MQ=255)
                                                   cAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCa                                        >  1:3023143/1‑93 (MQ=255)
                                                    aGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAg                                       >  1:3214978/1‑93 (MQ=255)
                                                     gCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGa                                      <  2:2154196/93‑1 (MQ=255)
                                                        tGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATAc                                   >  2:1003805/1‑93 (MQ=255)
                                                        tGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATAc                                   >  2:1403708/1‑93 (MQ=255)
                                                        tGATAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATAc                                   >  2:3061200/1‑93 (MQ=255)
                                                           tAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAg                                >  2:2415870/1‑93 (MQ=255)
                                                           tAATCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAg                                <  1:3939336/93‑1 (MQ=255)
                                                           tAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAg                                <  1:3668837/93‑1 (MQ=255)
                                                            aaTCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGt                               >  2:1537984/1‑93 (MQ=255)
                                                             aTCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTa                              >  1:1061071/1‑93 (MQ=255)
                                                             aTCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTa                              <  1:4332319/93‑1 (MQ=255)
                                                             aTCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTa                              >  2:542255/1‑93 (MQ=255)
                                                              tCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTAc                             >  1:1546430/1‑93 (MQ=255)
                                                              tCTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTAc                             >  1:5523367/1‑93 (MQ=255)
                                                              tCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTAc                             <  2:5840937/93‑1 (MQ=255)
                                                               cTTTCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACa                            >  2:4333557/1‑93 (MQ=255)
                                                                tttCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACAt                           <  2:3744843/93‑1 (MQ=255)
                                                                tttCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACAt                           <  2:4854742/93‑1 (MQ=255)
                                                                tttCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCACTGAGCGGCAGCCTTCAGATACCAGTACAt                           >  2:1126091/1‑93 (MQ=255)
                                                                 ttCTCAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATc                          <  1:5315851/93‑1 (MQ=255)
                                                                    tcAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACCTCgct                       >  1:203817/1‑93 (MQ=255)
                                                                    tcAACGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCgct                       >  2:3363176/1‑93 (MQ=255)
                                                                    tcAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCgct                       <  1:5843606/93‑1 (MQ=255)
                                                                    tcAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCgct                       <  2:6115283/93‑1 (MQ=255)
                                                                    tcAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCgct                       <  2:2180480/93‑1 (MQ=255)
                                                                     cAACGCCTTGTCCGTGTTTATGCATATACCCGATGTTGACGTAAGCGCCGAGGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCgctt                      <  2:6194459/93‑1 (MQ=255)
                                                                       aCGCCTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCgcttgc                    <  2:3018199/93‑1 (MQ=255)
                                                                          ccTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTta                 <  1:2115294/93‑1 (MQ=255)
                                                                          ccTTGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTta                 >  1:3616579/1‑93 (MQ=255)
                                                                          ccTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTta                 >  2:3848354/1‑93 (MQ=255)
                                                                            ttGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTtata               <  1:3508535/93‑1 (MQ=255)
                                                                            ttGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTtata               >  2:1926317/1‑93 (MQ=255)
                                                                            ttGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTtata               <  1:299311/93‑1 (MQ=255)
                                                                             tGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAg              >  2:4432140/1‑93 (MQ=255)
                                                                             tGTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAg              <  2:319690/93‑1 (MQ=255)
                                                                              ttCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGt             <  1:4704956/92‑1 (MQ=255)
                                                                              gTCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGt             <  2:2957198/93‑1 (MQ=255)
                                                                               tCCGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTc            <  1:2592717/93‑1 (MQ=255)
                                                                                ccGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTcc           <  2:4690755/93‑1 (MQ=255)
                                                                                ccGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTcc           <  1:5094564/93‑1 (MQ=255)
                                                                                 cGTGTTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCt          >  2:1101695/1‑93 (MQ=255)
                                                                                  gtgtTTATGCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTg         >  1:6568657/1‑93 (MQ=255)
                                                                                         taCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCAcc  >  2:6578598/1‑93 (MQ=255)
                                                                                          |                                                                                           
ACCAGGCAGTGGCGTCATTGCATTCAGCGGCTTTCGTAAACCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACC  >  NC_000913/674933‑675114

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