Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,599,241 T→A 8.4% K174* (AAA→TAA) ‡ dapA ← dihydrodipicolinate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,599,2410TA 8.4% 281.2 / 9.6 101K174* (AAA→TAA) ‡dapAdihydrodipicolinate synthase
Reads supporting (aligned to +/- strand):  ref base T (43/49);  new base A (7/2);  total (50/51)
Fisher's exact test for biased strand distribution p-value = 9.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

CTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAGGTCAGTAT  >  NC_000913/2599094‑2599387
                                                                                                                                                   |                                                                                                                                                  
ctgctgCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGa                                                                                                                                                  <  2:352394/150‑1 (MQ=255)
 tgctgcCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCGTTGAt                                                                                                                                                 >  2:584232/1‑150 (MQ=255)
  gctgcCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATc                                                                                                                                                >  2:582508/1‑150 (MQ=255)
  gctgcCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATc                                                                                                                                                >  2:610551/1‑150 (MQ=255)
   ctgcCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCt                                                                                                                                               >  2:1023954/1‑150 (MQ=255)
       cAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGtt                                                                                                                                           <  1:1551508/150‑1 (MQ=255)
       cAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGtt                                                                                                                                           <  2:1624558/150‑1 (MQ=255)
        aGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGTAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCATAACAAAATCATCTGAAACCAGCTCTTTGATCTGGttt                                                                                                                                          >  2:852742/1‑150 (MQ=255)
        aGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGttt                                                                                                                                          <  1:610551/150‑1 (MQ=255)
            tGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTacac                                                                                                                                      >  1:607703/1‑150 (MQ=255)
            tGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTacac                                                                                                                                      <  1:642339/150‑1 (MQ=255)
             gCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGGAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCAAAGAGCTGGtttcagat                                                                                                                                     <  1:702411/150‑6 (MQ=255)
              cacaTCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACAcgc                                                                                                                                    <  1:995153/150‑1 (MQ=255)
              cacaTCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACAAGCTCTTTGATCTGGTTTacactc                                                                                                                                    >  2:1216209/1‑148 (MQ=255)
               acaTCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATCACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACAcgcg                                                                                                                                   <  2:816647/150‑1 (MQ=255)
               acaTCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACAcgcg                                                                                                                                   <  2:371152/150‑1 (MQ=255)
               acaTCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACAcgcg                                                                                                                                   <  2:607703/150‑1 (MQ=255)
               acaTCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACAcgcg                                                                                                                                   >  1:1180375/1‑150 (MQ=255)
                     gggCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCCGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCCGCTCTGTGATCTGGTTATCACGCGTTAAGt                                                                                                                             >  2:1558171/1‑150 (MQ=255)
                      ggCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGtt                                                                                                                            <  1:657854/150‑1 (MQ=255)
                       gCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGAATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGACCAAAATCATCTGAAACCAGCTCTTTGATCAGGTGTTCACGCGTTAAGTTc                                                                                                                           >  1:260830/1‑150 (MQ=255)
                        ccATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCCTGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGATcc                                                                                                                          <  2:1040254/150‑1 (MQ=255)
                           tatCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCACTg                                                                                                                       >  2:596676/1‑150 (MQ=255)
                               aCGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGc                                                                                                                   >  2:634588/1‑150 (MQ=255)
                                cgcgCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGcc                                                                                                                  >  1:464081/1‑150 (MQ=255)
                                  cgcTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCctc                                                                                                                <  1:642022/150‑1 (MQ=255)
                                    cTGCGACGTTAGCCGTAACGGAAATAACCCCATGCCCGCCCAAATGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCctcat                                                                                                              >  2:183289/1‑148 (MQ=255)
                                    cTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCtt                                                                                                              <  1:511454/150‑1 (MQ=255)
                                     tGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCttt                                                                                                             <  2:470564/150‑1 (MQ=255)
                                       cGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGCACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGa                                                                                                           >  1:710846/1‑150 (MQ=255)
                                       cGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGa                                                                                                           <  1:504204/150‑1 (MQ=255)
                                         aCGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGAtt                                                                                                         <  2:1291987/150‑1 (MQ=255)
                                           gTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTcc                                                                                                       >  2:1346327/1‑150 (MQ=255)
                                           gTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTcc                                                                                                       <  2:2776/150‑1 (MQ=255)
                                           gTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGGTAAGTTCCCTGTTGCCTCTTTGATTcc                                                                                                       >  2:66461/1‑150 (MQ=255)
                                             tAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGa                                                                                                     <  1:596676/150‑1 (MQ=255)
                                              aGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCATCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGat                                                                                                    <  2:302335/150‑1 (MQ=255)
                                                  gTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGataata                                                                                                <  1:966254/150‑1 (MQ=255)
                                                      cGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATAtttt                                                                                            >  1:1324318/1‑150 (MQ=255)
                                                       ggAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATAttttt                                                                                           <  2:464081/150‑1 (MQ=255)
                                                       ggAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATAttttt                                                                                           <  2:1296255/150‑1 (MQ=255)
                                                        taaaTAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCATCTCTTTGATCTGGTTTACACGCGTTTAGTTCCCTGTTGCCTCTTTGATTCCGATAATAtttttt                                                                                          <  1:695908/149‑1 (MQ=255)
                                                        gAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCTCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATAtttttt                                                                                          <  1:94624/150‑1 (MQ=255)
                                                        gAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATAtttttt                                                                                          <  2:150331/150‑1 (MQ=255)
                                                            tAACCCCATGACCGCCCAATTGCATGAAGTCAAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCACTGTTGCCTCTTTGATTCCGATAATATTTTTTTCtt                                                                                      >  1:134622/1‑150 (MQ=255)
                                                               ccccATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCg                                                                                   <  1:192138/150‑1 (MQ=255)
                                                               ccccATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCg                                                                                   <  1:625036/150‑1 (MQ=255)
                                                                cccATGACCGCCCAATTGCATGAAGTGCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGc                                                                                  <  1:1097958/150‑1 (MQ=255)
                                                                   aTGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAg                                                                               <  1:1031215/150‑1 (MQ=255)
                                                                          cccAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGccc                                                                        >  2:426775/1‑150 (MQ=255)
                                                                          cccAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGccc                                                                        >  2:659329/1‑150 (MQ=255)
                                                                          cccAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGccc                                                                        <  2:452225/150‑1 (MQ=255)
                                                                            cAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCAc                                                                      <  1:1023954/150‑1 (MQ=255)
                                                                                 gCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGttt                                                                 >  1:1232802/1‑150 (MQ=255)
                                                                                 gCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAAACATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAACTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGttt                                                                 >  2:462226/1‑150 (MQ=255)
                                                                                   aTGAATTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAATATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTcc                                                               <  2:1505457/150‑1 (MQ=255)
                                                                                      aaGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCggg                                                            <  2:893006/150‑1 (MQ=255)
                                                                                       aGTCCAGCGCCCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTAGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTATCGCCAGACGGCCCACCGTTTCCGGGa                                                           >  2:891030/1‑150 (MQ=255)
                                                                                        gTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCggg                                                            >  1:1290517/1‑148 (MQ=255)
                                                                                        gTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCggg                                                            <  2:1290517/148‑1 (MQ=255)
                                                                                          ccAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCa                                                        >  1:1191380/1‑150 (MQ=255)
                                                                                                  cTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGAGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAg                                                >  2:733545/1‑150 (MQ=255)
                                                                                                    cGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCGTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCCCCGTTTCCGGGCGCAGATCGCAGcc                                              >  2:1254680/1‑150 (MQ=255)
                                                                                                      catcatCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAg                                            <  1:582508/150‑1 (MQ=255)
                                                                                                          atcGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTAGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCAAGACGGCACACCGTTTCAGGGAGCAGATAGCAGCCAGTACg                                        >  1:1233386/1‑150 (MQ=255)
                                                                                                           tcGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACgg                                       <  1:66461/150‑1 (MQ=255)
                                                                                                            cgccgcTCAGCAGCACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACggg                                      >  1:1339815/1‑150 (MQ=255)
                                                                                                                gcTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGATAAGTTACCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCTCCAGACGGCCCACAGTTTCAGGGAGCAGCTCCAACCAAGTACGGGACgg                                  >  2:336819/1‑150 (MQ=255)
                                                                                                                gcTCACCAGAACAAAATCATCTGAAACCAGCCCTTTGATCTGGTGTACACGCGTGAAGTACCCTGGTGCCTCTTTGATGCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTGCCGGGAGCAGATCGCAGCCAGTACGGGACgg                                  <  1:141110/150‑1 (MQ=255)
                                                                                                                   cagcagAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGcac                               >  2:437773/1‑150 (MQ=255)
                                                                                                                   cagcagAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGAGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGcac                               >  2:473182/1‑150 (MQ=255)
                                                                                                                     gcagAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACAt                             <  2:727335/150‑1 (MQ=255)
                                                                                                                     gcagAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACAt                             >  1:34709/1‑150 (MQ=255)
                                                                                                                        gAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCCGCCAGTACGGGACGGCACATtat                          >  1:1552642/1‑150 (MQ=255)
                                                                                                                               atcatcTGAAACCAGCTCTTTCATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAAt                   <  1:1120623/150‑1 (MQ=255)
                                                                                                                                tgatcTGAATCCAGTTCTTTGATGTGGTTTACTGGGGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCTGCCAGTACGGGACGGCACATTATACAGAttt                  <  1:1698807/148‑1 (MQ=255)
                                                                                                                                tcatcTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAAtt                  <  1:426775/150‑1 (MQ=255)
                                                                                                                                cgcgtgtaaaCCAGATCAAAGAGCTGGTTTACACGCGTTAAGTTCCCTGTTTCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGCACGGCACATTATACAGAAtt                  >  2:50906/8‑150 (MQ=255)
                                                                                                                                 gcgtgtaaaCCAGATCAAAGAGCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGCCGGCACATTATACAGAAttt                 >  1:1631496/7‑150 (MQ=255)
                                                                                                                                  cgtgtaaaCCAGATCAAAGAGCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTg                >  2:1497525/6‑150 (MQ=255)
                                                                                                                                  cgtgtaaaCCAGATCAAAGAGCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTg                >  2:457540/6‑150 (MQ=255)
                                                                                                                                  atcTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTg                >  2:216246/1‑150 (MQ=255)
                                                                                                                                   tcTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGc               >  1:1415632/1‑150 (MQ=255)
                                                                                                                                   tcTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGAGAATATTTTTTACTTTCGCCCGACGGCCCCCCGTGTCCGGGAGCAGATCGCGGCCAGAACGGGACGGCACATTAGACAGACGTTGc               >  2:722679/1‑150 (MQ=255)
                                                                                                                                   gtgtaaaCCAGATCAAAGAGCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGc               >  2:238713/5‑150 (MQ=255)
                                                                                                                                    cTGAAACCCGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCGCACCGTTTCCGGGAGCAGATCGCCGCCAGTACGCGACGGCCCATTATAGAGAAGTTGCg              >  2:861933/1‑150 (MQ=255)
                                                                                                                                    cTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCg              <  1:854588/150‑1 (MQ=255)
                                                                                                                                    cTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCAAGTACGGGACGGCACATTATACAGAATTTGCg              >  1:1179605/1‑150 (MQ=255)
                                                                                                                                    cTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGGTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAAGTTGCg              >  2:1021482/1‑150 (MQ=255)
                                                                                                                                     gtaaaCCAGATCAAAGAGCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGTCACATTATACAGAATTTGCgg             <  1:1497525/148‑1 (MQ=255)
                                                                                                                                      gAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGc            >  1:624441/1‑150 (MQ=255)
                                                                                                                                       ataCCATCTCTTTGATCTGTTTTACACGCGTTAAGTTCACTGTTGCCTCTTTTCTTCAGAAAATATTTTTAACTTTAGCCAGACTTCCCACAGTTTCCGGTAGCAGATCGCAGACAGTACGGGACGGCACATTATACAGAATTTGCGGCa           <  2:1018318/148‑1 (MQ=255)
                                                                                                                                       aaaCCAGATCAAAGAGCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGAAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCa           >  1:1700333/1‑150 (MQ=255)
                                                                                                                                        aaCCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAg          <  1:854338/150‑1 (MQ=255)
                                                                                                                                         accttctctTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTTCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCTCCGTTTCCGCGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAgg         <  1:447152/145‑1 (MQ=255)
                                                                                                                                           cAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAGGTc       <  1:1062418/150‑1 (MQ=255)
                                                                                                                                             gCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAGGTCAg     <  2:624441/150‑1 (MQ=255)
                                                                                                                                               tctTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCTATAATATTTTTTACTTTCGCCAGACGGCCCACCTTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATAAAGAATTTGCGGCAGGTCAGTa   <  1:852742/150‑1 (MQ=255)
                                                                                                                                               tctTTGATCTGGTTTACACGAGTTAATTTCCCTGTTTCATCTTTGTTTCCGATTTTATTTTTTACTTTCTCCAGACGGCCCACCGTTTCCGGGAGCAGATAGCAGCCAGTACGGGACGGCACATTATACTGAATTTGCGGCATGTCAGTa   <  2:1112393/150‑1 (MQ=255)
                                                                                                                                                ctTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAGGTCAGTAt  >  2:1168032/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAGGTCAGTAT  >  NC_000913/2599094‑2599387

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: