Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,180,694 A→G 6.4% T92A (ACG→GCG)  ygiC → putative enzyme with ATPase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,180,6940AG 6.4% 273.6 / 7.1 93T92A (ACG→GCG) ygiCputative enzyme with ATPase activity
Reads supporting (aligned to +/- strand):  ref base A (42/45);  new base G (2/4);  total (44/49)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

GCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTAACG  >  NC_000913/3180547‑3180841
                                                                                                                                                   |                                                                                                                                                   
gCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACg                                                                                                                                                   <  2:1291399/150‑1 (MQ=255)
gCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACg                                                                                                                                                   <  2:1538785/150‑1 (MQ=255)
 cccAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGc                                                                                                                                                  >  1:292232/1‑150 (MQ=255)
   cAGGTTGTAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGTCAGCCATGACTGGCGCAc                                                                                                                                                <  2:1310031/150‑1 (MQ=255)
   cAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCAc                                                                                                                                                <  2:962272/150‑1 (MQ=255)
   cAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCAc                                                                                                                                                <  2:792887/150‑1 (MQ=255)
      gTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAACACACACCTGGAGTTTTGTGCGTCAGCCATGACTGGCGCAcaaa                                                                                                                                             <  2:135078/150‑4 (MQ=255)
      gTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAg                                                                                                                                             >  2:1218168/1‑150 (MQ=255)
        tGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGcc                                                                                                                                           <  2:204897/150‑1 (MQ=255)
           aaaGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGGTGCGTCAGCCAtg                                                                                                                                        <  1:873638/150‑2 (MQ=255)
             aGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGc                                                                                                                                      <  1:401242/150‑1 (MQ=255)
               cTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCtt                                                                                                                                    <  2:82616/150‑1 (MQ=255)
                tGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCttt                                                                                                                                   >  2:11264/1‑150 (MQ=255)
                 ggAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTa                                                                                                                                  >  2:768020/1‑150 (MQ=255)
                 ggAAGAAGTCACCGCCGAACTGCACCAGATGTGACTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTa                                                                                                                                  >  2:275883/1‑150 (MQ=255)
                       agTCACCGGCTAACTTCACCATATGTGGCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATTAGCTGATGACCAAATTCCGCATTCAAAAACACACATTTAGTTTTGTGCGCCAGTCATGGCTGACGAACCAGCCATCTCTTTATTcgcg                                                                                                                            <  2:327779/150‑1 (MQ=255)
                       agTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTcgcg                                                                                                                            <  1:1332735/150‑1 (MQ=255)
                         tCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTc                                                                                                                          <  2:1114759/150‑1 (MQ=255)
                            ccgccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAACACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTg                                                                                                                       >  2:1524480/1‑150 (MQ=255)
                            ccgccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTg                                                                                                                       >  2:448807/1‑150 (MQ=255)
                            ccgccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTg                                                                                                                       >  2:122988/1‑150 (MQ=255)
                            ccgccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTg                                                                                                                       <  2:608036/150‑1 (MQ=255)
                              gccgAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGAt                                                                                                                     <  1:97745/150‑1 (MQ=255)
                                      gCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTg                                                                                                             >  1:167924/1‑150 (MQ=255)
                                       cACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTgg                                                                                                            >  2:604877/1‑150 (MQ=255)
                                        aCCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTggg                                                                                                           <  2:535147/150‑1 (MQ=255)
                                        aCCAGATGTGCCAGAAAGTGGTGGAAAAACTGATCGCAAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCATTATTCGAGACTAGATCTGGCGTggg                                                                                                           >  2:1237311/1‑150 (MQ=255)
                                         ccAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACCCCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACAAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGa                                                                                                          >  2:896211/1‑150 (MQ=255)
                                          cAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGAt                                                                                                         <  1:161228/150‑1 (MQ=255)
                                            gATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATgg                                                                                                       <  2:1143455/150‑1 (MQ=255)
                                                 gCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTg                                                                                                  <  2:447589/150‑1 (MQ=255)
                                                   cTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGt                                                                                                >  2:986555/1‑150 (MQ=255)
                                                   cTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCGGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGt                                                                                                <  2:1405160/150‑1 (MQ=255)
                                                          tggtggAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTa                                                                                         <  1:939112/150‑1 (MQ=255)
                                                          tggtggAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTa                                                                                         <  2:1332964/150‑1 (MQ=255)
                                                             tggAAAAAGTGATCGCCAGCGATGAGCTGATCACCAAATTCCGAATTCCAAAACACACCTGGAGTTTTCTGCGCCAGTCATTGCTGACGCACAAGCCATCGCTTTATACGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTacac                                                                                      >  2:766968/1‑150 (MQ=255)
                                                              ggAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACt                                                                                     <  2:961647/150‑1 (MQ=255)
                                                               gAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACtt                                                                                    >  1:1297527/1‑150 (MQ=255)
                                                               gAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACtt                                                                                    >  1:1401890/1‑150 (MQ=255)
                                                                aaaaaGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTc                                                                                   >  2:1600290/1‑150 (MQ=255)
                                                                 aaaaGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCt                                                                                  <  2:167924/150‑1 (MQ=255)
                                                                  aaaGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTg                                                                                 <  2:1323954/150‑1 (MQ=255)
                                                                   aaGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTgg                                                                                <  2:1243491/150‑1 (MQ=255)
                                                                   aaGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTgg                                                                                >  1:438840/1‑150 (MQ=255)
                                                                     gTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGGACCTAAACTTCTGGaa                                                                              >  2:1081456/1‑150 (MQ=255)
                                                                      tGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGaaa                                                                             <  1:1011704/150‑1 (MQ=255)
                                                                         tCGCCAGCGATGAGCTGATGACCAAATTCCGCTTTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATa                                                                          >  2:130605/1‑150 (MQ=255)
                                                                         tCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATa                                                                          <  1:276872/150‑1 (MQ=255)
                                                                           gCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAAc                                                                        <  1:1437960/150‑1 (MQ=255)
                                                                            ccAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACg                                                                       >  2:248433/1‑150 (MQ=255)
                                                                              aGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGcc                                                                     <  2:438840/150‑1 (MQ=255)
                                                                               gCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCg                                                                    <  1:1432746/150‑1 (MQ=255)
                                                                                 gATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATTTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGAt                                                                  >  2:1005510/1‑150 (MQ=255)
                                                                                       cTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCa                                                            >  2:1228358/1‑150 (MQ=255)
                                                                                         gatgaCCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAAc                                                          >  2:1556024/1‑150 (MQ=255)
                                                                                             aCCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCa                                                      >  1:746789/1‑150 (MQ=255)
                                                                                              ccAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCAc                                                     >  2:786718/1‑150 (MQ=255)
                                                                                              ccAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCAc                                                     >  2:471371/1‑150 (MQ=255)
                                                                                              ccAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCAc                                                     <  1:1523290/150‑1 (MQ=255)
                                                                                               cAAATTCCGAATTCCAAAACACACCTTGAGTTTTGTGCGCCAGTCATGTCTGACGCACCAGCCATAGCTTTATTCTCGTCTTGATCTGGCGTGTGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACt                                                    <  1:1134213/150‑1 (MQ=255)
                                                                                                      cGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAg                                             >  1:844710/1‑150 (MQ=255)
                                                                                                         aTTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGaa                                                                              >  2:1588822/1‑114 (MQ=255)
                                                                                                         aTTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGaa                                                                              <  1:1588822/114‑1 (MQ=255)
                                                                                                          ttCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAggcgg                                         >  2:264743/1‑150 (MQ=255)
                                                                                                          ttCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAggcgg                                         >  2:4037/1‑150 (MQ=255)
                                                                                                           tCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAggcggc                                        >  2:734693/1‑150 (MQ=255)
                                                                                                            ccAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGgcggcg                                       >  2:1614152/1‑150 (MQ=255)
                                                                                                            ccAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGgcggcg                                       >  1:1352137/1‑150 (MQ=255)
                                                                                                             cAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGt                                      <  2:316991/150‑1 (MQ=255)
                                                                                                              aaaaCACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGtt                                     >  2:973803/1‑150 (MQ=255)
                                                                                                               aaaCACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTc                                    <  1:1218168/150‑1 (MQ=255)
                                                                                                                aaCACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCt                                   <  2:1401890/150‑1 (MQ=255)
                                                                                                                aaCACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCt                                   >  1:226221/1‑150 (MQ=255)
                                                                                                                       cTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTg                            <  1:527587/150‑1 (MQ=255)
                                                                                                                       cTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTATTGATCTGGCGTGTGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGTCGGCGTTCTTTCAGTg                            <  2:762136/150‑1 (MQ=255)
                                                                                                                       cTGGAGTTTTGTGCGCCAGTCATGGCTGACACACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCACCGTCACTATACGAGGCGGCGTTCTTTCAGTg                            >  2:1115096/1‑150 (MQ=255)
                                                                                                                         ggAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTAGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATAAGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGa                          >  2:472773/1‑150 (MQ=255)
                                                                                                                           aGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATc                        <  1:448807/150‑1 (MQ=255)
                                                                                                                            gTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATGTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCATCGTCACTATACGAGGCGGCGTTCTTTCAGTGGATct                       <  1:1237311/150‑1 (MQ=255)
                                                                                                                            gTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATct                       <  1:248433/150‑1 (MQ=255)
                                                                                                                               tggtgCGTCAGCCATGACTGGCGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATctggc                    >  1:1203107/3‑150 (MQ=255)
                                                                                                                                tgtgCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATctggct                   <  1:604877/150‑1 (MQ=255)
                                                                                                                                 gtgCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATctggctg                  >  1:1252979/1‑150 (MQ=255)
                                                                                                                                   gcgtCAGCCATGACTGGCGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGa                >  2:1208898/1‑150 (MQ=255)
                                                                                                                                     gcCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAg              <  1:1600290/150‑1 (MQ=255)
                                                                                                                                          tCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAg         >  2:1139025/1‑150 (MQ=255)
                                                                                                                                          tCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAg         >  2:1144703/1‑150 (MQ=255)
                                                                                                                                           cATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGc        <  1:275883/150‑1 (MQ=255)
                                                                                                                                           cATGACTGGCGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGc        <  1:694299/150‑1 (MQ=255)
                                                                                                                                              ggCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTa     >  2:416946/1‑150 (MQ=255)
                                                                                                                                              ggCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTa     <  1:1228358/150‑1 (MQ=255)
                                                                                                                                                cTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTAAc   <  2:237699/150‑1 (MQ=255)
                                                                                                                                                 tGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTAACg  <  1:1498680/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCCCAGGTTGAAAAGCTGGAAGAAGTCACCGCCGAACTGCACCAGATGTGCCTGAAAGTGGTGGAAAAAGTGATCGCCAGCGATGAGCTGATGACCAAATTCCGCATTCCAAAACACACCTGGAGTTTTGTGCGCCAGTCATGGCTGACGCACCAGCCATCGCTTTATTCGCGTCTTGATCTGGCGTGGGATGGCACTGGTGAACCTAAACTTCTGGAAAATAACGCCGATACGCCAACGTCACTATACGAGGCGGCGTTCTTTCAGTGGATCTGGCTGGAAGATCAGCTTAACG  >  NC_000913/3180547‑3180841

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