Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 631,241 A→G intergenic (+19/‑164) cstA → / → ybdD carbon starvation protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913631,2410AG94.4% 320.6 / 11.0 110intergenic (+19/‑164)cstA/ybdDcarbon starvation protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  major base G (57/45);  minor base . (5/1);  total (63/47)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 123 total aligned reads displayed.

AAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAG  >  NC_000913/631152‑631332
                                                                                          |                                                                                           
aagaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAACTGCGCACACTACAGTCAGCGCGAGGGGGGGgtgt                                                                                          >  1:268640/1‑93 (MQ=255)
 agaaaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGtt                                                                                         >  2:2434208/1‑93 (MQ=255)
    aaCGCCGTATTAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGTCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGc                                                                                      <  1:795032/93‑1 (MQ=255)
    aaCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTggg                                                                                      >  2:2133097/1‑92 (MQ=255)
     aCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCg                                                                                     >  1:4023769/1‑93 (MQ=255)
      cGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGGGTTGGCGa                                                                                    >  1:2789677/1‑93 (MQ=255)
       gCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGaa                                                                                   >  1:843338/1‑93 (MQ=255)
       gCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGaa                                                                                   <  1:2900747/93‑1 (MQ=255)
        ccGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGaag                                                                                  >  2:1565719/1‑92 (MQ=255)
        ccGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAAt                                                                                  >  2:4231404/1‑93 (MQ=255)
         cGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGGTGTTGGCGAAt                                                                                  >  1:4700380/1‑93 (MQ=255)
         cGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATg                                                                                 <  1:2060183/93‑1 (MQ=255)
          gTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGt                                                                                >  1:693244/1‑93 (MQ=255)
           tATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGGTGTTGGGGAATGt                                                                                >  1:1065032/1‑93 (MQ=255)
           tATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGGGAATGtg                                                                               >  1:947278/1‑92 (MQ=255)
           tATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGtt                                                                               >  2:1843227/1‑93 (MQ=255)
           tATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGtt                                                                               <  1:2959374/93‑1 (MQ=255)
            aTGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAGAGTCAGAGTGAGGGGGGGGGGTTGGCGAATGTTg                                                                              >  2:4043904/1‑93 (MQ=255)
            aTGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTg                                                                              >  2:4495867/1‑93 (MQ=255)
               aaCCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGCGTTGGCGAATGTTGGCt                                                                           >  1:3440296/1‑93 (MQ=255)
               aaCCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAGTCAGAGGGAGGGGGGGGTGGTGGCGAATGTTGGCtt                                                                          >  1:221496/1‑93 (MQ=255)
                aCCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCtt                                                                          >  1:3475177/1‑93 (MQ=255)
                 ccAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTa                                                                         >  1:2449604/1‑93 (MQ=255)
                 ccAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGTGTTGGCGAATGTTGGCTTAg                                                                        >  1:2658214/1‑93 (MQ=255)
                  cAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAg                                                                        <  1:1532567/93‑1 (MQ=255)
                   aaTGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGt                                                                       >  2:1445440/1‑93 (MQ=255)
                   aaTGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGt                                                                       >  1:4190939/1‑93 (MQ=255)
                   aaTGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGt                                                                       >  1:2125006/1‑93 (MQ=255)
                    aTGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTg                                                                      >  1:702025/1‑93 (MQ=255)
                    aTGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTg                                                                      <  2:30969/93‑1 (MQ=255)
                     tGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGGGTTGGCGAATGTTGGCTTAGTGc                                                                     >  2:3257041/1‑93 (MQ=255)
                      gCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGcc                                                                    <  1:1541107/93‑1 (MQ=255)
                       ccGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGccc                                                                   <  1:1709129/93‑1 (MQ=255)
                         ggAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAg                                                                 <  2:727771/93‑1 (MQ=255)
                         ggAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAgg                                                                >  2:4669163/1‑93 (MQ=255)
                          gAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAgg                                                                >  1:3216496/1‑93 (MQ=255)
                          gAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAgg                                                                >  2:2773877/1‑93 (MQ=255)
                           aaaaTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAgg                                                                <  1:615149/93‑1 (MQ=255)
                           aaaaTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAggg                                                               >  1:791722/1‑93 (MQ=255)
                           aaaaTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAggg                                                               >  1:475449/1‑93 (MQ=255)
                           aaaaTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAggg                                                               <  1:4362098/93‑1 (MQ=255)
                            aaaTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGt                                                              >  2:2328360/1‑93 (MQ=255)
                               tGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTcc                                                           <  2:1103728/93‑1 (MQ=255)
                               tGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTcc                                                           <  2:390830/93‑1 (MQ=255)
                                gTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTccc                                                          >  2:520190/1‑93 (MQ=255)
                                 tCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCct                                                         <  2:1187892/93‑1 (MQ=255)
                                  cGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCctc                                                        >  1:4547980/1‑93 (MQ=255)
                                   gaggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCctct                                                       <  2:1167703/93‑1 (MQ=255)
                                    aggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCtctc                                                      <  1:4467967/93‑1 (MQ=255)
                                    aggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCtctc                                                      <  2:298344/93‑1 (MQ=255)
                                     ggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGGGGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCtctc                                                      >  2:3505412/1‑93 (MQ=255)
                                     ggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCa                                                     >  2:2056499/1‑93 (MQ=255)
                                      gagaTCGTGGCGCAGGCAAAACGCGCACACTAAAGTCAGAGGGAGGGGGGGGCGTTGGCGAATGTTTGCTTAGTGCCCAGGGTTCCCTCTCAc                                                    >  2:2595942/1‑93 (MQ=255)
                                         aTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCt                                                 <  1:2721296/93‑1 (MQ=255)
                                          tCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTa                                                <  1:3167409/93‑1 (MQ=255)
                                            gTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAAc                                              >  1:3532201/1‑93 (MQ=255)
                                            gTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAAc                                              >  2:2632537/1‑93 (MQ=255)
                                              ggCGCAGGCAAAAGGCGCACAGTAAAGTGAGAGTGAGGGGGCGGTGTAGGAGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAAccc                                            <  2:684819/93‑1 (MQ=255)
                                              ggCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAAccc                                            <  1:1362745/93‑1 (MQ=255)
                                              ggCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAAccc                                            >  2:3071212/1‑93 (MQ=255)
                                              ggCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAAccc                                            <  2:684818/93‑1 (MQ=255)
                                               gcgcAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCct                                           >  2:1753958/1‑93 (MQ=255)
                                                 gcAGTCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCctct                                         <  1:1895336/93‑1 (MQ=255)
                                                        aaaGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGt                                  <  2:771529/93‑1 (MQ=255)
                                                        aaaGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGt                                  >  1:1799859/1‑93 (MQ=255)
                                                          aGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTgg                                <  1:3336720/93‑1 (MQ=255)
                                                          aGGCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTgg                                <  1:691579/93‑1 (MQ=255)
                                                          aGGCGCACACTAAAGTCAGAGTGAGGGGGGGGCGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCGCGGTgg                                >  1:618672/1‑93 (MQ=255)
                                                           ggCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTggg                               >  1:519000/1‑93 (MQ=255)
                                                           ggCGCACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTggg                               <  1:4398434/93‑1 (MQ=255)
                                                             cgcACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGc                             <  2:1445809/93‑1 (MQ=255)
                                                              gcACACTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCg                            >  2:3364457/1‑93 (MQ=255)
                                                                acacTAAAGTCAGAGTGAGGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGa                           <  1:3179743/93‑1 (MQ=255)
                                                                acacTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAg                          >  2:1454388/1‑93 (MQ=255)
                                                                acacTAAAGTCAGAGTGAGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAgg                         >  1:292392/1‑93 (MQ=255)
                                                                 cacaAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAgg                         <  1:1128321/93‑1 (MQ=255)
                                                                  acTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAggg                        >  2:3388338/1‑93 (MQ=255)
                                                                  acTAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAggg                        <  1:1051149/93‑1 (MQ=255)
                                                                    tAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAggggg                      >  2:407140/1‑92 (MQ=255)
                                                                    tAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGa                      >  1:1194318/1‑93 (MQ=255)
                                                                    tAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGa                      >  2:4229936/1‑93 (MQ=255)
                                                                    tAAAGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGa                      >  1:983080/1‑93 (MQ=255)
                                                                     aaaGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGAc                     <  2:1058529/93‑1 (MQ=255)
                                                                      aaGTCAGAGTGAGGGGGGGGTGTTGGCGAATGTAGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCACCGGTGGGGCGAGGGGAca                    <  1:4202354/93‑2 (MQ=38)
                                                                         gcAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGAc                 <  1:3034662/92‑1 (MQ=255)
                                                                          cAGAGTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGAcc                >  1:4080869/1‑93 (MQ=255)
                                                                            gagTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGa              >  1:2353738/1‑93 (MQ=255)
                                                                             agTGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAg             <  1:184533/93‑1 (MQ=255)
                                                                             agTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAg             <  1:3728342/93‑1 (MQ=255)
                                                                               tGAGGGGGGGGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAgcg           >  1:411780/1‑93 (MQ=255)
                                                                                    ggggggGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGtt       <  1:1716121/93‑1 (MQ=255)
                                                                                    ggggggGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGtt       <  2:3251623/93‑1 (MQ=255)
                                                                                    ggggggGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGtt       <  1:2235288/93‑1 (MQ=255)
                                                                                    ggggggGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGtt       <  1:3995185/93‑1 (MQ=255)
                                                                                      ggggGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGa     <  1:3995364/93‑1 (MQ=255)
                                                                                      ggggGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGa     <  1:4517900/93‑1 (MQ=255)
                                                                                      ggggGTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGa     >  1:4740122/1‑93 (MQ=255)
                                                                                       ggggTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGAt    <  1:2538485/93‑1 (MQ=255)
                                                                                       ggggTGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGAt    >  2:635604/1‑93 (MQ=255)
                                                                                         gATGTTGGCGAATGTTGGCTTAGTGCCAAGGGTTCCCTCTCACCCTAACCCTCCCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAg  <  2:125825/93‑1 (MQ=255)
                                                                                          |                                                                                           
AAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAG  >  NC_000913/631152‑631332

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