Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,912,426 A→G 21.1% F409S (TTT→TCT)  relA ← (p)ppGpp synthetase I/GTP pyrophosphokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,912,4260AG21.1% 266.0 / 73.6 142F409S (TTT→TCT) relA(p)ppGpp synthetase I/GTP pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base A (55/57);  new base G (13/17);  total (68/74)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 154 total aligned reads displayed.

GCCCCGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCC  >  NC_000913/2912327‑2912525
                                                                                                    |                                                                                                   
gCCCCGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTaa                                                                                                     >  1:2889285/1‑100 (MQ=255)
 ccccGATTCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTaaa                                                                                                    >  2:2288014/1‑100 (MQ=255)
 ccccGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTaaa                                                                                                    >  1:2997500/1‑100 (MQ=255)
 ccccGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTaaa                                                                                                    <  1:1882656/100‑1 (MQ=255)
 ccccGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGAACCCGCAGGCAAATCAACGACATCACCTTTCGGCGTaaa                                                                                                    <  2:2400913/100‑1 (MQ=255)
  cccGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAg                                                                                                   >  1:1372702/1‑100 (MQ=255)
   ccGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCTAATTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTagag                                                                                                   >  2:1298310/1‑99 (MQ=255)
   ccGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGa                                                                                                  >  1:1811076/1‑100 (MQ=255)
    cGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGAc                                                                                                 <  2:2482207/100‑1 (MQ=255)
       tGCAGCGGTGTCTGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTa                                                                                              <  1:976163/100‑1 (MQ=255)
         cAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTaca                                                                                            >  1:163771/1‑100 (MQ=255)
          aGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTacac                                                                                           >  1:1832906/1‑100 (MQ=255)
             ggTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCg                                                                                        <  2:1476705/100‑1 (MQ=255)
               tgtCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGgt                                                                                      <  1:2833355/100‑1 (MQ=255)
                gtCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGgtc                                                                                     <  1:2554134/100‑1 (MQ=255)
                gtCCGACATCACTGTGGATGTGGTAAGCGAAATCCAGCGGCGTTGATCCCGCCGGCAAATCAACGACATAACCTTTCGGCGTAAAGACGTACACCCGgtc                                                                                     >  2:3139582/1‑100 (MQ=255)
                 tCCGATATCACTGTGGTTGTGGTTAGCTATGTTCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGgtcg                                                                                    <  1:1743955/100‑1 (MQ=255)
                 tCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGgtcg                                                                                    <  2:432118/100‑1 (MQ=255)
                 tCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGgtcg                                                                                    <  2:3569368/100‑1 (MQ=255)
                 tCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGgtcg                                                                                    <  2:300153/100‑1 (MQ=255)
                  ccGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGgtcct                                                                                   <  1:2651361/100‑3 (MQ=255)
                  ccGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGgtcgt                                                                                   >  2:3545568/1‑100 (MQ=255)
                    gACATCACTGTGGATGTGGTAAGCGAAGTCTAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGgtcgt                                                                                   <  1:587703/98‑1 (MQ=255)
                     aCATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCaa                                                                                <  2:1658615/100‑1 (MQ=255)
                       aTCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATTCCCGCAGGCAAATCAACGACATCACCCTTCGGCGTAAAGACGTACACCCGGTCGTCaa                                                                                <  2:833167/99‑1 (MQ=255)
                       aTCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGAGGTACACCCGGTCGTCaaa                                                                               <  2:2959639/99‑1 (MQ=255)
                       aTCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAg                                                                              <  2:3845963/100‑1 (MQ=255)
                        tCACTGTGGATGTGGTACGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGGAAAGACGTACACCCGGTCGTCAAAg                                                                              >  2:3271334/1‑99 (MQ=255)
                        tCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACCACATCACCTTTCGGCGTAAAGACGTAAACCCGGTCGTCAAAg                                                                              <  2:2536081/99‑1 (MQ=255)
                         cACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGAc                                                                            <  1:445685/100‑1 (MQ=255)
                          aCTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGAcc                                                                           >  2:471454/1‑100 (MQ=255)
                          aCTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGAcc                                                                           <  2:841443/100‑1 (MQ=255)
                             gtgGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGa                                                                        >  2:1031513/1‑100 (MQ=255)
                              tgGATGTGTTAAGCGAAGTCCAGCTGCGTTGATTTCGCAGGTAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGAc                                                                       <  1:1249862/100‑1 (MQ=255)
                               ggATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACt                                                                      <  1:1984001/100‑1 (MQ=255)
                                gATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTa                                                                     <  2:95776/100‑1 (MQ=255)
                                gATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTa                                                                     <  2:3709400/100‑1 (MQ=255)
                                gATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTa                                                                     <  1:2426416/100‑1 (MQ=255)
                                  tgtgGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACg                                                                   <  1:3609276/100‑1 (MQ=255)
                                      gTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACt                                                               >  2:575098/1‑100 (MQ=255)
                                        aaGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTtc                                                             >  1:2305475/1‑100 (MQ=255)
                                         aGCGAAGTCCAGCGGTGTTGATCCCGCAGGCAAATCAACGATATTACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTtcg                                                            <  1:783750/100‑1 (MQ=255)
                                         aGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTtcg                                                            >  1:1279425/1‑100 (MQ=255)
                                            gAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTtcgtcn                                                         <  1:3809401/100‑2 (MQ=255)
                                              aGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAg                                                       >  2:3689265/1‑100 (MQ=255)
                                               gTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGc                                                      >  1:3163734/1‑100 (MQ=255)
                                                tCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCa                                                     <  1:1125716/100‑1 (MQ=255)
                                                tCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTGCACCCGGTCGTCAAAGGCCTGACTACGTACTTCGTCGAGCa                                                     <  2:1133669/100‑1 (MQ=255)
                                                tCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCa                                                     >  1:2464723/1‑100 (MQ=255)
                                                tCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCa                                                     <  1:1768077/100‑1 (MQ=255)
                                                tCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCa                                                     <  2:2715917/100‑1 (MQ=255)
                                                tCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCa                                                     >  1:3211094/1‑100 (MQ=255)
                                                   aGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCAttt                                                  <  2:3881172/100‑1 (MQ=255)
                                                   aGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCAttt                                                  >  2:1828592/1‑100 (MQ=255)
                                                       gcgTTGATCCCGCAGGCAAATCAACGAAATCACCTTTCGGCGTAAAAAAGTACACCCAGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGcc                                              >  2:1937379/1‑100 (MQ=255)
                                                        cgTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCg                                             <  1:3155794/100‑1 (MQ=255)
                                                        cgTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCg                                             >  1:1399754/1‑100 (MQ=255)
                                                         gTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCgg                                            <  2:3063686/100‑1 (MQ=255)
                                                          ttGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGa                                           >  1:2743816/1‑100 (MQ=255)
                                                          ttGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGa                                           >  1:3890193/1‑100 (MQ=255)
                                                          ttGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGa                                           >  1:1582440/1‑100 (MQ=255)
                                                          ttGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGa                                           <  1:1469249/100‑1 (MQ=255)
                                                             aTCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATc                                        <  2:1645742/100‑1 (MQ=255)
                                                                cccGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAg                                      >  2:3401965/1‑100 (MQ=255)
                                                                cccGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCATCGAGCATTTCGCCGGAATCAg                                      >  1:1082680/1‑100 (MQ=255)
                                                                  cGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCa                                       >  1:683937/1‑97 (MQ=255)
                                                                  cGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGcc                                    >  2:64221/1‑100 (MQ=255)
                                                                   gCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCa                                   <  2:2193237/100‑1 (MQ=255)
                                                                    cAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCAt                                  >  2:3483897/1‑100 (MQ=255)
                                                                     aGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCAtc                                 <  1:2425293/100‑1 (MQ=255)
                                                                      ggCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCAtct                                >  1:189817/1‑100 (MQ=255)
                                                                       gCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCAtctc                               <  1:3302737/100‑1 (MQ=255)
                                                                        cAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGAGCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATctct                              >  1:1508694/1‑100 (MQ=255)
                                                                        cAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATctct                              >  1:3127115/1‑100 (MQ=255)
                                                                        cAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATctct                              >  2:2809805/1‑100 (MQ=255)
                                                                           aTCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTcc                           <  1:3618024/100‑1 (MQ=255)
                                                                           aTCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTcc                           <  1:2062060/100‑1 (MQ=255)
                                                                            tCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTcc                           <  1:3737489/99‑1 (MQ=255)
                                                                              aaCGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGc                        >  1:2159114/1‑100 (MQ=255)
                                                                              aaCGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGc                        <  2:481376/100‑1 (MQ=255)
                                                                              aaCGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGc                        <  1:3402003/100‑1 (MQ=255)
                                                                               aCGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGcc                       <  1:3009065/100‑1 (MQ=255)
                                                                                cGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCa                      <  2:3203394/100‑1 (MQ=255)
                                                                                  aCATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACg                    >  2:1145597/1‑100 (MQ=255)
                                                                                  aCATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACg                    <  2:1424282/100‑1 (MQ=255)
                                                                                   cATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGc                   <  2:507043/100‑1 (MQ=255)
                                                                                    aTCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCa                  <  2:3032609/100‑1 (MQ=255)
                                                                                     tCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCaa                 >  1:2247999/1‑100 (MQ=255)
                                                                                      cACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCAGTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAAt                >  1:2750068/1‑100 (MQ=255)
                                                                                       aCCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCATCGAGCATTTCACCGAAATCAGCCATCTCTTCCTGCCACGCaa                 >  1:3091219/1‑98 (MQ=255)
                                                                                       aCCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATc               >  1:2186720/1‑100 (MQ=255)
                                                                                         cTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAg             >  2:2232991/1‑100 (MQ=255)
                                                                                          tttCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGt            <  2:1311944/100‑1 (MQ=255)
                                                                                              ggCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTAc        >  2:2766886/1‑100 (MQ=255)
                                                                                              ggCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTAc        >  2:3811466/1‑100 (MQ=255)
                                                                                              ggCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTAc        >  1:1612105/1‑100 (MQ=255)
                                                                                                 gTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCa     >  2:1957988/1‑100 (MQ=255)
                                                                                                  tAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAg    >  2:216228/1‑100 (MQ=255)
                                                                                                   aaaGACGTACACCCGGTCGTCAAAGACCTGACTTCGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGc   <  1:367863/100‑1 (MQ=255)
                                                                                                    aaGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGcc  >  1:853931/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
GCCCCGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCC  >  NC_000913/2912327‑2912525

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