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.

GCCCCGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCC  >  NC_000913/2912327‑2912525
                                                                                                    |                                                                                                   
gCCCCGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTaa                                                                                                     >  1:2889285/1‑100 (MQ=255)
 ccccGATTCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTaaa                                                                                                    >  2:2288014/1‑100 (MQ=255)
 ccccGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTaaa                                                                                                    >  1:2997500/1‑100 (MQ=255)
 ccccGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTaaa                                                                                                    <  1:1882656/100‑1 (MQ=255)
 ccccGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑ACCCGCAGGCAAATCAACGACATCACCTTTCGGCGTaaa                                                                                                    <  2:2400913/100‑1 (MQ=255)
  cccGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAg                                                                                                   >  1:1372702/1‑100 (MQ=255)
   ccGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCTAATTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTagag                                                                                                   >  2:1298310/1‑99 (MQ=255)
   ccGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGa                                                                                                  >  1:1811076/1‑100 (MQ=255)
    cGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGAc                                                                                                 <  2:2482207/100‑1 (MQ=255)
       tGCAGCGGTGTCTGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTa                                                                                              <  1:976163/100‑1 (MQ=255)
         cAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTaca                                                                                            >  1:163771/1‑100 (MQ=255)
          aGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTacac                                                                                           >  1:1832906/1‑100 (MQ=255)
             ggTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCg                                                                                        <  2:1476705/100‑1 (MQ=255)
               tgtCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGgt                                                                                      <  1:2833355/100‑1 (MQ=255)
                gtCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGgtc                                                                                     <  1:2554134/100‑1 (MQ=255)
                gtCCGACATCACTGTGGATGTGGTAAGCGAAATCCAGCGGCGTTGA‑TCCCGCCGGCAAATCAACGACATAACCTTTCGGCGTAAAGACGTACACCCGgtc                                                                                     >  2:3139582/1‑100 (MQ=255)
                 tCCGATATCACTGTGGTTGTGGTTAGCTATGTTCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGgtcg                                                                                    <  1:1743955/100‑1 (MQ=255)
                 tCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGgtcg                                                                                    <  2:432118/100‑1 (MQ=255)
                 tCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGgtcg                                                                                    <  2:3569368/100‑1 (MQ=255)
                 tCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGgtcg                                                                                    <  2:300153/100‑1 (MQ=255)
                  ccGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGgtcct                                                                                   <  1:2651361/100‑3 (MQ=255)
                  ccGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGgtcgt                                                                                   >  2:3545568/1‑100 (MQ=255)
                    gACATCACTGTGGATGTGGTAAGCGAAGTCTAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGgtcgt                                                                                   <  1:587703/98‑1 (MQ=255)
                     aCATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCaa                                                                                <  2:1658615/100‑1 (MQ=255)
                       aTCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATTCCCGCAGGCAAATCAACGACATCACCCTTCGGCGTAAAGACGTACACCCGGTCGTCaa                                                                                <  2:833167/99‑1 (MQ=255)
                       aTCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGAGGTACACCCGGTCGTCaaa                                                                               <  2:2959639/99‑1 (MQ=255)
                       aTCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAg                                                                              <  2:3845963/100‑1 (MQ=255)
                        tCACTGTGGATGTGGTACGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGGAAAGACGTACACCCGGTCGTCAAAg                                                                              >  2:3271334/1‑99 (MQ=255)
                        tCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACCACATCACCTTTCGGCGTAAAGACGTAAACCCGGTCGTCAAAg                                                                              <  2:2536081/99‑1 (MQ=255)
                         cACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGAc                                                                            <  1:445685/100‑1 (MQ=255)
                          aCTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGAcc                                                                           >  2:471454/1‑100 (MQ=255)
                          aCTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGAcc                                                                           <  2:841443/100‑1 (MQ=255)
                             gtgGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGa                                                                        >  2:1031513/1‑100 (MQ=255)
                              tgGATGTGTTAAGCGAAGTCCAGCTGCGTTGA‑TTTCGCAGGTAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGAc                                                                       <  1:1249862/100‑1 (MQ=255)
                               ggATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACt                                                                      <  1:1984001/100‑1 (MQ=255)
                                gATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTa                                                                     <  2:95776/100‑1 (MQ=255)
                                gATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTa                                                                     <  2:3709400/100‑1 (MQ=255)
                                gATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTa                                                                     <  1:2426416/100‑1 (MQ=255)
                                  tgtgGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACg                                                                   <  1:3609276/100‑1 (MQ=255)
                                      gTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACt                                                               >  2:575098/1‑100 (MQ=255)
                                        aaGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTtc                                                             >  1:2305475/1‑100 (MQ=255)
                                         aGCGAAGTCCAGCGGTGTTGA‑TCCCGCAGGCAAATCAACGATATTACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTtcg                                                            <  1:783750/100‑1 (MQ=255)
                                         aGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTtcg                                                            >  1:1279425/1‑100 (MQ=255)
                                            gAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTtcgtcn                                                         <  1:3809401/100‑2 (MQ=255)
                                              aGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAg                                                       >  2:3689265/1‑100 (MQ=255)
                                               gTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGc                                                      >  1:3163734/1‑100 (MQ=255)
                                                tCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCa                                                     <  1:1125716/100‑1 (MQ=255)
                                                tCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTGCACCCGGTCGTCAAAGGCCTGACTACGTACTTCGTCGAGCa                                                     <  2:1133669/100‑1 (MQ=255)
                                                tCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCa                                                     >  1:2464723/1‑100 (MQ=255)
                                                tCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCa                                                     <  1:1768077/100‑1 (MQ=255)
                                                tCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCa                                                     <  2:2715917/100‑1 (MQ=255)
                                                tCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCa                                                     >  1:3211094/1‑100 (MQ=255)
                                                   aGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCAttt                                                  <  2:3881172/100‑1 (MQ=255)
                                                   aGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCAttt                                                  >  2:1828592/1‑100 (MQ=255)
                                                       gcgTTGA‑TCCCGCAGGCAAATCAACGAAATCACCTTTCGGCGTAAAAAAGTACACCCAGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGcc                                              >  2:1937379/1‑100 (MQ=255)
                                                        cgTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCg                                             <  1:3155794/100‑1 (MQ=255)
                                                        cgTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCg                                             >  1:1399754/1‑100 (MQ=255)
                                                         gTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCgg                                            <  2:3063686/100‑1 (MQ=255)
                                                          ttGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGa                                           >  1:2743816/1‑100 (MQ=255)
                                                          ttGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGa                                           >  1:3890193/1‑100 (MQ=255)
                                                          ttGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGa                                           >  1:1582440/1‑100 (MQ=255)
                                                          ttGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGa                                           <  1:1469249/100‑1 (MQ=255)
                                                             a‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATc                                        <  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)
                                                                                                    |                                                                                                   
GCCCCGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGA‑TCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGCAGCC  >  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: ‑