Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,692,300 G→C 7.2% L1082V (CTG→GTG)  purL ← phosphoribosylformylglycinamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,692,3000GC 7.2% 272.3 / 10.2 83L1082V (CTG→GTG) purLphosphoribosylformylglycinamide synthetase
Reads supporting (aligned to +/- strand):  ref base G (40/37);  new base C (3/3);  total (43/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

CGCCAGCGTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGTGGCAATATACG  >  NC_000913/2692151‑2692449
                                                                                                                                                     |                                                                                                                                                     
cgccagcgTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAg                                                                                                                                                       <  1:1282787/150‑1 (MQ=255)
   cagcgTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCCTGGAAGATTTCCACGcc                                                                                                                                                    <  1:1542174/150‑1 (MQ=255)
      cgTTTGCGGACGGTGGAAGAAGGTTGCAAACTCGTCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGt                                                                                                                                                 >  2:827871/1‑150 (MQ=255)
      cgTTTGCGGACGGTGGAAGAAGGGTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCCTGGAAGATTTCCACGccctg                                                                                                                                                 <  2:155517/150‑3 (MQ=255)
       gTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCCTGGAAGATTTCCACGccctgg                                                                                                                                                <  1:1592519/150‑4 (MQ=255)
          tGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGt                                                                                                                                             >  2:347118/1‑150 (MQ=255)
            cGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTcc                                                                                                                                           <  2:717539/150‑1 (MQ=255)
                cGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTc                                                                                                                                       <  1:31567/150‑1 (MQ=255)
                   tGGAAGAAGGTTGCAAACTCATCGCGTACACTGTCATTGAACAGGATTGACTTCGACCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAGATCTTCCAGGCCCGTGCGTCCGGTcagc                                                                                                                                    <  2:629080/150‑1 (MQ=255)
                   tGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTcagc                                                                                                                                    <  1:1082121/150‑1 (MQ=255)
                    ggAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTcagca                                                                                                                                   <  1:842485/150‑1 (MQ=255)
                         aagGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCa                                                                                                                              >  2:803493/1‑150 (MQ=255)
                           ggTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACt                                                                                                                            >  2:1468933/1‑150 (MQ=255)
                            ggtccAAAATCATCGCCTACACGGTAATTGAAGATGATTCACTTCGACCAAACTTCAACGGACACCAGCACATACCCGGAGGAGCACCCACCGAACGCGTCCAGGGCGTGGAAATCTTCCAGGCCCGTGAGGGCGGACAGCACGTCAcac                                                                                                                           >  2:727680/5‑148 (MQ=255)
                            gTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTc                                                                                                                           >  2:64148/1‑150 (MQ=255)
                              tGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCat                                                                                                                         >  1:546009/1‑150 (MQ=255)
                                cAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCCGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCatat                                                                                                                       >  2:1196535/1‑150 (MQ=255)
                                 aaaCTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATg                                                                                                                      <  2:252755/150‑1 (MQ=255)
                                 aaaCTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATg                                                                                                                      <  2:503865/150‑1 (MQ=255)
                                  aaCTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCCTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGc                                                                                                                     >  2:578685/1‑150 (MQ=255)
                                       aTCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTc                                                                                                                >  1:650939/1‑150 (MQ=255)
                                        tCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCg                                                                                                               <  1:512824/150‑1 (MQ=255)
                                         cgcgTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGa                                                                                                              >  1:1186419/1‑150 (MQ=255)
                                         cgcgTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGa                                                                                                              >  2:732800/1‑150 (MQ=255)
                                           cgTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGAtg                                                                                                            <  1:1104975/150‑2 (MQ=255)
                                           cgTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATa                                                                                                            >  1:1611759/1‑150 (MQ=255)
                                              acacGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCa                                                                                                         <  1:821429/150‑1 (MQ=255)
                                              acacGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCa                                                                                                         <  1:1160709/150‑1 (MQ=255)
                                               cacGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCAt                                                                                                        <  2:1186419/150‑1 (MQ=255)
                                                 cGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAgg                                                                                                                                 <  1:231079/123‑1 (MQ=255)
                                                 cGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAgg                                                                                                                                 >  2:231079/1‑123 (MQ=255)
                                                 cGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCa                                                                                                      >  1:581136/1‑150 (MQ=255)
                                                     cATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGc                                                                                                  >  1:417976/1‑150 (MQ=255)
                                                         gAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGc                                                                                              >  2:1277434/1‑150 (MQ=255)
                                                         gAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGc                                                                                              >  2:370923/1‑150 (MQ=255)
                                                         gAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGc                                                                                              >  1:5259/1‑150 (MQ=255)
                                                          aaCAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCa                                                                                             <  1:864896/150‑1 (MQ=255)
                                                          aaCAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCa                                                                                             >  2:306210/1‑150 (MQ=255)
                                                             aGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACgg                                                                                          <  2:1529270/150‑1 (MQ=255)
                                                             aGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACgg                                                                                          <  2:1523613/150‑1 (MQ=255)
                                                              ggATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCATGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAACCCTGCACggg                                                                                         >  1:1555203/1‑149 (MQ=255)
                                                                aTTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTgg                                                                                       <  1:829109/150‑1 (MQ=255)
                                                                       tCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTg                                                                                >  2:991805/1‑150 (MQ=255)
                                                                         gCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTgcc                                                                              >  1:396332/1‑150 (MQ=255)
                                                                           ccAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTgccgc                                                                            >  1:554129/1‑150 (MQ=255)
                                                                                cTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCAttt                                                                       <  1:732800/150‑1 (MQ=255)
                                                                                 ttCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTc                                                                      <  2:241256/150‑1 (MQ=255)
                                                                                    aCCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAAc                                                                   >  2:207476/1‑150 (MQ=255)
                                                                                      cGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAACGCCTGCACGGTGGAAAGCTGCCGCCATTTCACCAt                                                                 >  2:1142421/1‑150 (MQ=255)
                                                                                           cccAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAg                                                            >  2:1102249/1‑150 (MQ=255)
                                                                                           cccAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAg                                                            >  1:1134178/1‑150 (MQ=255)
                                                                                           cccAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAg                                                            <  1:1034367/150‑1 (MQ=255)
                                                                                           cccAGCACATCACCCTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAg                                                            <  1:578685/150‑1 (MQ=255)
                                                                                            ccAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGt                                                           <  1:803493/150‑1 (MQ=255)
                                                                                                 acaTCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACg                                                      <  1:891765/150‑1 (MQ=255)
                                                                                                       ccGTAGGAGAACCCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGGGCGTCCGGTCAGCAGGTCACTCATATGCCCGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGc                                                >  2:998034/1‑150 (MQ=255)
                                                                                                          tAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCa                                             <  2:1526294/150‑1 (MQ=255)
                                                                                                          tAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCa                                             <  1:727680/150‑1 (MQ=255)
                                                                                                           aGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCAc                                            <  1:778208/150‑1 (MQ=255)
                                                                                                            ggAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACg                                           >  2:743909/1‑150 (MQ=255)
                                                                                                            ggAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACg                                           >  1:807915/1‑150 (MQ=255)
                                                                                                             gagaAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCATGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGc                                          <  1:1388929/150‑1 (MQ=255)
                                                                                                              agaAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCa                                         <  1:548600/150‑1 (MQ=255)
                                                                                                              agaAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCa                                         >  2:504756/1‑150 (MQ=255)
                                                                                                                 aaCCCCCGCACGCGACCAGGGCATTGAAATCTTCCAGGCCCGTGCGTCGGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGGGGAAAGCTGCCGCCATTTCAACATGCGAGTTCCCGCCCTGCTCACGCAGtt                                      <  2:1555203/150‑2 (MQ=255)
                                                                                                                    caccGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACgg                                   >  1:1266682/1‑150 (MQ=255)
                                                                                                                     accGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGc                                  <  1:743909/150‑1 (MQ=255)
                                                                                                                            gcgACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTc                           <  1:60610/150‑1 (MQ=255)
                                                                                                                             cgACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACg                                                      >  2:592085/1‑122 (MQ=255)
                                                                                                                             cgACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACg                                                      <  1:592085/122‑1 (MQ=255)
                                                                                                                             cgACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATACCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCg                          >  2:682401/1‑150 (MQ=255)
                                                                                                                              cgcacgggCCTGGAAGATTTCCACGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCgg                         >  2:1385809/6‑150 (MQ=255)
                                                                                                                                cacgggCCTGGAAGATTTCCACGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGAc                       >  1:1575861/4‑150 (MQ=255)
                                                                                                                                 cAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACg                      <  1:133937/150‑1 (MQ=255)
                                                                                                                                 acgggCCTGGAAGATTTCCACGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACg                      >  1:434057/3‑150 (MQ=255)
                                                                                                                                  aGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGt                     >  2:119307/1‑150 (MQ=255)
                                                                                                                                       gTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGcc                <  1:64148/150‑1 (MQ=255)
                                                                                                                                          gAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGt             >  2:1190349/1‑150 (MQ=255)
                                                                                                                                               cTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGTGGCaa        >  2:1421666/1‑150 (MQ=255)
                                                                                                                                               cTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGTGGCaa        <  2:554129/150‑1 (MQ=255)
                                                                                                                                                ttCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGTg            >  2:1049822/1‑145 (MQ=255)
                                                                                                                                                ttCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGTg            <  1:1049822/145‑1 (MQ=255)
                                                                                                                                                ttCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCAGCCATTTCAACATGAGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTAGGAAGTGCGCCAGTGACAat       >  2:1335477/1‑150 (MQ=255)
                                                                                                                                                 tCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGTGGCAata      <  2:1611759/150‑1 (MQ=255)
                                                                                                                                                     ggCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGTGGCAATATACg  >  2:570130/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCCAGCGTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGTGGCAATATACG  >  NC_000913/2692151‑2692449

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